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CORRECTIONS 


Volume 35, 1957 
Page 572. equation (6), sin should read 


sin x”. 


Volume 38, 1960 
Page 332. line “(Wilson, Rose, and Pomerantz should read ‘‘(Wilson, 
Rose, and Pomerantz 


Page 333. line G., Rose, C., and Pomerantz, 1959. Can. 
Phys. 38, 328.” should read G., Rose, C., and Pomerantz, 1960. 
Can. Phys. 38, 328.” 
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SPECTROSCOPIC TEMPERATURE MEASUREMENTS 
SHOCK TUBE USING THERMOMETRIC MOLECULE! 


ABSTRACT 


Rotational intensity measurements the spectrum, excited through 
driven shock waves argon have been used infer 
the gas between 6350° and 8750° The measured values agree well with 
gas kinetic temperatures inferred from simple gas dynamic theory and shock- 
wave velocity measurements. 


INTRODUCTION 


Temperature parameter great importance gas dynamics, particu- 
larly the study transient luminous phenomena associated with shock 
waves. therefore became necessary determine, spectroscopically pos- 
sible, the temperature the luminous gas behind the reflected shock wave 
during series experiments upon shock excitation powdered solids 
(Nicholls 1959; Parkinson and Nicholls 1959a, Nicholls and Parkinson 
1957, 1958; Nicholls, Parkinson, and Van der Laan 1959; Nicholls, Watson, 
and Parkinson 1959). 

Spectroscopic temperature measurements using line reversal methods 
(Clouston, Gaydon, and Glass 1958; Clouston, Gaydon, and Hurle 1959) 
have been employed with success some workers for temperatures 
about 3600° the work described below, rotational temperatures (in 
the range radicals introduced into the reflected shock 
wave argon have been measured (Parkinson and Nicholls Nicholls 
1959) and compared with gas kinetic temperatures inferred, means 
one-dimensional theory (Glass 1958), from initial pressures, temperatures, 
and gas composition either side the diaphragm and also from velocity 
measurements. The gas behind the reflected rather than the incident shock 
wave was studied because the limitations the bursting strength the 
high pressure section the tube. 
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METHOD 


well known (Johnson 1949; Herzberg 1950; Gaydon and Wolfhard 
1953) that the intensity distribution lines molecular band 
given 
and further that Boltzmann energy distribution exists the rotational 
energy levels the molecules, i.e. 


then 
(3) = C/Qrot Sy exp (—Eror/RT rot). 


The notation conventional and defined the above references. The 
Boltzmann energy distribution rotational levels implies that thermal 
equilibrium exists between translational and rotational degrees freedom. 

From equation (3) seen that plot log versus for lines 
band should lead straight line slope thermal equilibrium 
obtains. The rotational temperature may thereby measured from 
intensity distributions within bands. The bands should preferably 
simple, well-resolved structure such the violet system. 

attempt introduce into the hot gas the reflected shock 
wave, series survey experiments involving shock excitation powdered 
metallic cyanides and other containing compounds were first performed 
(Parkinson and Nicholls all cases the violet system was easily 
excited, were the Swan bands, and underlying carbon-particle black- 
body continuum which prevented accurate intensity measurement the 
lines. Shock excitation solids appears ablate atoms (Nicholls 1959) which 
recombine into diatomic molecules the environment cool enough. This 
suggests that the carbon-particle continuum and bands can then reduced, 
relative bands, reducing the C:N ratio the gas. This was success- 
fully achieved shock-exciting solid mixture which gave rise hot 
luminous gas containing and with small C:N ratio. After series 
trials with various molar proportions the mixture 
was found behave very satisfactorily producing bands, few 
bands, and carbon-particle continuum. 


EXPERIMENTAL 


The powdered mixture was made into paste with acetone and painted 
the center mylar membrane. After evaporation the acetone, the 
membrane, coated with dry powder, was clamped near the window end 
the shock tube (Parkinson and Nicholls 1959a) with the powder towards the 
window intercept the reflected shock wave. had previously been 
established that the luminosity from shock excitation mylar involved only 
the lines sodium 5893 and the 4227-A line calcium. Spectra 
were recorded upon Kodak Royal Pan film using ARL Raman 
spectrograph with reciprocal dispersion A/mm 4200 and free 
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The response curve the emulsion was obtained using calibrated 
step wedge mounted the spectrograph slit. Homochromatic photometry 
was sufficient for the violet (0,1) band 4216 which only extends 
over about The relative intensities the best resolved lines 
(J’ 45-59 inclusive) the branch the band were measured incident 
Mach numbers 5.7, 5.9, and 6.3. Adequate spectra could not photo- 
graphed below Mach 5.7 with one exposure, band intensities decrease 
significantly with decrease Mach number. 


RESULTS 


Plots log vs. yielded straight lines each case. rotational 
temperature could inferred from the slope each line. The plots and best 
lines are shown respectively Figs. la, 1b, Table 


28 


2800 32 2000 2 


the inferred rotational temperature compared with gas kinetic tempera- 
ture derived from one-dimensional shock wave theory (Glass 1958) and the 
initial conditions pressure temperature and gas composition each side 
the membrane prior initiation the shock wave for number typical 


cases. 


TABLE 
Comparison between measured rotational temperature and calculated 
gas kinetic temperature 


Calculated gas kinetic Measured rotational 
Mach number temperature, temperature, 
6.3 9000 +5% 
7340 
5.9 8200 
6800 
5.7 7650 


q 
q 
116 1,06 \ 
t\ 
dq 
MAC 7 b) MA 5. 0.9 MACH 6. { 
0.98 0.96 
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DISCUSSION 


will noticed Table that the measured rotational temperatures are 
10-15% lower than the ideal calculated temperatures based one-dimen- 
sional theory and the assumption perfect gas. This consistent with the 
results other experiments (Parkinson and Nicholls which shock- 
wave velocities were measured while making time-resolved spectroscopic 
studies MgO. The measured velocities were about lower than the 
calculated velocities, which was equivalent about 10% reduction the 
effective gas kinetic temperature below that predicted one-dimensional 
perfect gas theory. Thus the measured rotational temperature probably 
reasonable measure the actual translational temperature the real gas 
the reflected shock wave. 

well known that most prevalent source error measurement 
rotational temperature self-absorption. This can especially significant 
when the source high pressures 100 p.s.i.a. the gas behind the 
reflected shock wave). The effect self-absorption the semilogarithmic 
plots reduce the slope the straight lines, which gives rise higher 
apparent temperature. Thus the obtaining straight line the conven- 
tional semilogarithmic plots such Fig. does not necessarily imply that 
self-absorption absent (Marr 1957). However, the isointensity method 
spectroscopic temperature measurement Knauss and McCay (1939) treats 
pairs band lines equal intensity and therefore independent self- 
absorption artifacts. This method was used obtain estimates rotational 
temperature for several pairs P-branch lines each case, and the tempera- 
tures obtained agreed well with those given Table derived from the slopes. 
Thus self-absorption did not appear play significant role affecting our 
measurements. 

therefore apparent that, under suitable conditions, the molecule 
good spectroscopic thermometer for indication the transient tem- 
perature the hot luminous gas behind reflected shock waves argon over 
the approximate temperature range 6000° 9000° Cooling the 
luminous gas ablation from the powder particles was probably minimized 
and compensated for the exothermic nature the and 
mixture used. 

While this paper was preparation, recent paper Sobolev al. (1959) 
became available. They describe measure temperatures shock 
waves air spectroscopically by: comparison relative intensities 
lines sodium (and also potassium) excited shock excitation following 
dusting the shock tube walls with the respective chlorides, (b) measurement 
spectral intensity distributions violet bands excited when dicyanogen 
was introduced into the shock-excited air. 

the experiments with alkali metal lines, which also require knowledge 
appropriate oscillator strengths, poor agreement between measured excita- 
tion temperatures and theoretically predicted gas kinetic temperatures was 
attributed self-absorption effects, and they conclude that the method 
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not useful for studying the the air behind shock wave’’. They also 
point out that agreement between their spectroscopic rotational temperatures 
from line intensities and gas kinetic temperatures inferred from velocity 
measurements was not good. They attribute this the effects finite 
vibrational relaxation times required dissociate into 2CN and 
equilibrate with the hot gas. 

Evidently, our experiments, the exothermic nature 
together with the most likely formation from atomic gas (Parkinson 
and Nicholls 1959; Nicholls, Watson, and Parkinson 1959) make the experi- 
mental environment more favorable for spectroscopic temperature measure- 
ment the hot argon with which the was mixed. 


REFERENCES 


G., Gaypon, G., and 1958. Proc. Roy. Soc. 248, 429. 

G., Gaypon, G., and 1959. Proc. Roy. Soc. 252, 143. 

and 1953. Flames, their structure, radiation and tem- 
perature (Chapman Hall, Ltd., London). 

1958. University Toronto Institute Aerophysics Review, No. 12, Part 

1950. Spectra diatomic molecules (D. Van Nostrand Co., Inc. New 

ork 

OHNSON, R.C. 1949. introduction molecular spectra (Methuen Co., Ltd., London). 

and McCay, 1939. Phys. Rev. 52, 1147. 

Marr, 1957. Can. Phys. 35, 1270. 

1959. Roy. Astron. Soc. Can. 53, 109. 

NICHOLLS, and 1957. Chem. Phys. 26, 423. 

1958. Can. Phys. 36, 625. 

W., H., and VAN DER Laan, 1959. Appl. Phys. 30, 
797. 

an. 53, 

Report No. Contract 19(604)-4560, Department Physics, University 
Western Ontario. 

1959b. Shock tube spectroscopy Scientific Report No. Contract 19(604)- 

4560, Department Physics, University Western Ontario. 

and ISAEF, 1959. Optika Spektroskopiya (Optics and 
Spectroscopy, Optical Society America, translation 185). 


q 
q 
4 


COMPTON-ELECTRON WITH 
TWO-DIRECTIONAL FOCUSING! 


LEE-WHITING 


ABSTRACT 


Improvements the design one type Compton-electron spectrometer 
for y-rays are proposed. The design requires magnetic field cylindrical 
symmetry and slow radial variation, simply curved radiator, and system 
apertures. Electrons are accepted only they are ejected from the radiator with 
small components momentum two orthogonal directions perpendicular the 
incident y-ray. Since the magnetic field can also used measure the momen- 
tum the selected electrons, the instrument can function y-ray spectro- 
meter. Higher-order aberrations are discussed, and method calculating the 
values the various spectrometer parameters corresponding maximum 
efficiency given. Calculations the intrinsic line-width, caused the motion 
the electron within the atom before collision with the photon, are carried out. 


INTRODUCTION 


Consider parallel beam finite cross section monoenergetic y-rays 
incident upon thin sheet some light material such cellophane. Measure- 
ment both the energy electron which has been set motion 
Compton scattering and its angle ejection with respect the y-beam 
enables one calculate the energy the incident photon. apparatus 
which selects the electrons recoiling particular direction and measures 
their energies will thus serve y-ray spectrometer. the sort device 
considered this paper the electrons ejected the forward direction 
are utilized; their momenta are measured the type 
which the electron motion predominantly perpendicular the magnetic 
field, i.e. the type Siegbahn (1955, Chap. 

example the class instrument discussed, let take the 
spectrometer described Groshev al. (1956). well-collimated beam 
neutron-capture y-rays from extended source reactor incident upon 
curved sheet polystyrene (to called the radiator) marked Fig. 
The magnetic field the 6-spectrometer axially symmetric about per- 
pendicular raised and has mirror symmetry the plane the figure; 
the field this plane perpendicular it. Within the circle radius the 
magnetic field uniform. The radiator, which lies entirely the uniform 
field, shaped that forward Compton electrons produced photons 
the correct energy will move circular orbits which pass through the focus 
after turning through approximately 90°. After passing the electrons, 
now outside the circle radius are region where the field falls off 
radially such manner produce improved radial focus 180° 
beyond Though Groshev al. not say so, the field presumably 
shaped accordance with the theory Beiduk and Konopinski (1948); 
this case there also small amount axial focusing outside the circle 
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Fic. the principle the type instrument used Groshev al. 
(1956). The thin radiator, lies the uniform field region within the circle radius ro. 
The orbits electrons ejected the extreme rays the well-collimated y-beam are labelled 
with one two arrow-heads. 


radius ro. should clear that the combination curved radiator and two 
slits selects those electrons the desired energy having initially zero com- 
ponent momentum the direction perpendicular the y-beam and 
the plane the figure. The only limitation placed upon the component 
initial momentum perpendicular this plane results from the finite heights 
the radiator and the slits and order reduce the back- 
ground, coincidences are recorded between the counters and placed 
behind the slits and respectively. 

The principle the extended radiator used Groshev al. was first 
suggested Dzhelepov and Orbeli (1948), and was employed the design 
the (see Dzhelepov 1954). discussion the use curved 
radiator uniform magnetic field has been given Dzhelepov (1951). 

Great improvement was effected the performance 
the type when was realized that radially varying magnetic fields 
could used and that the focus need not 180° (see Siegbahn (1955, 
Chap. III) for discussion and references). The improvement resulted from 
first-order focus the axial direction, allowing the use shorter detectors, 
and from higher-order focusing the radial direction, producing sharper 
line given transmission. Can one obtain similar benefits with spectro- 
meter the general type illustrated Fig. permitting the use 
inhomogeneous magnetic field and allowing the designer choose the 
relative positions the radiator and slits? turns out, shall show 
Section that with this extra freedom one can arrange select electrons with 
zero initial momentum the axial direction; the length detector needed 
also decreased. also possible employ higher-order radial focusing 
sharpen the spectral line the instrument. 
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Before going into the detailed discussion the various components the 
have some estimate the intrinsic width caused the motion the 
electrons within the target atoms before the scattering takes place. Section 
calculations this effect are presented; they are more complete than those 
which have already been given Dzhelepov (1955). few observations 
the shape the pole pieces iron-cored magnet designed produce the 
required field are made Section 

discussing the Compton process shall take formulae from the Chapter 
Davisson Seigbahn’s handbook (1955). Energies will measured 
units momenta units moc, being the electron’s rest mass. Let 
the energy the photon and let the energy electron 
recoiling the forward direction. Then 


(1) 


Now since the relative half-width the momentum line 8-spectrometer 
independent the momentum, usually more convenient work with 
momentum instead energy. the momentum forward recoiling 
electron, then 


(2) 


will also useful have expression for the cross section per target 
electron for the ejection Compton electron per unit solid angle the 
forward direction, viz. 


here cm? the square the classical radius the electron. 
Note that although the total cross section for Compton scattering decreases 
the photon energy increased, the differential cross section the forward 
direction increases sharply. increased from 20, rises from 
2.97 Thus, for energies which the half-angle 
the cone electrons accepted Compton spectrometer small compared 
the width the forward lobe the differential cross section, the efficiency 
the instrument will improve the photon energy increased. 

Although customary quote the relative half-width the line profile 
plotted against electron momentum measure the resolution Compton 
spectrometer, the quantity which relevant the resolving y-rays the 
relative half-width the line profile plotted against the corresponding 
photon energy. differentiation (2) obtain 


(4a) Aa/a U(a)Ap/p 


The function falls monotonically from 1.2 increased from 
Hence may use the electron momentum resolution instead the 
photon energy resolution and not seriously misled. 


LEE-WHITING: COMPTON-ELECTRON 7-SPECTROMETER 723 


PRINCIPLES FOCUSING 


The general arrangement the spectrometer under consideration this 
section, shown Fig. very similar that Fig. The magnetic field 


Fic, the instrument discussed this paper. The radiator lies astride the 
optic circle. The magnetic field varies radially all regions. 


axially symmetric about the perpendicular erected and has mirror 
symmetry the plane the page. The two slits and are situated 
circle (to called the optic circle) radius ro, centered and lying 
the plane symmetry. describing the position electron con- 
venient use the radial and axial departures from the optic circle, expressed 
units 


The third co-ordinate will the angle cylindrical co-ordinates, 
measured clockwise sense from the line OR. 

well known, the properties the magnetic field the neighborhood 
the optic circle are completely determined the radial variation the 
magnetic field the plane symmetry; thus the coefficients the ex- 
pansion 


n=0 


specify the field. full discussion the expansion parameters the field 
components and method calculating the electron orbits given 
Lee-Whiting and Taylor (1957). 

Fig. well-collimated beam y-rays incident upon radiator 
placed roughly symmetrically astride the optic circle. The beam forward- 
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second focus for The symmetric arrangement the radiator 
gives the minimum aberration (and hence the sharpest line) for given 
width radiator. discussing the general principles focusing, 
this section, permissible imagine perfectly parallel y-beams and 
slits zero width. The effects finite slit-widths resolution are discussed 
Section 

The first problem show how design the radiator surface ensure 
that particles the desired momentum pass through the slit Since 
would extremely difficult construct rigid warped surface out thin 
cellophane, the radiator assumed portion the right cylinder 
erected the curve the plane symmetry. The calculation performed 
only for those photons incident the plane symmetry. Estimates the 
aberrations caused this approximation are made Section 

Let assume that f(@) describes orbit passing through and 
let (see Fig. point this orbit. the angle between the tangent 


Fic. Illustration the angles involved the derivation the equation the radiator 
surface, RT; general point the orbit 


the orbit and the line PO, the condition that the tangent parallel 
the incident y-direction 
(7) 
Now from elementary geometry have 
tan 


Thus the co-ordinate the point common the orbit and the 
radiator may found solving the equation 


(8) 


\ 
4 \ } 
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the co-ordinate may obtained, course, from 


Now let consider one-parameter family orbits passing through for 
the parameter convenient take the maximum departure the orbit 
question from the optic circle, Then, replacing obtain 
parametrization the required radiator curve from the equations (9). 

good first approximation the radiator shape may obtained using 
the first-order solutions the orbit equation. adapting Lee-Whiting 
(1957) find 


(see Fig. for the definition The resulting radiator curve easily 
shown 


practical cases the arc the curve defined (11) nearly straight 


line. The position this line can specified giving the value the 
angle, between the tangent the curve and the radius OR. 


Calculations the radiator shape can, course, carried out using orbit 
equations higher order approximation, though formula simple 
(11) should not expected. The second-order term which should 
added the (10a) can calculated the methods used Lee-Whiting 
(1957); 


Calculation the radiator shape for line source considered briefly 
Appendix. 

Having got the electrons far the slit may now think the 
second focus The condition for the first-order radial focus 
electrons emitted point source well known; the present notation 


The two radial foci have been arranged the cost only one constraint, 
equation (13), the three disposable parameters and may 
thus approach the question focusing the axial direction with the 
edge that there still considerable freedom available. the spectro- 
meter particles emitted from one point the source different angles 
the plane symmetry pass through one point (in the first approximation) 
the exit slit. This the opposite the behavior required Compton 


3 
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spectrometer. Here wish select only those orbits which are initially 
parallel the plane symmetry. the axial variation the spectro- 
meter described that required for the Compton spectrometer 
Note also that the selection orbits should independent 
the axial co-ordinate the point ejection the radiator. 

For the moment shall treat only those points the radiator lying upon 
the generator erected upon the optic circle; the defocusing for other radiator 
points will treated aberration Section the first order the orbits 
under consideration lie the cylinder containing the optic circle. From 
Lee-Whiting (1957) one obtains 


for this case; the value the point ejection, proportional 
the initial slope the orbit. bring about the desired selection orbits 
one should satisfy the equation 


(15) => 


and place aperture infinitesimal height, well width, Then 
only those orbits for which zero pass through the final aperture, and, 
furthermore, the value does not affect the selection. Sketches the 
development the cylinder containing the orbits are given Fig, 


T<O 
tt 


2 


Fic. Development showing orbits originating from points the central generator 
the radiator. Electrons ejected horizontally different heights, i.e. and pass through 
not pass through 


Remembering that both and are functions alone, see that 
the desired two-directional focusing imposes two constraints, (13) and (15), 
the three disposable parameters and ay. One can find values 
and satisfying (13) and (15) for range values from the equations 


and 
(17) 
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and one gets the greatest radiator width for given maximum radial depar- 


ture from the optic circle. 


When decreased below 1/10, approaches 


and the detector too near the y-beam. There can little point using 
values less than 1/10. Curves and well for the useful 
range are given Fig. 


SCALE FOR AND ,DEGREES 


SCALE FOR 8, , DEGREES 


focusing. 


versus for Compton spectrometer with two-directional 


For the 8-spectrometer the two-directional first-order focusing requires 


the fixing the paramete 


used eliminating second-order aberrations. the present case, however, 
approach the question second-order focusing with two free parameters, 


and 
Formulae for the import 
the work Lee-Whiting 


ant aberrations may adapted from 


(1957). Since that paper the symbols and 


refer the value and slope the orbit the source the 


(18a) 
and 


fication, one finds: 


(20a) 


them and respectively, where 


sin 


Thus, after some simpli- 
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(22) 

(23a) 


The terms depending upon the width the slit have been omitted, since 
this width should small enough importance the first order only. 

Our next task try choose the four (19) such way 
minimize the total second-order radial aberration. There point elimina- 
ting the term since, shall see Section larger term proportional 
arises from the dependence the electron energy angle ejection. 
possible choose value for any value the range —0.2 
the specified range, however, can one make more than one the three 
coefficients vanish. one arranges have equal zero, then (since 
occurs only the first term), larger value and hence wider radiator, 
can used. principle one could continue cancel the coefficients 
higher and higher powers increasing the permissible width the 
radiator each stage. theoretical limit the degree useful cancellation 
arises when terms such and HtT become appreciable. Let 
assume that the first powers have zero coefficients; 
small enough that the power series rapidly convergent, then the only 
H-dependent term that must retained CH". 

the case the was possible, arranging cancel 
either two coefficients, have wide-aperture high-aperture instru- 
ment (Lee-Whiting and Taylor 1957). For the Compton spectrometer, 
have just seen, the wide-radiator type instrument corresponds the 
choice the high-radiator type, which would require the simultaneous 
vanishing both and not possible. 

Returning the wide-radiator spectrometer, see that the condition 
that should zero enables calculate C3, and for given 
value a;. Curves the numerical values calculated for the range 
—0.1 are given Fig. 

Since has been decided design the spectrometer such way that 
the coefficients the second-order term zero, knowledge the 
H-dependent third-order aberrations useful. The most important these 
terms, has been evaluated: 


Though might not possible design the pole pieces produce desired 
value still useful able estimate the size this aberration. 


2 
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-.2 -18 -16 -.10 


Fic. a2, C2, C3, and against for spectrometer. 


the Compton-electron spectrometers which have been built study 
y-rays the resolution achieved (2.5% 5%) has not been good enough 
make the failure restrict the axial momentum the electrons selected 
important source line broadening. Compton spectrometers are 
compete with other types instrument, however, improvement the 
resolution order magnitude desirable. such resolutions the 
improved form focusing becomes almost necessity. 

For resolutions the order 0.3% will not always possible use 
counter behind the first slit, because the energy loss and change direction 
caused multiple scattering the counter windows. counter could prob- 
ably constructed give tolerable spread for the electrons ejected 
photons energy Mev, though the deterioration the axial focusing caused 
the change electron direction motion might serious. Because the 
amount multiple scattering decreases rapidly the energy increased, 
there should difficulty photon energies greater than Mev. When 
the two-counter mode operation cannot used, however, one the great 
advantages the Compton spectrometer, its low background counting rate, 
danger being lost. alternative method minimizing the back- 
ground which might effective use counter telescope behind the final 
aperture the system. Since, the first-order theory orbits, the desired 
electrons are focused point the exit aperture the spectrometer 
are discussing, the size this aperture will relatively small. The smaller 
exit aperture admits fewer stray electrons and permits the use smaller 
detector, which would less affected environmental background. The 
background can certainly much decreased recording coincidences be- 
tween the selected electrons and the backscattered photons (Motz 1956; Motz 
al. 1959), but the counting rate will reduced appreciably. 
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INTRINSIC LINE-WIDTH 


the previous sections have assumed that the electrons taking part 
the Compton scattering are initially rest; this case all electrons ejected 
the forward direction have precisely the same value energy momentum. 
When one takes into account the motion the electron within the target 
atom before the collision, the sharp spectral line replaced one finite 
width. This intrinsic width is, course, the minimum width attainable with 
any Compton spectrometer using that particular material the radiator. 

The idea that the atomic motion the electrons could cause broadening 
the energy line the scattered photons was first suggested Jauncey 
(1925), who performed his calculation using the classical Bohr orbits. Wave- 
mechanical calculations this same effect were carried out and compared 
with experiment Mond (1933). Dzhelepov (1955) calculated the width 
the energy line forward recoiling electrons for two values the energy 
the incident photon. The purpose the present computation extend 
the range photon energies, and examine the effect using better atomic 
wave functions. 

For the case electrons scattered the forward direction the conservation 
equations are much simpler than those the general case treated 
Mond. Let work, for the moment, with the electrons one atomic shell, 
and let their binding energy units The conservation energy 
the equation 


(25) 


here and are the photon energies before and after scattering and 
again the kinetic energy the ejected electron. Let the magnitude the 
linear momentum vector the electron before the collision (in units 
moc), and let the angle between this vector and the direction the incident 
photon Then conservation momentum gives 


here the momentum the electron after collision (entirely 
direction the incident y-ray), and the angle scattering the photon. 
The equations (25) and (26) determine the quantities a’, and functions 

The distribution values can calculated with the aid quantum 
mechanics. Following Dzhelepov select nitrogen suitable 
cellulose atom. will assumed that the system consists isolated 
atom; the effects molecular binding should considerably smaller than 
those atomic binding. For each atomic orbital shall assume spatial 
wave function the form 


| 
730 
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The corresponding momentum wave function 


assumed that the unit length the first Bohr radius; the fine- 
structure constant. Now the quantity 


clearly the momentum density appropriate full shell electrons. For 
the incomplete p-shell nitrogen shall use Then the angular 
dependence the momentum density averaged out, i.e. the probability 
photon encountering electron with momentum specified 
made independent shall use the symbol for the total momentum 
density function for all the electrons the nitrogen atom. 

Radial wave functions for the nitrogen atom have been calculated 
Hartree and Hartree (1948) the method the self-consistent field with 
exchange. was convenient fit these wave functions with analytical 
formulae the form suggested Lowdin and Appel (1956), viz. 


i=l 
i=1 
(30c) 
i=1 


set parameters yielding fit (to about 0.1% where the wave functions 
were appreciable) was found with the aid program for Burroughs 205 
computer developed Pope and Kennedy (1959) for another problem. The 
numerical values are given Table 


TABLE 


Parameters for Lowdin fit Hartree wave functions for nitrogen 


ro A, As A; a; as 

21.56 13.41 6.016 8.358 
0.3181 3.191 6.752 1.657 6.465 
1.295 2.411 5.468 


1.195 


When hydrogenic wave functions are used the integration (28) and the 
summation (29) are easily carried out for the states required; indeed, 
formula for the general case has been given Podolsky and Pauling (1929). 
For wave functions the Lowdin the calculation more tedious, but 
still elementary. Fig. the radial momentum distribution 


| 
| 
~ 
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Fic. Radial distributions electron momenta for the nitrogen atom; the abscissa 
momentum units moc. Curve (a) computed for Hartree wave functions and curve for 
unscreened hydrogenic wave functions; curve (c) taken from the paper Dzhelepov (1955). 


calculated from the Hartree nitrogen wave functions given. For the sake 
comparison the same distribution calculated for unscreened hydrogenic 
wave functions also plotted; because the stronger resultant binding forces, 
the latter distribution falls off more slowly for large values than does the 
former. The satisfaction reasonable accuracy the normalization condition 


was used check the calculated distributions. his calculation 
intrinsic line-widths Dzhelepov used the screened hydrogenic wave functions 
for carbon computed Mond (1933) (obtained via Compton and Allison 
1934), and made adjustments for the extra electron nitrogen. The values 
the function obtained are also shown Fig. 

Having seen how compute the initial momentum probability density, 
turn our attention again the solution the conservation equations (25) 
and (26). The dependence upon the small quantities and may 
represented with sufficient accuracy McLaurin expansions far the 
the pair simultaneous equations arising from the differentiation (25) 
and (26a); recall that for Thus 


this section convenient use the symbol for the momentum after 
collision electron initially rest, i.e. for the expression (2). The 
truncated series for may now written 


(35) 
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Because the peak the distribution Fig. occurs about 0.01, 
the error resulting from the omission higher powers should not 
large. 

The magnitude for the electrons nitrogen can calculated from 
the critical absorption wavelength for rays; using value 31.1 for 
this wavelength (Compton and Allison 1934) one gets 0.78 moc? for 
the accuracy required this problem there need distinguish 
between the values for the and electrons. From the value 14.5 
for the first ionization potential nitrogen (Handbook Chemistry and 
Physics 1946) one finds 0.28 moc? representative value for the 
shell. First note that the values for both shells are sufficiently small 
justify the retention only the first power (34). The effect the 
presence shift the spectral line toward smaller values momentum. 
This shift would have effect the line-width were not for the fact that 
the magnitude the shift different for different shells. Compared the 
effect values the order 0.01, the shift completely negligible for 
the shell. For the electrons the shift about 10% the effect the 
important values (34). However, since only two the seven electrons 
are the shell, the effect may ignored without fear serious 
error, i.e. shall put equal zero equation (34). his treatment 
the conservation energy Dzhelepov assumed that the electron was initially 
free; this equivalent replacing (34) negative function Since 
Dzhelepov also retained terms order not easy distinguish the 
effect his different treatment the conservation energy. The combined 
effect his additional terms appears affect the line-width only very slightly; 
the line made asymmetric with the steeper side the lower-momentum 
region. more correct calculation, i.e. one using suitable values for 
the different shells, the steeper side would lie the region higher momentum. 
The calculations this paper should more correct with regard the 
conservation energy than those Dzhelepov, though the difference 
probably negligible. 

Now are ready derive expression for the relative number, 
electrons with momenta after collision the range p+dp. Let 
assume for the present that the probability scattering independent 
the values and Then the number collisions which the range 
having their momenta initially this range, sin wdg dw. 
shall find V(p) transforming one the variables integration 


et 


The integration over will eliminated favor integration over 
sufficiently good approximation the relation between and obtained 
range Let and the maximum and minimum values respectively 
attained for fixed and variable 


a 
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(37) 
Noting that dp/d cos find that 


suitable parameter use describing the line-shape the fractional 


departure the momentum from its nominal value, i.e. 
convenient also define the function 


y 
Upon changing the order integration (38) one finds 


Thus the line shape, the dependence the relative number electrons 

The function F(y) has been evaluated for the Hartree and for the unscreened 
hydrogenic momentum densities numerical integration the Runge-Kutta 
process. The results are shown Fig. also possible, small number 


Q 4 4 


Fic. function F(y) for Hartree, (a), and for unscreened hydrogenic, (b), wave 
functions. This function gives the number electrons against fractional departure from the 
mean momentum for incident photon energy 0.420 


elementary integrations, evaluate F(y) analytically for the hydrogenic 
density: 

fine-structure constant. The part (41) depending upon arises from the 


» 
60} | 
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shell, the part depending upon from the shell; different effective values 
could used for the two shells. Similar, but much more complex, results 
could obtained for the momentum density calculated from the Lowdin 
formulae fitted the Hartree wave functions. Instead the method numerical 
integration was tested comparing its values F(y) for the hydrogenic 
case with those obtained from (41). Another check afforded the normal- 
ization condition inherent (40), viz. 


0 
All checks were satisfied three four significant figures the regions 
importance. 


There are two qualitative features the spectral function which 
should noted. Firstly, the line symmetric about i.e. about the 
momentum value corresponding electron initially rest. Secondly, the 
line-shape the same all photon energies, suitable changes scale 
the abscissa are made. These two properties the line are valid, should 
understood, only because the assumptions that the final electron momen- 
tum linear function initial electron momentum and that the binding 
energy can ignored. 

The graph the function F(y) given Fig. is, course, half the line 
profile for the photon energy which equal unity, i.e. fora 0.420; 
the relative width half the maximum height for the Hartree case 1.30%. 
The percentage line-width any other sufficiently large photon energy 
may calculated from formula (43). 


(43) 1.30¢(a)% 
The function defined plotted Fig. Also shown are the two 


Le] 2 4 6 8 10 12 4 16 18 20 a 


Fic. The intrinsic half-width the momentum line, plotted against the energy 
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points computed Dzhelepov (1955). Note that because Dzhelepov quotes 
the relative width the line plotted against the electron kinetic energy, w’, 
our results are not directly comparable; one must use 


(44) 


Because the momentum density used Dzhelepov similar the one 
calculated from the Hartree wave functions, not surprising that the 
results are good agreement. 

Before drawing any conclusions from the results shown Fig. 
desirable discuss the validity the approximation that the probability 
the Compton process independent the initial state the electron. The 
effect the atomic binding forces the electron during the scattering has 
been discussed Heisenberg and Bewilogua (1931); numerical results have 
been supplied Nelms (1953). The scattering cross section differs appre- 
ciably from the free electron value only for small values the scattering 
angle Thus our case the electron can considered free when 
the cross section being calculated. The Klein—Nishina formula for the Comp- 
ton cross section usually quoted applicable the case which the 
electron initially rest. Before using this formula must, therefore, carry 
out Lorentz transformation from the laboratory co-ordinate system one 
which the electron initially rest. the new co-ordinate system the 
energy the photon, a*, will differ slightly from and the path the electron 
after collision will inclined small angle, the incident y-direction. 
Since the changes and are small, shall calculate the effects these 
changes independently. the initial velocity the electron units 
one finds 


the first order (45) becomes 


The fractional change the forward scattering cross section for small 
change may estimated with the aid the logarithmic derivative 


expression (3a) 


The modulus less than unity for all positive values Now, 
since the most probable value cos zero, the average value cos 
for value near the peak the distribution should considerably less 
than 0.01. Hence the change cross section induced the change the 
value under the Lorentz transformation negligible. 

estimating the magnitude the departure from collinearity the ejected 
electron and the incident y-ray, convenient introduce the angle 


( 
— 
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scattering the transformed co-ordinate system, 6*. Starting from formulae 
quoted Davisson (Siegbahn 1955), one can show that 


Because the simple relation between energy and momentum for photon, 
the conservation the scalar product the four vectors 
for the incident and scattered photons yields 


The relation between and has already been given (45); similarly 


The expansion cos powers can obtained from equations (25) 
and (26) differentiation; shall need only the first two terms. 


Substitution (51) and (52) (49) leads 


0.4 degree and degree respectively. see that the differential 
cross sections these values energy and differ from the values 
negligible amounts, one has only inspect the curves computed 
values given Nelms (1953). 

The only remaining errors importance the values Fig. arise 
from the approximate nature the Hartree wave functions. From the curves 
Fig. one learns that the width for unscreened hydrogenic wave functions 
about double that for Hartree wave functions. (The greater width the 
hydrogenic case occurs because the stronger resultant binding forces give rise 
larger kinetic energies.) Since the Hartree wave functions are believed 
considerably more accurate than unscreened hydrogenic wave functions, 
would surprising the error the calculated values were greater 
than 20%. Such accuracy certainly adequate for discussion the role 
intrinsic width determining the resolution existing spectrometers. 

The measured relative half-widths are considerably greater than the 
intrinsic half-width all energy regions. Below Mev there rapid deterior- 
ation the resolution Compton spectrometer, attributable multiple 
scattering electrons the radiator, and, some cases, scattering 
the residual gas the instrument; this deterioration masks the increase 
below higher energies the measured resolution roughly constant, 
e.g. the spectrometer Groshev (1956) the measured resolution for 
photon energies between 1.5 and 8.5 Mev was about compared with 
intrinsic width the order 


TRANSMISSION AND RESOLUTION 


this section examine the problem choosing the parameters 
Compton spectrometer such way that its efficiency maximized for 
operation given resolution. order handle the analysis with reasonable 


3 
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facility necessary make sufficiently bold approximations that the 
various functions involved can obtained simple algebraic form. Thus 
the relationship derived among efficiency, resolution, and photon energy, 
though believed accurate enough for design purposes, not intended 

Let the transmission, defined the fraction the photons incident 
upon the radiator which eject electron which passes through the slit 
(see Fig. and through the axial aperture the effect the radial 
width the slit will considered separately. The transmission can 
approximated the product three factors. The first these, the 
differential cross section for the production Compton electron the 
forward direction, has been given formula (3). The second the average 
solid angle, defined point the radiator the width the slit 
2sro, and the height the aperture F2, Let the angle 
between the y-direction and the projection the plane symmetry the 
initial velocity vector electron whose orbit passes through either edge 
the slit let the angle between the y-direction and the projection 
(on plane perpendicular the plane symmetry and containing the 
y-direction) the initial velocity electron whose orbit touches the 
upper lower edge the aperture Then 


and 

(55) 

Also 


The third factor the projected superficial density radiator electrons 
plane perpendicular the y-beam. calculate this density shall assume 
that mean radiator atom has the atomic number and the atomic weight 
14. Let represent the thickness, measured perpendicularly its surface, 
the radiator, The projected electron density then approxi- 
mately sec electrons per the curvature the radiator 
not too great; the angle between the central normal the radiator and 
the y-direction (see relation (12)). Collecting all the factors one finds 


Now the efficiency the spectrometer also depends upon the component, 
the area the radiator projected the y-direction. well-collimated 
beam y-rays sufficient cross-sectional area available, the counting rate 
will proportional the product and quantity which shall 
term the luminosity, the maximum value the parameter corre- 
sponding irradiated point the radiator, then the effective height 
the radiator shall assume that the effective width the radiator 
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defined baffles between the two slits which limit the radial departure 
orbits from the optic circle; let the magnitude the maximum permitted 
departure the value not too large, and the angle not 
too large, good approximation the effective width the radiator 
sin this distance the width the electron beam the direction 


and 


The central problem spectrometer design choose values the five 
width. may seen, the line-width increases when any the five parameters 
increased. Our next task derive simple approximate formula relating 
line-width these parameters. most this section line-width may 
assumed mean instrumental line-width, the intrinsic line-width Section 
being added the end the calculation. Imagine detector height 
and infinitesimal width placed the magnetic field varied 
while the instrument receiving well-collimated beam monoenergetic 
photons, then the total under the observed line equal the product 
the luminosity and the flux density photons the radiator. For con- 
venience shall refer the line observed with infinitesimal exit slit 
the image, reserving the term line for the case exit slit finite dimensions. 
The symbol will used for the total basal width the image, relative 
the magnetic field which causes the electrons maximum momentum 

The sources image broadening fall into two classes: the ordinary aberra- 
tions and those resulting from the greater curvature 
the orbits the lower-energy electrons emitted angle the y-direction. 
shall begin calculating the contribution the former class the basal 
width the image; relative momentum is, course, equivalent relative 
magnetic field strength. The contribution from the ordinary aberrations 


The factor relates relative radial displacement the focal plane relative 
changes momentum (Lee-Whiting 1957). The first term the bracket 
(60) the result the finite width the slit The remaining terms 
correspond higher-order aberrations; formulae for the coefficients are 


listed Section 

The multiplying the term present because, with symmetrically 
situated radiator, the contribution the basal width aberration which 
either 2C, depending whether even odd, where the co- 


— 


— 


a 
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manipulating formulae given Davisson (Siegbahn 1955), one can 
show that for electron emitted small angle the y-direction the 
fractional change momentum given 


recall that the kinetic energy the electron units terms 
the angles and defined earlier this section, the maximum value 
permitted the geometry the system Since the term 
may discard the Note that depends only upon the height the 
aperture and not upon the height the radiator; this result con- 
sequence the first-order orbit theory. 

The geometrical limit described above can exceeded for those 
electrons which undergo appreciable scattering the radiator. This effect 
the principal reason for limiting the thickness the radiator. The distribution 
scattering angles calculated the simpler theories multiple scattering 
proportional the scattering angle and Gaussian function this 
angle (see Bethe and Ashkin Segre (1953)). The mean value the square 
the scattering angle for electron which has passed through mg/cm? 


deriving (62) slowly varying logarithmic factor has been approximated 
mean value appropriate mg/cm? and photon energies the 
range Mev. Since the electrons can assumed produced uni- 
formly throughout the radiator, and since the normal the radiator makes 
angle approximately with the initial electron direction, have 


(63) 


order get effective basal width add (60), shall approximate 
the shape the distribution function for scattering angle triangle with 
the same initial slope, the same area, and the same value the mean the 
square the scattering angle. For such distribution function the square 
length the base, related the mean the square the scattering 


angle 

Using the calculated from (62) and (64), shall assume that the 
expression (61) may approximated The total contribution 


the basal image width resulting from the acceptance electrons ejected 
angle the y-direction thus 


The thickness the radiator also affects the image width through the 
reduction the energy the electron multiple scattering. Using approxima- 
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tions similar those the preceding paragraphs, one can show that the 
fractional reduction momentum the most probable energy loss (Segre 


1953, page 256) given 

Our estimate the importance this term probably too small, because 
the most probable energy loss has been used instead some effective maximum 
value. Since (66) always much smaller than the term (65), the error 


should not serious. 
The discussion the preceding paragraphs enables express the basal 
width the image, function the five parameters which the 


luminosity proportional. Thus 


with 

(68a) 

(68c) 

(68d) 

(68f) 


first step our program maximizing the counting rate for given 
resolution, shall maximize the luminosity (59), holding the basal width 
the image (67) fixed. not difficult show that maximum attained 
for the following values the parameters: 


(69a) 


The corresponding formula for the maximized luminosity, which shall 
call Lo, 


this point let return consideration the best choice the field 


parameter The product varies only from 0.3 0.4 
increased from 0.1 0.2. Figure shows one that and and 
hence and not vary enormously rapidly with The two terms 
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(68e) vary opposite ways; low values their effects tend cancel; 
very high values the resulting variation tends counteract that 
The coefficient does not depend upon The only remaining factor, 
depends upon higher fieid-expansion coefficients well upon and 
not likely the designer’s disposal; the fact that root involved 
decreases the probability any great variation. Our conclusion then that, 
far known factors the luminosity are concerned, there significantly 
best value 

One should bear clearly mind that the conclusion the last paragraph 
depends upon the assumption that fixed and that well-collimated beams 
y-rays are available with cross sections large required. For given 
value true, the width the radiator increases when increased. 
But, the other hand, the solid angle which electrons are accepted decreases 
aberration. These two effects cancel approximately the circumstances 
have assumed. 

for given instrument the cross section the available beam smaller 
than the optimum size for given value then advantageous 
increase the solid angle which electrons are collected decreasing the 
value i.e. increasing One must not, however, make small 
that the second-order aberration caused the acceptance lower-energy 
electrons ejected angle the y-direction greater than the first-order 
aberration due finite slit-width There are, however, other reasons 
more pressing for preventing from becoming too small. 

becomes very small the angle between the central radiator normal and 
the y-direction approaches 90°; fact, when small (see 
formula (12)). Since the optimum value the radiator thickness proportional 
cos the radiator will become very thin, and hence difficult make, 
approaches More careful shaping and positioning the radiator would 
also required for smaller values 

Another point considered the choice the necessity keep 
the background from scattering low possible. One way this 
keep the detector far from the y-beam possible; this accomplished 
allowing close possible. the other hand, one must ensure 
reasonable clearance between the y-beam and the slit reduce scat- 
tering; this necessity argues against too small values 

reasonable compromise the face all the effects discussed, the value 


suggested. The corresponding value 44.8°. For values 
closer —0.20, the optimum radiator thickness begins decrease 
rapidly. 


Our discussion the criteria for fixing the various spectrometer parameters 
based upon the assumption that the beam y-rays perfectly collimated 
and adequate cross section. practice the idealized beam simulated 
inserting system longitudinal between the radiator and broad 
source; the cross section the source should cover the projection the 
radiator. Let suppose that only those y-rays whose directions lie within 
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cone half-angle about the nominal y-direction are permitted fall 
upon the radiator. Then the fraction the photons leaving the source which 
are accepted the spectrometer approximately Thus, ignoring the 
absorption y-rays the source, the current electrons arriving the 
final slit the spectrometer the specific activity the 
source photons per unit area. The divergence the y-beam will also affect 
the resolution. One can approximate this effect adding the term 
the right side equation (67). the design the collimating system 
included that the spectrometer, then one should maximize instead 
and add the list parameters chosen. Because and appear 
the function being maximized and the constraint the same way, the 
best value may obtained from the relation 


(71) 


The ratios the terms (67) may still calculated from the relations (69), 
but the normalization must altered. 

far have ignored the influence the radial aperture Let the 
width this slit Using our knowledge the dispersion 
meter, can compute the relative change momentum, 2e, corresponding 
change 2s’ the value the focal plane. 


, 


should clear that increase the momentum width the exit slit 
beyond the basal width the image (i.e. making does not result any 
increase the counting rate. The best value for given value depends, 
course, upon the distribution intensity the image, i.e. upon the line- 
shape observed with detector infinitesimal width. the distribution 
intensity the image uniform, the best slit-width that corresponding 
the width the line half the peak height and the peak counting 
rate for unit incident flux density equal the luminosity. actual 
spectrometer the distribution the image almost certainly strongly peaked, 
since most the broadening arises from quadratic terms. 

For present purposes not unreasonable assume that the graph the 
intensity distribution the image isosceles triangle, with base length 
Let the half-width the line observed with exit slit equivalent 
fractional change momentum and let the counting rate 
the peak the line for y-beam unit flux density. shall use the symbol 
for the ratio elementary calculation leads 


maximizing for fixed assume that has already been maximized 
with 3>A>7/2. Now, held constant, (72a) defines function 


a | 
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examining dC/dx, that the region 0<x<2/5 the function C(x) increases 
monotonically for all values interest. maximum exists 
for value the range this value say xo, may obtained 
solving the equation 


Recall that and that With these restrictions the solution 
(73), accuracy better than 1%, 


Combining (69a) with the definitions and find that the best value 
for the ratio the widths the slits and 


(75) 


the corresponding value may calculated substituting (74) 
(726). 

example the use the various criteria that have been developed, 
let work out the parameters for spectrometer with —0.14 having 
instrumental half-width 0.2% for perfectly collimated y-beam. From 
equations (16) and (17) one can compute the angles specifying the positions 
the slits: 46.42° and 240.5°. The value for which 
zero —0.505. Let assume that other field expansion coefficients are 
the disposal the designer, i.e. that The optimum value the 
ratio 0.431, and the corresponding value the parameter 
The half-widths the first and second slits may found from 

The optimum value for depends upon which turn depends upon 
the unknown parameters a3. Making use (24) one finds 
shall assume that the field variation essentially quadratic the region 
the optic circle, i.e. that small, and put equal unity; since the 
cube root that enters into the calculation, probable that large 
error has been made here. The resulting value 0.069. 

The best values and depend upon the photon energy which the 
optimization carried out; let assume photon energy about 
1.5 assumed for the material the radiator, thickness 2.5 
obtained; since radiator thin has been used Compton spectro- 
meter (Dzhelepov and Shestopalova 1956), our calculated value not 
impossibly small. 

The values the effective radiator area and the transmission for the 
instrumental half-width 0.2% are: 


The instrumental half-width depends upon photon energy, because the 
dependence and upon the electron energy. curve (a) Fig. 
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Fic. variation with photon energy the resolution and efficiency Compton 
spectrometer optimized Curve (a) shows the relative half-width 
the momentum line, curve the sum the and intrinsic half-widths. (c) 
shown the efficiency the number electrons passing through the exit slit per incident 
photon the peak the line. 


the variation the instrumental half-width with displayed. Also shown 
the sum the intrinsic and instrumental half-widths; this sum upper 
bound the combined effect the two types broadening. Note that the 
range Mev photon energy the effective resolution falls from 0.35% 
0.2%; the intrinsic width makes appreciable contribution. The sharp 
rise curve (a) decreased below chiefly caused the increasing 
amount change direction for electrons multiply scattered the radiator; 
the slow rise (a) for larger values originate the increasing relative 
departure from the maximum momentum for electrons ejected fixed angle 
the y-direction. 

The quantity usually quoted measure the efficiency, Compton 
spectrometer the number electrons counted the peak the line per 
incident photon. one assumes that the electron detector 100% efficient, 
spectrometer parameters are maintained their optimum values for 
will vary with because the variation the forward Compton cross section, 
and because the dependence line-shape upon energy; curve (c) Fig. 
shows the dependence upon photon energy. 

The curves Fig. will somewhat modified when the divergence the 
incident y-beam taken into account. estimate the optimum value 
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the half-angle the cone divergence can made from equation (71). For 
our example 1.2°, and the fraction y-rays leaving the source which get 
through the collimating slits about The curve half-width against 
photon energy shown Fig. should moved slightly upward, diver- 
gence large assumed. The neglect this effect least partly 
compensated the fact that the intrinsic and instrumental half-widths have 
been combined additively get the result Fig. 10. 

The normal the radiator the central point is, the example have 
been discussing, inclined angle 44.8° the y-beam. very great 
curvature the radiator needed; the normal the inner end the radiator 
more inclined the y-beam than the central normal about 13°. The 
exact value the effective width which the radiator presents the y-beam 
agrees within with the value this width calculated means the 
approximation, sin used deriving the expression for the luminosity. 

One the approximations made Section was apply the axial focusing 
condition only points lying the central generator the radiator. Electrons 
originating other points will arrive after having been turned through 
angle slightly different from Electrons emitted from such points with 
zero axial component momentum will not have this error 
equivalent effective increase the height the aperture Now 
the radiator lies within the sector defined the lines +3°. The corres- 
ponding increase effective slit height sin 0.0012. Since the 
increase less than the height the aperture should not lead 
any serious additional line broadening. 

The second-order axial aberrations also contribute increased effective 
aperture height Substituting numerical values the formulae (23), 
one finds 


(76) =... —2.98 HT+.... 


This effect equivalent 50% increase the effective aperture height, 
but only for one quarter the electrons emitted from the extreme ends the 
upper and lower edges the radiator the maximum permitted angle the 
plane symmetry. 

Another approximation Section was that the radiator surface could 
section cylinder, instead piece warped surface. Because this 
approximation, orbits originating above below the plane symmetry, will 
not pass through Another aberration occurs for those electrons 
which are ejected small angle the plane symmetry; for the moment 
shall ignore the effect the change energy. The relevant terms 
calculated the methods used Lee-Whiting (1957), are: 
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Inserting the numerical values for the special case that have been discussing, 
and the change only less than 10% the width 
the slit The effect this aberration cause small decrease the 
probability that electrons emitted from the extreme upper and lower edges 
the radiator will pass through the slit 

There yet another aberration affecting the focus electron leaving 
the radiator with its momentum decreased fraction below the nominal 
value will have additional aberration given 


For the special case under consideration the aberration The quantity 
radial shift about 14% the width the slit Note that this effect 
applies almost entirely electrons which have undergone the maximum 
amount multiple scattering, and may therefore beneficial. 

deriving the equation for the shape the radiator the first-order orbit 
equation was used. The effect the second-order correction term, (10c), 
the width the slit would probably desirable include 
(10c) the f(@) equation (9a). 

Albeit that the design that have been considering intended for high 
resolution, where the axial focusing necessary, comparison with the 
instrument Groshev al. (1956) not without interest. These authors 
achieved peak counting rate electrons per incident photon 
resolution 2.7% for photon 1.33 Mev; the area the radiator was 
cm? and the radius the spectrometer was in. Scaling from the 
numerical example that has been discussed and using in. find 
cm? and electrons per photon 2.7% resolution. 
Thus the estimated efficiency instrument the design suggested this 
paper about three times the measured efficiency the spectrometer 
Groshev under comparable conditions; large part the improvement 
probably caused the symmetric position the radiator the optic 
circle (compare Fig. and Fig. 2). 


REALIZATION THE MAGNETIC FIELD 


The value for the proposed Compton spectrometer, —0.14, lies between 
the values for the 180° and the 8-spectrometers, and respectively. 
Since fields for both the latter spectrometers have been obtained without much 
difficulty, there apparent reason why the field for the Compton spectro- 
meter should not realized. Indeed, the problem should even easier than 
that the spectrometer, since any value fairly wide range 
satisfactory. 

The method used Lee-Whiting and Taylor (1956) for calculating the 
parameters system pairs air-cored coils producing desired field 
variation should also applicable this case. The use air-cored coils has 
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the advantage that calculations high accuracy are possible, and the dis- 
advantage that the field expansion coefficients sometimes 
desirably large for higher values the case Compton spectrometer 
which used the neighborhood reactor, there the further 
difficulty that the earth’s magnetic field, distorted nearby iron, would 
difficult cancel. 

designing the pole pieces iron spectrometer one can usually assume 
that the permeability the iron sufficiently large that the magnetostatic 
potential constant the pole faces. the exciting coil and the outer edge 
magnet are sufficiently remote from the region interest, the problem 
reduces one finding the shape two symmetrically placed, equipotential 
surfaces giving desired field distribution the plane symmetry. the 
remainder this section discuss simple method finding these equi- 
potential surfaces. 

Consider finite magnetic dipole lying the axis symmetry the 
spectrometer, with one pole strength distance above the plane 
symmetry, the other strength distance below. The potential 


(79) 
The magnetic field the plane symmetry given 


(80) 


The quantities which one wishes fix spectrometer design are the expansion 
coefficients about the point ro, defined (6); the first few these 
coefficients are: 


(82) 


and only are have given values, sufficient use only two 
dipoles; let stand for the ratio have three variables, and 
choose such way satisfy two equations, (8la) and (816). Let 
hold fixed and seek solutions for and Define 


(83) Y = (3+4a2)B,? —2a28,']. 


i 
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After eliminating one gets quadratic equation 83. Using the fact that 
satisfies the same equation, one finds that 


The corresponding value 


Solutions the problem exist for all values and for which the right 
side (84) positive. good choice for one which does not give large 
values higher derivatives. One can guard against this danger calculating 
the value and plotting the resulting field the region interest. 

obtain the equipotential surface return (79) and calculate some 
point through which desired have the surface pass. Other points may 
found the numerical solution algebraic equation; the process 
simple and rapid, though somewhat tedious. necessary check that the 
free poles the axis not lie between the equipotential surfaces chosen. 

good fit the field form the example that discussed Section 
—0.14 and —0.505—is obtained with: 


By = 0.7 Bo = 2.448 Q2/Q, = — 18.84 


The corresponding value 0.47. 

probably necessary make approximate correction for the fact 
that the permeability the iron not infinite. The magnitude this correc- 
tion will depend upon the dimensions the magnet, well upon the 
magnitude the permeability. Because the uncertainty the calculated 
value probably unwise rely the possibility cancelling out 
the third-order aberration proportional 


SUMMARY 


Section the principle two-directional focusing Compton-electron 
spectrometer developed, and the required form the magnetic field 
derived. equation obtained for the shape the radiator. The intrinsic 
line-width caused the initial motion the atomic electrons, which limits 
the ultimate resolution attainable with Compton spectrometer, calculated 
function photon energy Section The relationship between trans- 
mission and resolution considered Section and attempt made 
calculate the optimum values various spectrometer parameters. Section 
method calculating the shapes iron pole pieces which would produce 
the desired magnetic field described. 
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APPENDIX 
possible set equation for the radiator shape for the case 
point source located large distance, from the center the radiator. 
Then the y-ray incident radiator (see Fig. will make angle with 
the direction the y-ray incident 


and are the co-ordinates The condition for the selection forward- 
ejected electrons which replaces (7) 


(A2) 
leading the equations 


for the co-ordinate point the radiator. Equations (A3) must solved 
numerically, even for the first-order form 

The heights the radiator and the source may increased considerably, 
system slits provided prevent the y-rays from striking the radiator 
too large angle the plane symmetry. 
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THE THERMAL NEUTRON CAPTURE CROSS SECTION AND 
RESONANCE CAPTURE INTEGRAL PROTACTINIUM-233! 


ABSTRACT 


The thermal neutron capture cross section and resonance capture integral 
form the 1.2-minute and 6.7-hour isomeric states have been 
measured using activation methods. The results are follows: 


(1.2 minute) (6.7 hour) 

Thermal neutron capture 

cross section, barns 20+4 19+3 
Resonance integral, including 

1/v part, 0.5 barns 470+90 
Effective cross section 

position, 

0.023), barns 3245 
Total effective cross section 


Irradiations made with different thicknesses cadmium cover show the absence 
major resonance for the formation the 6.7-hour state the region the 
cadmium cutoff. The results are compared with earlier values and some reasons 
for discrepancies have been studied experimentally. 


INTRODUCTION 


neutron irradiation. well known that this conversion has great significance 
the development nuclear power because this means natural thorium 
converted useful reactor fuel, The important nuclear reactions 
which occur during neutron irradiation thorium are given Fig. and the 
relevant nuclear properties the nuclides are given Table Figure shows 
that consumed two processes nuclear reactor. The first, 8-decay, 
leads fissile and governed the half-life for radioactive decay. The 
second process neutron capture which leads and thence 8-decay 
non-fissile isotope uranium. governed the capture cross 
section and the conditions neutron irradiation. The first process 
leads useful product, while the second not only leads non-fissile 
product but also wastes neutrons. knowledge the cross section for 
neutron capture therefore important the reactor designer. 

The significance this cross section was recognized when was shown 
fissile and several measurements have been reported. The first these 
was Katzin and Hagemann (1952). 1946 they compared uranium pro- 
duced from thorium irradiated high and low neutron fluxes and observed 
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(o, = 530b) 


i 


(oe, < 0.7b) 


Fic nuclear reactions the chain nuclear reactor. 


TABLE 


Radiations protactinium and thorium isotopes summarized 
Strominger al. (1958) 


Half-life 


Max. 
energies (Mev) 
and abundances (%) 


energies (Mev) 


and relative 
abundances 


27.0 days 


0.15 (37%) 
0.257 (58%) 
0.568 (5%) 


0.016 
0.028 
0.040 
0.058 
0.075 


0.301 (22) 
0.313 (150) 
0.341 (6) 


1.175 min 
0.58 (1%) 
1.50 (9%) 
2.31 (90%) 
0.043 (2.5) 
0.230 
0.255 (0.10) 
0.770 (0.87) 
0.803 (0.16) 
807 
(0.09) 


6.66 


0.32 (32%) 
1.13 (13%) 


0.043 (20) 
0.099 (14) 
0.125 (15) 
0.153 (7) 
0.225 (13) 
0.293 
0.333 
0.368 (3.5) 
0.566 (7) 
0.732 (9) 
0.803 (0.2) 
0.887 (14) 
0.924 (25) 
1.24 
1.43 (3.5) 
1.68 (2.5) 


24.10 days 


0.100 (35%) 
0.191 (65%) 


752 
232 
23.3 min 
) Pa 234m 
6.66 
0.029 
0.063 
0.091 
0.140 
0.272 
0.376 
0.400 
0.417 
0.476 
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the ratio a-particle emission the fission rate thermal neutron beam. 
The a-particle emission rate reflects the content and the fission rate gives 
the content. They obtained effective neutron capture cross section 
37+14 barns for graphite-moderated reactor. According Halperin 
al. (1955; 1956) and Stoughton and Halperin (1959) the same uranium 
product was mass-analyzed later date and value 55+6 barns was 
obtained. Apparently the integrated neutron flux for this experiment was 
obtained from the reactor operators and the uncertainty this datum not 
included the error given for the effective cross section. 

Halperin a/. (1955; 1956) irradiated thorium for varying lengths time 
the Low Intensity Test Reactor and also unspecified graphite reactor 
and deduced the effective cross section from the 233-234 mass ratio 
the uranium produced. The yield gave the integrated neutron flux and 
the effective capture cross section was taken 8.0 barns. The 
values obtained were 151 and 129 for the LITR and graphite reactor respec- 
tively. The average effective cross section for thermal reactor neutrons was 
given 140+20 barns. 

This completes the information available the time when the present 
investigation was begun. Recently the results another investigation have 
become available. Smith and co-workers the MTR have determined the 
cross section form both isomeric states and their measure- 
ment described U.S. Atomic Energy Commission Report (1955). 
Their method very similar the one used the present investigation. They 
used activation method which samples were irradiated for 
suitable time the reactor and the yields the isomeric states were 
determined counting methods. Irradiations were made with and without 
0.040-in. cadmium covers that epithermal neutron effects could separated 
from thermal effects. The thermal cross sections for formation the 1.2-minute 
and 6.7-hour states were found and 25+4 barns respectively, 
and the corresponding resonance capture integrals were 400 and 270 barns. 

Smith and co-workers employed 0.03 in. diameter alloy 
wires containing 0.52% cobalt measure the cadmium ratio for cobalt the 
irradiation positions they used. The thermal neutron capture cross section 
and resonance capture integral were taken 37.0 barns and 
barns respectively. Recent investigations Johnston (1960) and East- 
wood and Werner (1959) have shown that there large amount self- 
shielding cobalt for neutrons the resonance energy. appears that 
adequate allowance for this effect was not made early measurements the 
cobalt resonance integral and therefore the results were too low. more 
accurate value for dilute cobalt alloys the type used Smith al. has 
been found barns. would appear that the resonance integrals 
obtained Smith should increased the factor 75/41. Self-shielding 
should also considered reactor neutron irradiations made without 
cadmium covers but Smith made allowance for this that their thermal 
neutron values are unaltered. The adjusted values Smith have been 
summarized Table 
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TABLE 


Summary the results Smith (1955), 
corrected for the change the resonance integral cobalt 


Reaction Reaction 
(1.2 min) (6.7 hr) 
Thermal neutron capture cross section, 43+5 
barns 
Resonance capture integral, 
including part, barns 740 470 
EXPERIMENTAL 


pointed out earlier the method used the present investigation was very 
similar that used Smith and his associates. The rates the 
(1.2 minute) and (6.7 hour) reactions were measured 
determining the disintegration-rate ratio samples irradiated 
for measured periods. This method has the advantage that the chemical 
processing the irradiated protactinium need not quantitative because 
using suitable counting techniques the disintegration rates both and 
can determined the same sample. Therefore chemical yields not 
enter the calculation the cross section. order separate the thermal 
neutron and resonance neutron effects some irradiations were made under 
cadmium. Cobalt was used flux monitor. 


All irradiations were made empty lattice position the NRX reactor 
Chalk River (Hurst 1955) using pneumatic move the samples 
into and out the reactor. About 0.01 disintegrations per second) 
the form protactinium nitrate was evaporated the bottom 
small sample capsule which was enclosed iron carrier 
capsule. cobalt wire monitor, 0.005-in. diameter, 0.5-in. long, wrapped 
aluminum foil, was placed the carrier capsule beside the sample capsule. 
Subsidiary experiments showed that the spatial variation flux the irradia- 
tion position was less than 

the 1.2-minute experiments irradiations without cadmium covers 
were made nylon sample capsules because, short irradiation followed 
short decay period, nylon produces less interfering radioactivity than most 
container materials. Nylon capsules could not used for irradiations under 
cadmium owing the high temperature reached these irradiations. Quartz 
capsules were therefore used and minimize breakage they were packed 
the carrier capsule with little powdered graphite. Aluminum sample capsules 
were used for the irradiations. 

The irradiations for the measurement the isomer were about 
two minutes long and the uncertainty the length irradiation had 
the capsule entered the reactor was detected photocell 
which started electric its return journey passed the photocell 
again and stopped the timer. the design the apparatus and the speed 
which the capsule travels was estimated that about one second elapses 
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between the starting the timer and the arrival the capsule the high 
flux position. similar period elapses the return trip. The length the 
irradiation was therefore taken seconds less than the reading the 
timer and uncertainty one second the length irradiation was 
included the estimates the errors. 

The decay time following the irradiation was measured timer which 
started when the capsule passed the photocell the return journey. 

The irradiations for the measurement the 6.7-hour isomer were about 
hours long and timing was not important source error. 


Preparation 

The used for the irradiation was prepared the neutron irradiation 
thorium and was purified ion exchange. Dried powder was 
irradiated for about two days and then after suitable decay period 
was dissolved boiling concentrated After conversion chloride 
was dissolved and applied anion exchange column. The 
thorium was removed washing with and, obtain separation 
from uranium, the protactinium was eluted with redistilled 3.8 After 
evaporation dryness the product was metathesized nitrate form with 
concentrated, redistilled HNO; and stock solution was 
prepared. Suitable amounts the stock solution were taken for irradiation 
and evaporated dryness carefully washed sample capsules. 


Purification 1.2-minute 

The irradiated protactinium was separated from fission products and other 
impurities emitting solvent extraction using 2-thenoyltrifluoro- 
acetone (TTA). The separation was made the sample capsule using 0.5 
0.4 TTA benzene and 0.5 Preliminary experiments showed 
that some and was extracted with the protactinium into the organic 
phase but the experiments reported here this was kept minimum 
adding aluminum and zirconium holdback carriers the The TTA- 
mixture was agitated for seconds using pipette. The mixture was 
then transferred small test-tube and, after the phases had separated, the 
organic layer was transferred preheated standard aluminum counting 
tray. The solution was evaporated dryness and the tray was heated until 
decomposition the organic matter was complete. After cooling the source 
was ready for counting. The time from the end the irradiation the 
beginning counting was minutes. 

When irradiation was made quartz sample capsule specially 
constructed stainless steel mortar and pestle were used crush the capsule 
after irradiation. The extraction with was made the mortar. 
The mixture was then transferred test tube for processing described 
above. 


The protactinium from the irradiations was also purified solvent 
extraction with TTA, but because the the product was longer 
greater care could taken obtain good purification. The protactinium 
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was dissolved from the sample capsule using HNO; containing 
Na, Zr, and holdback carriers. The protactinium was then extracted into 
0.4 TTA benzene mixing the two solutions vigorously with 
pipette for minutes. The phases were allowed separate and the organic 
phase was transferred standard aluminum counting tray where the benzene 
was evaporated. The tray was heated destroy the organic material. After 
cooling was covered with paper tape and counted. 


Measurement 
outline, the method consisted counting the after the had 


decayed, using end-window proportional the efficiency which 
was determined using sources made from solution standardized 

Sources used for the measurement the cross section for formation the 
1.2-minute isomer were covered with paper tape after the 1.2-minute isomer 
had decayed and were allowed decay for least hours before the first 
measurement They were then counted with end-window methane- 
flow proportional The decay was observed for half-lives 
sure that was radiochemically pure and the decay rate the end the 
irradiation was obtained extrapolation the decay curve. allow for 
minor day-to-day variations the efficiency the counter standard was 
counted after the and the count was corrected accordingly. 

The efficiency the end-window counter was determined using sources 
made from solution prepared described earlier section and 
counted methane-flow proportional counter. Redistilled reagents were 
used minimize source thickness. Four aliquots the solution, ranging 
size from wl, were evaporated gold-coated VYNS films with 
total thickness about (Pate and Yaffe 1955). The aliquots were 
made cover about the same area that the source thickness varied, 
although the variation would not factor because dissolved material 
agglomerates evaporation. was found that the counting rate per unit 
volume solution was independent sample size. This not unusual 
observation and indicates that the sources are thin. Merritt (1959), for 
example, found that with the maximum energy which close 
that the main group the amount self-absorption constant 
mean superficial thicknesses below 

obtain the disintegration rate from the counter data, self-absorption 
corrections must considered. This subject has been discussed several 
papers, the most pertinent which Merritt They have established 
relationship between the amount self-absorption and the maximum 
energy for sources about mean superficial density prepared 
evaporation. The correction obtained for the self-absorption the 
radiations about 17% using this relationship. Considerable internal con- 
version follows the decay and the conversion electrons have 
sufficient energy detected with high efficiency the counter. Also, 
they are coincidence with the particles. The presence these conversion 
electrons reduces the self-absorption correction about 4%. The reduction 


a 
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cannot estimated exactly because the decay poorly understood. 
Fortunately, the uncertainty the amount self-absorption not carried 
directly into the cross section result because the disintegration 
rate ratio which enters the calculation. Nevertheless, will indicated 
later section which all the errors are considered, the uncertainties the 
B-counting account for most the error the rate reaction shown 
Tables and 

Absorption the the VYNS film was studied 
adding second layer gold-coated VYNS the back each source. 
change counting rate was observed and was concluded that absorption 
the source mount was negligible. 

Five aliquots the solution were also evaporated standard 
aluminum counting trays. The sources were covered with paper tape just 
were the irradiated protactinium sources, and the counting rates were measured 
with the end-window counter. Because total about mg/cm? absorber 
was present this counting arrangement differences the thicknesses the 
sources used for calibration and the irradiated protactinium sources themselves 
were considered unimportant. 

The efficiency the end-window counter for detecting was thus found 
0.231%. The standard deviation this result was 2%. 


Measurement 1.2-minute 

The content the irradiated protactinium was determined using 
anthracene crystal detector. The bias the amplifier was adjusted that 
pulses from the high energy were accepted and smaller 
pulses from were rejected. The efficiency the detector was determined 
using sources 1.2-minute equilibrium with 24-day prepared 
from solutions standardized 8-counting. 

The detector was anthracene crystal in. diameter, 3/16 in. thick 
coupled DuMont 6292 photomultiplier tube using Silicone oil. The crystal 
and tube were covered with 0.0003 in. thick aluminum foil. order 
certain that the source, crystal, and photomultiplier tube were coaxial, 
brass positioning ring was placed over the assembly. 

Standard aluminum trays were used for counting; they are in. 
diameter and have turned-up lip about 3/32 in. high. The lip served 
confine the protactinium circle diameter in. The sources after 
they had been prepared described earlier were turned upside down top 
the detector. The lip determined the source-to-anthracene distance and 
positioning ring placed the source over the center the anthracene. few 
tests showed that the efficiency the detector was independent the location 
the protactinium the tray, probably because the anthracene was larger 
than the area the counting tray. This essential because the protactinium 
could not confined particular area the tray owing time limitations. 

The output the photomultiplier, after passing cathode follower, went 
Franklin Electronic Company DD-2 linear amplifier and pulses above 
certain bias were recorded scaler. The counting intervals were measured 
with stop watch and most experiments 15-second intervals were used 
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followed seconds for recording the counts and resetting the apparatus. 
Considerable care was taken select the bias used discriminate between 
the high energy pulses and the less energetic ones from that 
changes gain the apparatus did not have large effect the efficiency. 
Corrections for minor day-to-day variations efficiency were made 
counting standard source 1.2-minute equilibrium with 24-day 
before and after each sample irradiated protactinium. 

The efficiency the anthracene counter was determined comparison 
with its parent were used. The was separated from natural uranium 
using solvent extraction and ion exchange methods (Dyrsseen 1950) and care 
was taken minimize the solid content the solution. The counter 
detected from the decay both and but the anthracene 
counter recorded only the pulses from decay; the electronic bias was such 
that the low energy pulses from were rejected. Furthermore, less than 
the 1.2-minute isomeric state decays 6.7-hour (Zijp al. 1954). 
Therefore the disintegration rate was taken one-half the corrected 
counting rate observed with the 

Four aliquots solution were counted with the anthracene 
counter. set aliquots was also evaporated gold-coated VYNS films 
and counted methane-flow proportional counter. The aliquots the 
set were different sizes varying more than factor and were 
made cover about the same area film. They were therefore different 
superficial thickness although, pointed out earlier, the actual thickness 
source material would not expected vary factor 20. The disintegra- 
tion rate per unit volume solution was found independent sample 
size and was concluded that the sources were thin. The correction for self- 
absorption given Merritt al. for such sources about 15% conversion 
electrons from are ignored about 11% the rather inadequate informa- 
tion conversion electron abundance used. This uncertainty considered 
the estimation errors (Table VI). 

Loss absorption the VYNS film was investigated the 
manner described for and was found that tripling the thickness the 
film had detectable effect the counting rate and correction was 
therefore made. 

second solution was made and equal-size aliquots were 
taken for and for measurement with the anthracene detector. 
The average value the two efficiency determinations was The 
error shown the deviation from the mean. 

compare this method efficiency determination with that employed 
Smith and his associates, their method was also used. Weighed amounts 
finely ground natural with equilibrium were spread aluminum 
and platinum counting trays with the aid little ether. The superficial 
thickness ranged from mg/cm?. These were measured with the anthracene 
counter and comparison the results will given later section. 


EASTWOOD AND WERNER: PROTACTINIUM-233 759 


Measurement 6.7-hour 

The decay the 6.7-hour isomer yields several y-rays energy 
greater than the highest y-ray energy 476 kev (Table I). was 
therefore possible use Nal crystal spectrometer differentiate 
between and 6.7-hour the irradiated protactinium. Because the 
large amounts less energetic radiation from saturate unshielded 
Nal crystal lead absorber was interposed between the source and the 
crystal. The efficiency the spectrometer for was determined using 
sources prepared from solution 6.7-hour which had been standardized 

Nal detector in. diameter was used conjunction with 
RCA 5819 photomultiplier tube. After amplification the pulses from the 
photomultiplier were displayed 24-channel pulse height analyzer. The 
counts obtained measured counting period were plotted give spectrum 
and the sum the counts the energy range 0.80 1.03 Mev was taken 
the count. correct for small changes gain and efficiency the 
apparatus standard was counted before and after each measurement 
irradiated protactinium. The 6.7-hour isomeric state equilibrium 
with could not used standard because the gamma rays from the 
decay the 1.2-minute isomeric state are the same energy range those 
from the 6.7-hour state and are much more abundant. Therefore source 
was used standard. has 83.9-day half-life and gamma rays 0.885 
and 1.119 Mev. The position the y-rays was used determine the 
energy scale and the counting rate the energy range 0.80 1.25 Mev 
was used monitor the detector efficiency. 

The major y-activity the irradiated protactinium after chemical purifica- 
tion was has prominent y-ray 313 kev energy and the source 
strengths were such that this radiation saturated the detector. Lead 
absorbers, usually g/cm? (5/16 5/8 in. thick) were therefore used 
attenuate the y-rays. The y-ray peak was still prominent feature 
the spectrum but contributed less than the initial count rate 
the 0.83-1.03 Mev region. Allowance for this and the counter background was 
made from the decay curve. 

The efficiency the Nal crystal detector for 6.7-hour was determined 
using sources standardized the method. Solutions 
were separated from parent Th™ anion exchange using method based 
the work Kraus and Moore (1950). The was obtained from natural 
uranium described for the 1.2-minute isomer. known fraction each 
solution was evaporated gold-coated VYNS film and the counting rate 
was measured with methane-flow proportional Care was taken 
minimize the solid content the sources and they were assumed thin. 
From Merritt the self-absorption correction would about the 
presence conversion electrons ignored about the scanty informa- 
tion the decay used. The uncertainty thus introduced into the 
rate reaction was included with the other errors (Table XII). Auxiliary 
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experiments similar those described for and 1.2-minute were 
made and was found that absorption 8-particles the film was negligible. 
The remainder each solution was evaporated dryness standard 
aluminum trays and counted with the Nal detector. The standard was 
counted with each measurement and the decay was observed with 
both the and the Nal counter order correct for the presence few 
per cent contamination. 

The efficiency determination was made triplicate and the average 
efficiency obtained using the g/cm? absorber was found (3.89+0.11) 
The error given the maximum deviation from the mean. 


Flux Measurements 

Cobalt wires, 0.005 in. diameter, 1/2 in. long were used flux monitors for 
these experiments. The yield from the short irradiations used for the 
1.2-minute experiments was measured with Nal crystal counter 
and the yield from the longer irradiations the 6.7-hour experiments 
was measured with high pressure ionization chamber. The thermal neutron 
capture cross section was taken 36.65 barns, the average the 
activation and absorption values given the summary Hughes and 
Schwartz (1958) and the resonance integral was taken barns (Johnston 
al. 1960; Eastwood and Werner 1959). The half-life used for decay 
corrections was 5.28 years. 

After decay period about week the cobalt flux monitors used the 
1.2-minute irradiations were weighed and the content was measured. 
well-type Nal crystal coupled DuMont 6292 photomultiplier tube 
served the detector and conventional cathode follower, linear amplifier, 
and scaler completed the apparatus. The bias the amplifier was adjusted 
that only pulses y-rays energy greater than 570 kev were accepted 
for counting. The efficiency the apparatus for Co® was determined using two 
sources prepared from solution which had been standardized 
B-, y-coincidence counting. 

The cobalt monitors from the 6.7-hour irradiations were also allowed 
decay for least week before their weight and content were deter- 
mined. The content these wires was too high measured with the 
Nal counter and high-pressure ionization chamber was used instead. The 
measurements were made placing the cobalt monitor well which entered 
the chamber and the current produced was compared with that from radium 
standard. The equivalent the radium standard was determined using 
sources standardized and y-coincidence counting techniques. 
Further details this method measurement are given Jervis (1957). 

The self-shielding corrections applied the observed 
reaction rate 0.005 in. diameter cobalt wire have been determined 
Eastwood and Werner (1959). They are 1.06+0.03 and 1.64+0.05 for thermal 
and epithermal neutrons respectively. These corrections have been applied 
the observed rates and give the corrected rates which appear Tables 
and VIII. 
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RESULTS 


The thermal neutron capture cross section and the resonance integral were 
caiculated from experimental observations using the convention proposed 
Westcott (1958) and Westcott (1958). The application this convention 
activation cross section methods has been discussed detail elsewhere 
(Eastwood al. 1958; Roy and Roy 1959). The symbols used these references 
have been retained the present work. The rate the 
reaction is: 


(dN 


where the ratio disintegration rates and 
the end the protactinium irradiation (destruction during the 
irradiation negligible these experiments), and are the decay 
constants and respectively, and the length irradiation. The 
effective cross section, defined the equation: 


(1) 


(2) 
and related the thermal (2200 m/sec) cross section the equation: 


The value for infinitely dilute cobalt was taken 1.77 (Eastwood and 
Werner 1959). cadmium ratio experiments the factor related 
the cadmium ratio, follows: 


(4) rJ/(1 /To) = 
The coefficient depends the thickness cadmium used and the appro- 


priate values were taken from Westcott (1958). 
The resonance integral excluding the part given as: 


and obtain the total resonance integral which usually quoted, the 

these experiments assumed that the cross section follows the 
law the thermal neutron region and therefore Furthermore, 
assumed that the epithermal neutron flux per unit energy interval propor- 
tional the irradiation position used. 


Cross Section for the Production 1.2-minute 

Four groups measurements were made and the reactor was steady 
operation during each group. There may have been minor changes operating 
conditions between the groups. Each group consisted one irradiation with 
cadmium cover and one without cadmium. The results these experiments are 
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given Tables III The estimated errors are given Table and 
summary the important conclusions appears Table XIII. Comments 
certain aspects the results follow. 


TABLE III 
Flux measurements with cobalt for the (1.2 min) reaction 


Uncovered Cadmium 


reaction rate reaction rate for cobalt 
TABLE 


Rates (1.2 min) reaction 


(1.2 min) 
disint. rate disint. rate Rate 
Uncovered 

7.833 4.47, 0.6862 

Covered 

24.3; 5.18, 9.24, 
46.1, 10.6; 10.0; 

TABLE 
Cross sections for the (1.2 min) reaction 
for Cadmium for for for 

30.3 2.68 29.5 18.8 491 499 
34.2 2.66 28.1 21.1 527 536 
3.55 20.0 22.1 392 402 
28.3 2.60 29.3 448 455 
Averages 31.0 26.7 19.8 464 473 


The decay irradiated protactinium from uncovered irradiations was 
resolved graphically into three components with half-lives the ranges 
minutes, minutes, and minutes. The first period was 
taken and the counting rate the end the irradiation was 
obtained extrapolation this component. accounted for 70-80% the 
activity the end the irradiation. The other periods could not assigned 
definitely particular nuclides. 
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The counting rate the end the irradiation was also obtained 
fitting three components, one which had the accepted half-life (1.175 
minutes), the decay data using the method least squares. Datatron 
digital computer was used for this calculation. The results obtained agreed 
with those obtained graphically within 10% and there was systematic 
variation the two methods. The average the results obtained the two 
methods was used. 

The decay the protactinium from irradiations under cadmium showed 
only two components, one with half-life the range 1.08 1.25 minutes 
and the other longer than 150 minutes. Because the decay data were simpler 
only the graphical method analysis was used. The short-lived component 
which was taken accounted for 22-88% the total activity. 

summary the estimated errors the measurements given Table VI. 
Errors are considered two types: first, random non-systematic 
errors which are reflected the deviation the individual measurements 


TABLE 
Errors the cross sections the (1.2 min) reaction 


(a) Errors 


Non-systematic errors dev. average 11% 
Systematic errors rate reaction 13% 
Systematic errors nvo 


Quadrature sum 18% 


(6) Errors 


Non-systematic errors dev. average 2’) 
Systematic errors 14% 


Quadrature sum 20% 


from the mean and second, systematic errors which may affect all measure- 
ments the same direction. have taken the non-systematic error the 
average value and twice the standard deviation calculated from 
the results the individual measurements given Tables and 

The most important systematic error the determination the efficiency 
the counters for and These determinations rest upon the 
mentioned earlier, contributions from conversion electrons make the self- 
absorption corrections uncertain but they are about the same for and 
They have therefore been taken equal and, because the uncertainty 
thus introduced, error 10% has been assigned the disintegration rate 
ratio. This, together with systematic errors and the reproduci- 
bility the efficiencies the and counters and the saturation 
factor, when added quadrature gives total systematic 
error 13% the rate reaction. 

The systematic error includes 1.5-barn uncertainty the 
cross section, uncertainty the self-shielding corrections for cobalt, and 
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error measuring the duration the irradiation. The largest part the 
systematic error and listed Table connected with the 
uncertainty about the direction incidence neutrons the cadmium- 
shielded samples and the proper value Westcott’s coefficient use for 
the calculations. 

will observed that the capture cross section and resonance integral for 
the formation the 1.2-minute state are about one half the values 
reported Smith summarized Table Two possible reasons for this 
discrepancy have been studied. The first the loss recoil into, 
“hot reaction with, the nylon walls the sample capsule. order 
determine the magnitude this effect four additional measurements were made 
with the protactinium different environment. the first these 500 
was dissolved the solution before irradiation. The solution 
was evaporated dryness nylon sample capsule and then irradiated. The 
second and third irradiations were the same except that these the 
solution was evaporated and then irradiated 
nylon capsule. the fourth irradiation solution, with 
was evaporated quartz capsule and sealed. After the irradiation the 
capsule was crushed the stainless steel mortar and pestle and the fragments 
were treated with the usual mixture. these experiments recoil 
protactinium atoms would expected dissolve with the bulk the pro- 
tactinium and calcium salts. Isotope fractionation should minimized 
although the possibility that the chemical form which 
extracted TTA under the conditions our experiments not elimi- 
nated. was observed that protactinium from these irradiations was less 
pure than protactinium from the irradiations nylon mentioned earlier. 
Nevertheless, component with half-life minutes was resolved each 
experiment and accounted for 17-60% the total activity. The results 
obtained are shown Table VII. Since the effective cross sections agree 


TABLE VII 


Cross section for the (1.2 min) reaction various substrates 


nvo from 
cobalt monitor (1.2 min) Rate 
10-8 disint. rate disint. rate Pa reaction for Pa 
Substrate (cm~? sec™!) (sec) 10-3 (sect) (sec™!) 10% (cm?) 

Ca(NOs)2 7.78 3.57 1.37 0.686 16.6 21.4 
Ca(NOs)2+ CaCOs 7.62 46.9 26.8 0.686 24.8 32.5 
Ca(NOs)2+CaCOs 7.67 41.6 21.1 0.686 21.9 28.6 
Ca(NOs): in quartz 7.40 86.9 18.7 0.686 9.29 12.6 


reasonably well with the values Table and are any rate larger was 
concluded that preferential loss during irradiation was not significant 
source error. The data from these auxiliary experiments using different 
substrates were not included the calculation the final results. 

The other reason for the discrepancy may the measurement the 
efficiency the detector. Smith a/. used the 1.2-minute 


& 
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equilibrium with weighed amount natural uranium standard whereas 
equilibrium mixture was used the present work. The 
sources used Smith a/. were thick compared the irradiated protactinium 
sources. However, they argued that scattering and absorption effects should 
small for the energetic electrons emitted 1.2-minute Platinum 
mounts were used Smith al. and differences the number back- 
scattered electrons detected from sources different thicknesses might 
expected more important when platinum the backing than when 
aluminum used, although here again the difference might small owing 
the high average energy the Smith allow for 10% 
uncertainty the determination the disintegration rate but 
experimental foundation for this allowance seemed desirable 
intercomparison the methods was made. 

Measurements the efficiency our anthracene counter were made using 
and equilibrium sources prepared and standardized 
described early section. was found that for sources mounted 
aluminum the detection efficiency obtained using sources was 
about higher than that obtained using equilibrium sources 
the type used Smith a/. was also found that the increase the contribu- 
tion backscattered electrons the counting rate going from aluminum- 
platinum-backed sources was greater for the thinner sources 
about 2%. These observations confirm the estimates Smith and show 
that these factors are not responsible for the discrepancy. 


Cross Section for the Production 

Nine experiments were made, five without cadmium and four with 0.025-in. 
cadmium covers. The results obtained from individual experiments are given 
Tables VIII The method obtaining and shown together with 
the results Table XI. summary the estimated errors given Table 
XII and, finally, the important conclusions are presented Table XIII. 
Comments about certain aspects the measurements are given the 
paragraphs which follow. 


TABLE VIII 


Flux measurements with cobalt for (6.7 hr) reaction 


Uncovered Cadmium- 
Cadmium covered Cadmium for 
Expt. rate thickness rate ratiofor ry(T/To) cobalt nvo 


X 101° (sec!) cobalt 102 X 10% (cm?) 10713 (cm~? sec™!) 


number X10! (sec™!) (in.) 


1 26.92 0.025 1.450 18.57 2.476 38.26 7.036 
2 28.84 0.025 1.42; 20.2: 2.204 38.0s 7.574 
3 25.09 0.025 1.352 18.56 2.475 38.26 6.55s 
4 25.97 0.025 1.247 20.83 2.138 38.04 6.827 
5 24.69 0.025 1.415 17.41 2.563 38.30 6.446 
6 27.05 0.025 1.397 19.36 2.370 38.19 7.083 
7 27.98 0.025 1.490 18.76 2.376 38.19 7.320 
8 25.54 0.025 1.407 18.15 2.536 38.30 6.668 
9 26.45 0.030 1.331 19.88 2.250 38.12 6.940 
10 24.9% 0.025 1.338 18.64 2.392 38.19 6.532 
11 23.90 0.040 1.203 19.87 2.340 38.17 6.261 


& 
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TABLE 
Effective cross section for (6.7 hr) reaction 
Uncovered 
Expt. rate disint. rate rate 
5.11, 0.898, 31.4 
7.285 3.802 17.23 26.7 
Average 30.4 
TABLE 
Cadmium ratios and values for the (6.7 hr) reaction 
Cadmium- Calculated 
Paws (6.7 hr) covered uncovered 
Expt. disint. rate disint. rate Pa rate rate Cadmium 1+r 


number X10-* (sec!) 1072 1010 (sec!) X10!9 _ratio So v(T/To)se 


2 1.105 3.092 0.8987 9.207 23.0 2.50 29.4 1.65 
4 86.30 205.8 0.8987 7.879 20.8 2.64 29.8 1.64 
6 2.48: 4.568 0.7573 7.219 21.5 2.98 20.5 1.49 
7 20.64 44.60 0.7064 9.08: 22.3 2.45 28.2 1.67 

Averages 27.0 1.61 
9 10.92 21.85 0.8997 6.601 21.1 3.20 19.5 1.44 
11 9.491 19.63 0.8868 6.926 19.0 2.75 23.7 1.55 


was found that the radiochemical purity the was excellent. all 
but one irradiation without cadmium, was the only short-lived nuclide 
present and accounted for more than 99% the counting rate observed 
the end the irradiation. the exception, was present the extent 
60%. irradiations under cadmium 90-99% the activity observed 
the end the irradiation was three the four experiments. 
unidentified short-lived activity was present the extent 45% one 
experiment. graphical analysis the decay data was therefore entirely 
satisfactory. 

The irradiations were usually hours long and required careful planning 
fit them into the irradiation program the Branch. Furthermore observa- 
tions the decay occupied the counting apparatus for several days 
after each irradiation. For these reasons the irradiations were made intervals 
week more and the reactor operating conditions were usually different 
for each experiment. This necessitated somewhat different approach the 
analysis the results than that used for the 1.2-minute measurements. 
order characterize the neutron spectrum for each experiment auxiliary 
cobalt monitor was irradiated after each protactinium irradiation. the 
protactinium sample and the monitor which accompanied were cadmium- 
covered the auxiliary monitor was not covered during its irradiation and 
vice versa. The rates the reaction uncovered irradiations were calculated 
using equation (1) and from this the effective cross section for each individual 
experiment and the average for all five was obtained (Tables and 
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order arrive the thermal neutron cross section and the resonance 
integral knowledge the cadmium ratio for protactinium required. 
Normally this would obtained from irradiations made under similar reactor 
conditions but since this was not possible the observed rate the 
(6.7-hour) reaction under cadmium was compared with the rate 
without cadmium calculated using the observed and average value the 
effective cross section. The factors and for each covered irradiation 
were then obtained (Table X). average value was then deduced 
dividing the average the average and similarly average was 
calculated from the average and average so. These manipulations are shown 
Table XI. 


TABLE 


Calculation the average and values for the 
(6.7 hr) reaction 


Average value (Table IX) and std. dev. average 30.4+1.0 barns 
Average value (Table and std. dev. average 27.0+2.2 
Average value and std. dev. average 1.61+0.04 
integral, excess over 1/v part, from 0.5 

(average So) (average oo) 452+41 barns 
(resonance integral including part, from 0.5 

0.44 460+41 barns 


second approximation the rate the uncovered reaction the 
individual experiments was calculated from the measured and 
values and the values obtained from the first approximation. The average 
and values deduced from the second approximation agreed with those 
from the first. 

summary the estimated errors given Table XII. Here again 
seen that major source systematic error the determination the rate 


TABLE 


Errors the cross section the (6.7 hr) reaction 


Non-systematic errors (2Xstd. dev. average 
Non-systematic errors 1+rs (2Xstd. dev. average 1+rs) 
Systematic errors rate reaction 


Quadrature sum 16% 


Non-systematic errors dev. average 16% 


Quadrature sum 23% 


reaction. pointed out earlier section large fraction this error 
associated with the determination the disintegration rate and 
The systematic error includes the uncertainty the capture cross 


— = 
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section and the self-shielding corrections for cobalt. The error 
connected with the uncertainty the choice just was the case with 
the 1.2-minute measurements. 

order investigate the position the resonance resonances responsible 
for the large resonance integral two additional experiments were made with 
0.03 and 0.04 in. thick cadmium covers, the effective cutoff energies which 
are about 0.5 and 0.55 respectively. The results are shown for experiments 
and Tables VIII and The values for and determined from these 
experiments agree, within experimental error, with the values obtained using 
0.025-in. cadmium. concluded therefore that there large resonance 
near 0.5 ev. 

SUMMARY AND DISCUSSION 


The results obtained for the neutron capture cross section and resonance 
capture integral the reactions (1.2 minute) and 
(6.7 hour) are summarized Table XIII. The convention 


TABLE 


Summary results 


Reaction Reaction 
(1.2 min) (6.7 hr) 
Thermal neutron capture 20+4 19+3 
cross section, oo, barns 
Resonance capture integral, excess 460+90 450+ 100 


barns 

Resonance capture integral, 470+90 460+ 100 
barns 

Effective cross section NRX 
lattice position 
0.023), barns 

Effective cross section 45+9 
0.05), barns 


proposed Westcott has been used analyze the experimental data. 
this treatment assumed that 1/v absorber the thermal energy 
region and that the epithermal neutron flux proportional 
1/E per unit energy interval. 

Experiments using different thicknesses cadmium cover indicate that 
there large resonance for the formation the 6.7-hour isomeric state 
near the cadmium cutoff (0.5 ev). 

The results for the 6.7-hour agree with those Smith but differ 
from theirs factor for the 1.2-minute isomer. has been shown experi- 
mentally that the discrepancy between our results and those Smith 
cannot ascribed differences the method measuring the produced 
nor the loss recoil ‘‘hot reactions although the possibility 
that the chemical form which not extracted TTA has not 
been eliminated. The reason for the discrepancy remains unresolved. 


=. 
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comparison our results with the earlier results Halperin 
Katzin and Hagemann difficult because the neutron spectra which they 
made their irradiations are not known. assume value 0.05 for 
for their experiments, which not unreasonable for graphite- 
moderated reactor, the effective cross section for formation both states 
calculated from the results Table XIII 92+13 barns. This result 
lower than the value obtained Halperin al. (140+20 barns) but higher 
than the value Katzin and Hagemann (55+6 barns). The presence such 
large resonance means that the effective cross section very dependent 
the neutron spectrum and therefore the lack agreement between the various 
effective cross sections not surprising. 
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THE DISTRIBUTION RADIO-AURORA CENTRAL CANADA! 


ABSTRACT 


Five bistatic v.h.f. radio systems were operated central Canada during 
the I.G.Y. for the purpose detecting auroral ionization. Consistent records 
were obtained for period months and these records have now been analyzed. 
Two types events were detected. The nighttime occur most fre- 
quently the auroral zone and characteristically are observed simultaneously 
two points separated about 300 kilometers. The daytime (S) events occur 
simultaneously over much larger area. The time maximum occurrence 
events becomes later with decreasing latitude whereas the reverse true for 
events. The variation with latitude the occurrence events similar 
that other auroral phenomena. 


INTRODUCTION 


The term has been used indicate the ionospheric ionization, 
associated with auroral disturbances, that gives rise radio reflections the 
very-high and ultra-high frequency ranges. The degree association between 
the radio-aurora and the visible optical aurora still uncertain. Some the 
early work (see, for example, Currie, Forsyth, and Vawter 1953) did suggest 
very close association for certain specific instances but recently Gadsden 
(1959) has concluded, the basis statistical studies, that such associa- 
tion exists. clearly some importance establish the statistical charac- 
teristics the radio-aurora order compare them with those the optical 
aurora. While most observations radio-aurora have been carried out with 
radar (backscatter) systems which the radio transmitter and receiver are 
situated the same station, some studies have been made recently using 
bistatic systems which the transmitter and receiver are separated 
distance about 1000 kilometers (Forsyth and Vogan 1957; Collins and 
Forsyth 1959). These studies were started effort provide simple 
method studying radio-aurora that might free from some the observa- 
tional bias that may operate the case backscatter observations. 
part the Canadian program for the International Geophysical Year the 
existing network bistatic systems was extended include five new circuits 
operating frequency 38.07 Mc/s and five operating frequency 
49.99 Mc/s. The transmitters were located Yellowknife, N.W.T., and the 
receivers were located Baker Lake, Churchill, The Pas, Saskatoon, and 
Sulphur Mountain (near Banff, Alberta). The technical details these 
circuits, together with examples the records and discussion the method 
analysis have been given the earlier paper (Collins and Forsyth 1959), 
but for convenience map central Canada showing the location the 
various stations given Fig. and some details the paths are given 
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BAKER LAKE 


YELLOWKNIFE 


MOUNTAIN 


Fic. Map showing the location the transmission paths central Canada. 


Table Several changes were made the equipment during the some 
because instrumental difficulties and some effort examine more 
carefully particular features the scattered signals. result, the longest 


TABLE 
The radio circuits used for the I.G.Y. study 
Path length, 
Path No. location 
Baker Lake 96.2° 900 
The Pas 53.8° N., 101.3° 1180 
Saskatoon 52.2° N., 106.6° 1200 
Sulphur Mountain §1.2° N., 115.6° 1200 


Note: Common transmitters located at Yellowknife (62.5° N., 114.5° W.); frequencies, 
38.07 Mc/s and 49.99 Mc/s; radiated power, 700 watts. 


continuous period for which consistent recordings are available one about 
months starting December, 1957. The purpose this paper record 
some the results revealed analysis these recordings. 


ANALYSIS THE RECORDINGS 


The earlier work (Collins and Forsyth 1959) established the occurrence 
two classes events, designated and ‘S’ events, that appeared increase 
frequency occurrence the auroral zone was approached. Each these 
events was indicated the signal record period strong enhancement 
the received signal strength. The events occurred primarily during the 
daylight hours and were not associated with magnetic disturbances. The 
principal reason for suggesting that they might properly included with the 
other forms radio-aurora was the increase occurrence with increasing 
geomagnetic latitude. The events occurred primarily night and were 
associated with magnetic disturbance. Three subclasses events, designated 


SULPHUR 
SASKATOON 
7 
SUFFIELD CWINNIPEG 
4 
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and were distinguishable the records, which events seemed 
most closely associated with visible aurora, events, always short- 
lived, with the earlier stages, and events with the later stages the auroral 
occurrence. 

For the 5-month period covered the present analysis the bistatic systems 
were intentionally operated very low sensitivity order that only the 
major events would recorded. preliminary analysis revealed that for 
this mode operation very few events were detected. was expected 
from the earlier work, only one two, doubtful, events were detected. 
consequence, was decided retain only the and designations the 
analysis although should borne mind that nearly all, not all, the 
events were actually the type. addition, presumably because the 
reduced sensitivity, very few events were detected the higher frequency. 
The detailed analysis was restricted therefore the lower frequency 
(38.07 Mc/s). 

The available records were divided into 15-minute intervals and the presence 
absence event for each path for each interval was noted. This initial 
record scaling was carried out for the most part one the authors who had 
previous experience with records this type. encouraging note that 
working only from the descriptions and illustrations the earlier paper, 
was able recognize the various events and classify them without difficulty. 
The total number recording intervals for which the records could scaled 
given for each path Table together with the percentage occurrence 
and events. may noted from the table that while the first four 
paths each yielded some 11,000 useful recording intervals, the fifth path 
Mountain) yielded only about 5000 intervals. This 
was because local difficulties with equipment and interference the 
receiving station. 


VARIATION WITH GEOMAGNETIC LATITUDE 


The variation occurrence events with latitude indicated Table 
There slight decrease occurrence for path 


TABLE 
The relative occurrence events for each the paths 


occurrence 


Path No. Recording intervals 
11,712 7.8 16.9 
11,808 12.1 12.7 
11,712 9.3 18.7 
11,808 17.9 
5,232 4.0 


the mid-point which lies close the maximum the auroral zone. 
hard decide the basis the present evidence this decrease significant. 
may that the high occurrence both and events for this path has 
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caused some the events masked events but this seems unlikely. 
Alternatively, there may higher incidence daytime absorption for 
path than for the other paths which would reduce the apparent rate 
occurrence events. Until further confirmation this latitudinal distribu- 
tion obtained probably not wise attempt any further interpretation 
the observed variation. 

The events are more likely behave manner that similar other 
auroral phenomena. interesting compare the observed latitudinal 
variation occurrence events with that derived Vestine (1944) from 
optical and magnetic measurements. The relative occurrence auroral 
phenomena found Vestine for longitudes close that Saskatoon 
plotted against geomagnetic latitude the solid line Fig. the same 


“Radio-aurora 
(A Events) 


RELATIVE OCCURRENCE ——> 


GEOMAGNETIC LATITUDE 


Fic. The relative occurrence radio-aurora events) and visible aurora (after 
Vestine) function geomagnetic latitude. 


figure the relative occurrence events for each the four most northerly 
paths plotted (as cross) against the geomagnetic latitude the path mid- 
point. The value for path (Yellowknife-Sulphur Mountain) omitted 
because its doubtful validity. order fill the diagram for lower 
latitudes the relative occurrence for three the paths used the earlier 
investigation, Churchill-Winnipeg (6), Suffield-Winnipeg (7), and Green- 
(8) are included. Because the different sensitivities and the 
different periods observation the absolute values occurrence for these 
paths can not compared directly with those observed the northern paths. 
For this reason scaling factor was introduced adjust the occurrence rates 
paths and such that the value for path was reasonable agreement 
with the values for paths and 

There seems doubt that the general trend the occurrence 
events function geomagnetic latitude similar that deduced 
Vestine. The observations were made over period only months and the 


4 
(2) 
Aurora (Vestine) \ 
43) 
/ \ 
/ 
(6) 
(7) 
(8) 
| 
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values plotted Fig. might changed substantially any long-time 
average. addition, the previous work has shown that the scattered signals 
not originate exclusively from ionization near the path mid-points. Particu- 
larly for east-west paths, the scattering may take place regions located 
two three hundred kilometers the north the path mid-point. This 
argument might used justify the movement the graph Fig. 
the points (1), (7), and (8) much two three degrees toward higher 
latitudes. this were done the agreement with Vestine’s curve would 
much better. any case, these particular measurements could hardly 
used establish, verify, the precise position the auroral zone; however, 
does seem likely that the zone maximum occurrence radio-aurora 
closely associated with that for optical aurora. 


DIURNAL VARIATIONS 


The diurnal variation occurrence for each the five paths shown 
separately for the and events Fig. Also included are the histograms 


OCCURRENCE TIME) 


LOCAL TIME 


Fic. showing the diurnal variation occurrence events and events 
for each seven paths arranged from top bottom order decreasing latitude. 


found earlier Collins and Forsyth for paths and Each histogram 
plotted with respect local (mean) time the path mid-point, but the 
left-hand column the time occurrence magnetic midnight indicated 
arrow. Clearly, for all paths the events tend concentrate near mid- 
night and the events near noon. For each class event there seems 
systematic shift time maximum occurrence with latitude. For the 
events the maximum occurs progressively later times with decreasing 
latitude and for the events, progressively earlier times with decreasing 
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latitude. Unfortunately, for the range latitudes and longitudes involved, 
not clear the variation time maximum occurrence would more 
systematic measured with respect magnetic time rather than solar time. 
the earlier work, Collins and Forsyth were led the symmetrical distribu- 
tion events about noon (particularly for path believe that the 
events were subject strong solar control. Although solar control may 
operate, some further mechanism required explain the high-latitude 
observations. 

The observed variation time maximum occurrence with latitude 
reminiscent that found other workers for magnetic and ionospheric 
phenomena. Fig. the times maximum occurrence for the and 


LATITUDE 


GEOMAGNETIC 


LOCAL TIME 


Fic. showing the time maximum occurrence events and events 
function geomagnetic latitude. Also shown the latitudinal dependence time maximum 
occurrence for each three other investigations (by Cox and Davies, and Meek) 
possibly related events. 


events are plotted against geomagnetic latitude. Instead the time maxi- 
mum occurrence, the time corresponding the mass’’ the 
distribution may used. When this done the latitudinal variation seems 
somewhat more systematic but hardly enough justify this refinement 
presentation. studying the time maximum occurrence ionospheric 
Cox and Davies (1954) found the latitudinal variation indicated 
the diagram Fig. Nikol’ski (1956) suggested latitudinal variation 
time maximum occurrence magnetic disturbance. The approximate 
variation that found for longitude 105° west also indicated Fig. 
Meek (1955) found that when times maximum magnetic disturbance were 
plotted polar diagram geomagnetic time vs. geomagnetic latitude, two 
opposing spirals were produced. Fig. the relevant portions these two 
spirals have been replotted without correction for differences between solar 
and magnetic time. 
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The points representing the and the events Fig. may fairly 
represented the two straight broken lines. drawing the lines, may 
appropriate attach less importance the points representing the data from 
path because this path located more than degrees longitude the 
east the other paths. While there seems little relation between the 
present results and those the other workers mentioned, interesting 
note that the line representing the events and that representing the results 
Cox and Davies, although displaced about five hours, have nearly the 
same slope. 


GEOGRAPHICAL EXTENT THE DISTURBANCES 


any consideration the nature radio-aurora important have 
some indication the geographical extent single occurrence. The data 
from the I.G.Y. network may used estimate the size the area affected 
each class event and S). assuming once again that the path mid- 
points are representative the regions which the signals are scattered, the 
distances between the various scattering regions can calculated. Leaving 
out path (because the very poor statistics the data from this path) the 
remaining four paths provide range separations from 200 730 kilo- 
meters. For each type event the percentage time simultanous occurrence 
the event each pair adjacent paths was determined. Then the per- 
centage time simultaneous occurrence was determined for three adjacent 
paths and finally for all four paths. These figures were taken measure 
the occurrence the event simultaneously over the whole area represented 
the distance between the corresponding path mid-points. order 
calculate the probability, that event would occur simultaneously 
two paths having mid-points kilometers apart, given that the event 
occurs one them, the percentage times simultaneous occurrence were 
divided the mean percentage time occurrence the event for the paths 


involved. 
Fig. the probability, P(0,x) plotted function the dimension 
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Fic. The probability simultaneous occurrence event two points separated 
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for the and events independently. reasonable estimate the size the 
area affected single ‘occurrence each type event may obtained 
from the value where falls 0.5. this basis single event charac- 
teristically affects area having diameter about 300 kilometers. Similarly, 
the events affect much larger area, possibly having diameter greater than 
800 kilometers. There little justification for assuming circular area for 
each event but, least for the events, the apparent consistency the 
plotted Fig. for widely differing directions mid-point separation 
suggests that there great systematic departure from circular symmetry. 
The use the path mid-point separation measure the diameter the 
affected area may somewhat misleading. Since the area sensitivity 
each path quite large the actual disturbance area may somewhat smaller 
than that indicated Fig. This effect should more pronounced for the 
events, but even for these events would largely offset the aspect 
sensitivity the scatterers. Thus, spite the approximations involved, this 
method measurement should give quite good estimate the size the 
affected area. 

Fig. there much more scatter the points representing the events 
than those representing the events. This may due part the very 
large size the area affected event. The events may fact tend 
more widespread higher latitudes and more localized lower latitudes. 
Certainly, Fig. the two points that lie above the broken line were derived 
from the most northerly combinations paths and the two points that lie 
below the line were derived from the most southerly combinations paths. 


DISCUSSION 


The very limited period time covered the present analysis makes 
impossible draw firm conclusions relating the long-time behavior 
radio-aurora. the other hand, large number events were observed and 
the coverage was virtually continuous both space and time for the area and 
period concerned. seems likely, therefore, that the results give accurate 
picture the actual behavior radio-aurora over this limited time span. 

unfortunate that there not available, the present time, optical 
evidence similar nature covering the same period. The characteristics 
the events seem show marked dissimilarities those optical aurora 
generally. Certainly, the diurnal variations and distributions latitude are 
similar. The geographic extent single optical auroral display still the 
subject some disagreement. There some suggestion that optical displays 
tend occur simultaneously over very large areas. The alternative view 
that while the very largest displays occur simultaneously over large areas, 
most displays are quite localized and may even move systematically latitude 
longitude both. The evidence presented Fig. seems support the 
latter view that, the average, the radio-auroral events are quite 
localized. 

the earlier paper (Collins and Forsyth 1959) the events were tentatively 
associated with the recurrent daytime absorption observed high latitudes 
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and which now known “polar-cap absorption”. The present evidence 
concerning the diurnal variation and geographic extent the events 
accord with this tentative association. interesting note that Little 
and Leinbach (1958) found that daytime absorption measurements stations 
separated 800 kilometers produced correlation coefficient 0.57 while 
nighttime absorption measurements produced lower value. 

The latitudinal variation the time maximum occurrence for both 
events and events seems real but not easily understood. 

For the events the operation solar control not completely excluded 
the present measurements. Assuming solar control, there are least two 
possible explanations the high-latitude lag the time maximum occur- 
rence. the morning absorption observed Cox and Davies fairly 
regular occurrence (see Fig. then high latitudes the morning occurrences 
events would masked and the apparent time maximum occurrence 
would displaced the afternoon. Alternatively, the delay maximum 
occurrence after local noon may due increase high latitudes 
what Appleton (1953) has called the the ionosphere. Accord- 
ing this latter suggestion the observed results could come about through 
decrease recombination coefficient with increasing latitude or, more 
probably, slight decrease the height the scattering region with increasing 
latitude. Appleton has shown that the ‘‘relaxation time’’ the ionosphere 
falls off very rapidly with height the lower ionosphere. Any description 
the basic nature the auroral phenomenon must include explanation 
the diurnal and latitudinal variation occurrence well explanation 
the variation the time maximum occurrence with latitude. would 
valuable, therefore, extend the present study another period time and 
perhaps over larger geographic area. 
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THE GROUP VELOCITY PLANE SURFACE 


MUNGALL AND Morris 


ABSTRACT 


The velocity transmission modulated microwave signal over dielectric- 
covered plane conductor has been measured over range dielectric thicknesses 
for which the first order surface wave mode known propagate. exact 
theoretical treatment plane surface wave group velocity for both and 
surface wave modes any order has been developed and calculations carried 
out for the first two orders. The experimentally determined signal velocities are 
good agreement with the group velocities calculated for the first order 
mode. Effects possible surface wave propagation the accuracy microwave 
distance-measuring methods are discussed. 


INTRODUCTION 


During the past few years number papers dealing with the propagation 
plane surface waves over dielectric-covered metal surfaces have been 
written (Attwood 1951; Barlow and Cullen 1953; Brown 1953; Cullen 1954; 
Rich 1955; Mungall and Morris 1959). However, apart from theoretical 
discussion Smorgonskiy (1955), very little attention seems have been 
paid the question the speed which information can transmitted 
means modulated surface waves. Smorgonskiy has shown theoretically 
that the group velocity surface waves, which defines the speed which 
modulated wave travels, can very different from the phase velocity. This 
paper will describe experimental measurements the velocity transmission 
frequency-modulated waves over surface wave structure, and will show 
that this velocity good agreement with the theoretically predicted group 
velocity. 

Practically, the group velocity can importance the field precise 
radio distance measurement. this application, the velocity which signal 
travels between the two measuring stations must known few parts per 
million the distance between them determined similar accuracy. 
The experimental case dealt with this paper, that layer sand covering 
conducting surface, one which might encountered practical field 
conditions. 


THEORETICAL 


previous paper (Mungall and Morris 1959) dealt with the experimental 
measurement surface wave phase velocity over sand-covered conducting 
surface and outlined the theoretical treatment the problem. The relation 
between the group velocity and the phase velocity over such structure will 
now discussed. For the and cases, the phase velocity surface 
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wave propagation over dielectric-covered perfectly conducting metal surface 
was given the following expressions: 


(1) and 
(2) TM Bo tan Bol 


(3) cot Bol 


where all symbols are the same the previous paper; 
angular frequency, 
velocity electromagnetic waves free space, 
dielectric constant the dielectric, assumed loss-less, 


thickness the dielectric, 
phase constant wave propagation the dielectric. 


The group velocity may general (Slater 1950) defined 
(4) 
where the phase constant. The phase velocity given 
(5) 


The group velocity can determined approximately the manner 
followed Smorgonskiy for the first mode. This involves replacing the 
exact expression for the relations between and simpler, approximate 
algebraic expression giving explicitly function and The derivation 
then involves differentiation this expression. possible, however, 
derive exact solution for all modes, both and follows. 

From (4) and (5) follows that 


(6) 


This can expressed dimensionless terms, using the relative surface wave 
phase velocity and the relative dielectric depth wavelengths 
where the free space wavelength. Thus equation (6) becomes 


where the relative surface wave group velocity. Substitution (1) 
equations (2) and (3) yields the following expressions for terms 


2 
Cc 
, 2 
(7) 
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differentiation dL/dV may determined and hence obtained, using 
equation (7). 
For modes: 


(10) 


where 
Ke KoV -1 


and for the modes 


where 


Since these equations express terms rather than necessary 
(10), (9) and (11) for various values 

the previous paper, calculations were carried out for both and 
modes for 3.0, value representative dry sand. These results are 
shown Fig. where both the relative phase and group velocities are plotted 


or 


Group Velocity 
—-—-—Phase Velocity 


c 


RELATIVE SURFACE WAVE VELOCITY, 


SAND DEPTH WAVELENGTHS 


Fic. Relative surface wave phase velocity v/c and group velocity functions 
the sand depth wavelengths. Only the first two and modes are shown. 
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functions the relative dielectric depth, For the sake clarity, only 
the first two and modes are shown. 

Examination this figure shows that the group velocity always less than 
the phase velocity. For example, the case the mode assumes 
value only 0.51 the phase velocity for value 0.59. For all modes 
the group velocity passes through minimum before increasing slightly 
approach asymptotically the value Both group and phase velocities 
tend toward this same value for large values 


EXPERIMENTAL 


The group velocity was measured determining the transit time 
frequency-modulated waves over the surface wave structure. This consisted 
copper gauze sheet about feet long and feet wide covered with 
layer dry sand, the thickness which was varied required. This structure 
was completely surrounded walls microwave absorber feet high 
eliminate reflections from objects the laboratory. each end openings 
the walls were made for the transmitting and receiving antennas. Each con- 
sisted in. diameter parabolic reflector with double dipole feed. 

The transit time was measured pair Tellurometer sets. These instru- 
ments employ both frequency modulation and pulse techniques carrier 
frequencies near 3000 Mc/sec. signal transmitted from one, the master, 
the other, the remote, where received and retransmitted slightly 
different carrier frequency back the master. Directions polarization 
the outgoing and returning microwave beams are 90° each other and 
45° the horizontal. Accuracy measurement approximately sec. 

Successful launching the surface waves using the parabolic reflectors 
required adjustment the antennas relative the surface and shaping the 
surface their vicinity. was not possible launch surface wave whose 
amplitude exceeded that the free space component which always accom- 
panies unless the sand layer was tapered gradually from the full thickness 
several feet from the antenna zero thickness immediately front it. 
was also found, mentioned the previous paper, that the launching 
efficiency (power transmitted surface wave relative total power emitted) 
decreased with decreasing surface wave velocity. 

The surface waves were launched two different methods. The first and 
simplest required only that the reflector tilted direct the beam 
downward onto the surface. angle incidence about 60° the surface 
normal was found satisfactory. This technique was only applicable 
for surface waves whose group velocity did not differ from more than 
about 10%. For slower waves the surface wave launching efficiency deterio- 
rated and transmission between the antennas was due predominantly free 
space propagation. 

second, somewhat more complicated technique 2b) was found 
excite surface waves such way that the free space component was relatively 
small even for very slow waves This involved surrounding the 
antennas completely with microwave absorber both sides, back and top, 
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Microwave 


Transmitting and 
Antennas 


Sand 


Metal Sheet 
(a) (b) 


Fic. Diagram showing positioning antennas relative the surface wave structure 
and placing absorbing panels. 


and allowing radiation leave this enclosure only through narrow horizontal 
rectangular slot aperture above the sand surface several feet front the 
antenna. This aperture was bounded below the sand surface, above the 
lower edge the absorbing sheet, shown Fig. and each side 
other absorbing sheets. The slot width was varied raising lowering sheet 
was found that both the slot width and the distance between antenna and 
slot affected the relative power the surface wave and free space wave 
components. the slot was made sufficiently narrow, less than wavelength, 
the free space component was apparently absorbed, leaving the surface wave 
component effectively undisturbed. This did not occur the absorbing panel 
was placed too close the antenna. 

most the experiments described the present paper the slot technique 
was used with slot width about 7.5 This width appeared 
satisfactory compromise between the ideal situation with the entire absence 
free space component, attainable only for very narrow slot, and the 
presence relatively strong one resulting from unlimited reflector 
aperture. 

The measurement the transit time between the two antennas means 
the Tellurometer requires further explanation. This instrument measures 
effectively the phase shift frequency modulation, super- 
imposed carrier, due the time taken for the signal 
travel from one set, the master, the other, the remote, and back. The exact 
system measurement has been discussed another paper (Wadley 1958) 
and will not described detail here. essentially phase comparison 
system, and result, interference between the wave train interest and 
any other results unwanted phase shift. can shown that this altera- 
tion phase due retlected waves, any other interfering waves whose 
amplitude independent frequency over small frequency range, varies 
cyclically with carrier frequency alternately positive and negative. first 
approximation, the resultant phase measured over several cycles variation 
can averaged value close the true phase. 

the present experiment, interference occurs between the surface waves 
and free space waves, whether direct reflected. This results variation 
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the measured transit time; the variation may small shown the 
solid curve Fig. for only slight interference, quite large the dashed 
curve when there free space signal the same order the surface wave. 


Free space transit time 


TRANSIT TIME 


FREQUENCY~ TELLUROMETER DIAL SETTINGS 


Fic. Transit time function carrier frequency: solid curve, amplitude surface 
wave greatly exceeds that the free space wave; dashed curve, amplitude surface wave 


about equal that free space wave. 


Both curves this figure represent transit times for sand depth about 
11.5 mm. The first was obtained with 7.5-cm slot and the second with 
15-cm slot. The solid curve follows very closely the calculated surface wave 
transit time, which increases slightly with frequency over the range measure- 
ment. The dashed curve very erratic, but tends oscillate near the free 
space value. 

Measurements transit time were carried out over sand depths 
mm, the range over which previous work the phase velocity had shown 
that the first mode was easily measurable. The surface wave signal 
strengths for thin sand layers tended similar those occurring direct 
free space transmission over the bare metal sheet. the layer was increased 
thickness, however, the signals received each antenna decreased sharply 
until thickness they were near noise level. The velocity signal 
propagation was calculated from the transit times measured and the values 
are shown Fig. relation the theoretical curve for the group velocity 
the first order mode. All experimental transit times incorporate approxi- 
mate corrections for the different velocities over the short tapered portion 
the sand. The experimental points tend fall below the curve the upper 
portion and above the lower. These deviations result from the fact 
that the depth plotted simply the arithmetic mean measurements 
taken along line between the two antennas. Since the upper end the 
curve depths greater than the average tend affect the velocity more than 
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0.5 © Experiment 


— Theory 


RELATIVE WAVE GROUP VELOCITY 


16 


SAND DEPTH - mm 


Fic. Comparison between theoretical and experimental variation surface wave group 
velocity with sand depth. 


those less than the average, variation sand depth will tend make the 
experimental points drop below the curve. The reverse holds true for the 
lower part the curve. Despite the irregularity the sand depth, however, the 
experimental points fall quite close the theoretical curve and demonstrate 
that the velocities measured approximate the theoretically predicted group 
velocities. 

Although was difficult ascertain the depth with any degree accuracy, 
several measurements for sand depths the vicinity were 
made. These indicated that surface wave with group velocity about 
propagated between the antennas power level similar that 
found for depth 11.5 mm. Since the power transmitted the first 
mode had decreased nearly noise level for depths greater than mm, the 
relatively large power transmission depth suggests the existence 
the second order mode. Equipment designed provide greater uni- 
formity sand depth now being constructed order make possible 
investigation the higher order and modes. 

The fact that was possible excite surface waves simply directing the 
microwave antennas downward toward the sand-covered metal surface 
indicates that these waves may also excited any microwave methods 
distance measurement the field. such waves are excited, clear that 
the transit times arising from surface wave propagation can very different 
from those due free space propagation, indicated Fig. Large errors 
the measured distance could result. 
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IMPROVED ELECTROSTATIC ELECTRON 


MARMET AND LARKIN KERWIN 


ABSTRACT 


cylindrical electrostatic electron selector described which provides 
ampere beam electrons whose energy may varied from and whose 
energy spread less than 100 mv. Space-charge problems involved the con- 
struction the device were overcome the use grids for focusing electrodes 
with exterior electron traps, and non-reflecting surface for electrons made from 
network tiny tubes. The selector provided with energy analyzer. When 
used determine the appearance potential curve argon, the selector resolved 
the and states separated 0.18 ev, and indicated the formation 


INTRODUCTION 


few years ago, Father Ernest Clarke (1954) reported from this laboratory 
the use electron selector measure ionization probability curves. Some 
details the instrument are shown Fig. 


Fic. Electrostatic electron selector (Clarke 1954). 


Electrons from filament are accelerated about 3-volt energy and 
enter the selector through slit plate There the electrostatic field main- 
tained between (1.5 and (4.5 caused the 3-volt electrons describe 
circular path emerge from the exit slit They were then further 
accelerated any desired energy means voltage applied the ioniza- 
tion chamber Here they were used bombard gas molecules entering the 
chamber from The resulting ions were then analyzed mass spectro- 
meter whose accelerating slit system shown and 
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The 127° selector used will not pass electrons other energies than those 
for which adjusted, good approximation. also produces direction 
focusing (Hughes and McMillan 1929), which necessary for useful beam 
intensity. Fig. shown the energy distribution obtained from this selector. 


TRAP CURRENT amp 


12,0 13.0 14,0 
ACCELERATION POTENTIAL (VOLTS) 


Fic. Energy spread electron beam selector. Outer deflection plate potential: 
1.1 volt (6) 1.2 volt (c) 1.3 volt —O—; 1.4 volt —A—. 


Under typical operating conditions effective energy distribution 
300 half-height was obtained. This stems from the measured distribution 
(Fig. 2c) modified the fact that the low-energy electrons are less effective 
producing ions. The electron beam current was about amp, and could 
used from about 7-ev 60-ev energy. With this instrument measurements 
were made leading the determination the dissociation energy nitrogen 
(9.79 ev), the electron affinity oxygen (1.5 ev), and other data. 

would have been desirable increase the electron beam current and 
reduce its lower energy limit well the energy spread. These are opposed 
desiderata and the above characteristics represent the optimum conditions 
obtainable under the circumstances. The principal block improvement was 
the formation space charge reflection from surfaces and scattering the 
gas atmosphere the mass spectrometer ion source. 


IMPROVEMENTS 


Grid Construction 

Since the space charge results part from the reflection electrons from 
surfaces, was thought prevent this making the surfaces transparent 
and capturing the space-charge-creating electrons charged plates when 
they left the selector region. This was tried large scale selector model, 
where the long paths favored space-charge formation. shown Fig. 
the selector was constructed 90% transparent wolfram mesh grids and 
beyond which were outer plates (P; and one experiment the elec- 
trodes were connected circuit shown with connected and 
The following readings were obtained: 
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SETTING SETTING SETTING 


Total emission, 500 500 500 


Grid voltage, 


Comparing settings and seen that the first case appreciable 
current reached the grids, the electrons which passed through the slit being 
reflected back the slit the space charge created the second case 
the voltage the plates resulted very little field penetration through the 
grids may seen from the small grid current. setting the grid voltage 
attracts electrons, which are captured the plates they pass through the 
grids, preventing space-charge build-up. The current previously 
reflected the slit now continues the plates. 


Fic. Grid construction selector electrodes. 


With the plates grounded, electron beam could passed through the 
selector. With the plates volts, the removal space charge resulted 
current about amp with energy spread about 200 mv. was also 
found that the selector could operated very low electron energies (about 
0.5 ev), which important, since the theoretical spread energy propor- 
tional the energy (AE/E AR/R (Hughes and McMillan 1929)). 


Focusing 

The cylindrical electrostatic lens provides direction focusing the planes 
normal the cylindrical walls. certain focusing the surfaces concentric 
with these walls was achieved adding electrodes the top and bottom 
the selector. These were wolfram mesh and were provided with exterior 
plates for the removal space charge the case the selector grids. The 
use this focusing always increased the beam intensity least factor 
These electrodes and plates may seen the photograph the complete 
selector Fig. 


790 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


The large selector constructed for these experiments would not operate well 
high Hg) pressures because the long electron trajectories. Its 
relatively small height also limited the beam strength, which was found 
experimentally vary inversely the volume the selector. 


Electron Velvet 

When small electron selector was made using the above-mentioned 
improvements, its characteristics were measured means 
(the and beam about amp was obtained. 
However, when the mass spectrometer ionization chamber was located the 
selector exit slit, space-charge limitation was encountered anew region 
(Fig. the electrons scattered from the ionization chamber walls. The 
solution grid walls could not applied here, ions would formed 
regions between the grid walls and plates electrons whose energy was 
undetermined. Another possibility lay the finding surface which would 
absorb electrons rather than reflect them. practice these surfaces would 
not physically ‘‘clean’’, as, the atmosphere the mass spectrometer 
source, surface layers gas molecules and products catalysis the 
would always present. 

Measurements were carried out number surfaces shown Fig. 
Electrons from the selector were directed onto surface whence, reflected, 


Fic. Apparatus for measuring reflection coefficients. 


they proceeded surface then back etc. Currents and 
were measured shown. The experiment was made low pressure 
Hg) avoid the effects gas-scattering. 
the percentage electrons reflected from and that from 
then when the current from the selector: 


whence 


(1) 
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When and are identical, this simplifies 
(2) 


result these measurements, the following results were obtained for 
electrons ev: 


COEFFICIENT REFLECTION 


SURFACE +5% 
Iron, polished 
Platinum 
Copper, etched 
Iron 
Gold 
Cadmium 
Palladium 
Copper 
Iron, etched 
Soot metal 
Finned copper 


The high reflectivities obtained here are not variance with the low values 
obtained other workers (e.g. Fowler and Farnsworth 1958) for similar 
metals when the surface conditions are taken into account. Thus attempt 
was made degas the surfaces for our measurements, and obvious that 
the reflectivities obtained are characteristic the gross texture and surface 
conditions ordinarily met high pressures mass spectrometer source 
rather than the metals used. 

From these results was clear that simple surface lent itself the 
absorption slow electrons. Soot (deposited from flame) was relatively good. 
stainless steel ionization chamber coated with soot was used with the 
selector, and the space-charge problem (Fig. immediately reduced 
practical proportions. However, this surface could not tolerated the 
spectrometer. The low reflectivity surface bristling with many small fins 
suggested that this seriously considered. 

After many attempts, good surface was obtained follows. Fine aluminum 
wire (0.5 mm) cleaned and plated with 0.01 copper. Thousands 
6-in. lengths plated wire are inserted aluminum cylinder the 
same length and about diameter. The whole compressed approxi- 
mately oval sections with pressure about 600 kg/cm?, squeezing the wires 
tightly together. The end the cylinder then worked flat and polished, 
when the fine network copper plating may seen the aluminum surface 
formed the ends the squeezed wires. This surface placed for few 
seconds sodium hydroxide, which dissolves little aluminum and leaves 
the copper plating (tubes) slight relief. The surface then heavily plated 
with copper (0.5 mm). This end the cylinder then sliced off form 
plate about thick which ground and polished both sides 
thickness about mm. then immersed sodium hydroxide until all 
the aluminum has been dissolved, leaving the fine copper tubes, about 
long and 0.5 diameter anchored their plated copper base. The sheet 
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Fic. Electron velvet ionization chamber. 


then plated with gold. The appearance this material (named 
Dr. Clarke) shown Fig. where has been hand-worked into 
the shape ionization chamber. Its coefficient reflection measured 


Fic. selector and analyzer. 


= 
Ce XY 
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Fig. seen photograph complete selector, together with 
analyzer. this picture the ionization chamber has been replaced analog 


Fic. Complete electron selector and analyzer. 


mesh structure. the final model the selector described below the space- 
charge-capturing plates outside the mesh grids were also constructed 
electron velvet. 


CHARACTERISTICS 


All potentials are referred mid-filament potential, and should corrected 
for contact potential indicated measurements made the appearance 
potentials gases and the use negative potentials some electrodes. 
strip filament placed about one millimeter from the selector entrance was used. 
Shields were placed about the filament, the selector, the ion chamber, the 
analyzer, and the beam detector eliminate stray-reflected electrons. 


je 
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Because the difficulty knowing the precise value (in addition 
contact potential there was voltage drop along the filament and penetration 
fields), was calibrated practice. 
The electron energy was set ev, not mask the energy spread 
the selector beam the resolution the analyzer. The current the 


TABLE 

Electron selector data 

Selector Analyzer 

Entrance slit width 0.5mm 0.5mm 

Exit slit width 0.5mm 0.5 

Height slits 4mm 4mm 

Initial electron energy, lev 

Potential ionization chamber Variable 


ionization chamber was about amp, and that emerging from the analyzer 
was about amp. sweep the analyzer voltage produced the curve 


Fig. 


w 


The measured energy distribution seen about 0.05 half-height. 

portion this width contributed the resolving power the analyzer, 
that may infer that the electron beam traversing the ionization chamber 
had energy spread half-height not more than about 0.04 ev. This 


I 
ame 4, i 
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the order the thermal agitation gas molecules ordinary source tempera- 
tures, and probably represents practical lower limit beyond which 
would not very useful proceed. 

Measurements appearance potentials, described the next section, 
confirm the attainment this narrow energy spread. 


APPLICATION 


The electron selector was placed bell jar which was evacuated, out- 
gassed, and filled with argon various pressures. The ions formed this gas 
the energy the electron beam provided the selector was increased 
were collected electrode placed few millimeters from the ionization 
chamber and maintained few volts positive with respect it. The resulting 
are shown Fig. which abstracts the salient features 


many runs. 


Fic. curves for argon. 


relatively low pressure Hg) curves indicated were 
obtained. The break following the appearance potential appeared 
average value for runs +0.01 ev. This agrees well with the 0.178-ev 
difference levels between the and the states argon known from 
spectroscopy (Moore 1949). This break has been mentioned Fox (1956) 
and recently Fineman (1960). The resolution these two levels the 
electron selector considered creditable performance.* 

relatively higher pressures Hg), the appearance potential 
curve appeared ascribe the process appearing 0.8+0.1 
below the threshold for (average runs) the formation Horn- 


*Following visit our laboratory inspect the selector, Dr. Conrod, the 
Department Aeronautics and Astronautics, has just advised that they have successfully 


duplicated it. 
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beck and Molnar (1951) have reported the formation this ion these 
pressures considerable abundance, and report its appearance potential 


being below the threshold for At. 
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ANALYSE BANDE 1-0 SYSTEME PH! 


ABSTRACT 


The (1,0) band the system has been photographed absorp- 
tion under high dispersion. Values for the rotational constants, the spin-coupling 
constant, and the vibrational interval AG} the state have been obtained. 
study the A-doubling shows the existence perturbation the state; 
not possible, however, determine unambiguously the nature the perturbing 
state. 


INTRODUCTION 


depuis longtemps été photographiée émission haute dispersion par 
Pearse (1930), puis par Ishaque Pearse (1936, 1939). Ces derniers ont 
deux bandes leur permis déterminer les constantes des deux états 
transition pour plus bas niveau vibration. Jusqu’a présent possédait 
aucune information relative aux autres niveaux vibration. 

présent travail est consacré bande 1-0 3200 qui 
étre photographiée absorption avec une intensité Les con- 
été mis évidence. 


TECHNIQUE EXPERIMENTALE 

Les bandes 0-0 ont été observées par Dr. Johns 
Dr. Ramsay cours d’expériences photolyse par éclair 

Les conditions expérimentales peuvent décrire briévement comme suit: 

Longueur parcours d’absorption: traversées d’une cuve métre. 

Lampe photolyse: métre long, alimentée par une batterie con- 

Lampe source: lampe capillaire quartz alimentée par condensateur 

Pression phosphure d’hydrogéne: Hg. Gaz renouvelé aprés 
chaque éclair. 

Les bandes ont été photographiées dans troisiéme ordre 
spectrographe réseau concave pieds, sur des plaques Eastman Kodak 
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creuse fer dont les longueurs d’onde ont été prises dans les tables 
(Harrison 1939) rapportées vide des corrections d’Edlén (1953). 
Les positions des raies sur plaque photographique bande 1-0 ont été 
mesurées d’un comparateur équipé d’un viseur photoélectrique 
type décrit par Tomkins Fred (1951). n’a pas été jugé utile mesurer 
nouveau bande 0-0. 


ANALYSE ROTATIONNELLE 


des raies bande 1-0 présente aucune difficulté 
Ishaque Pearse. Les nombres d’ondes numération des raies des dif- 
branches sont portés dans Tableau une reproduction spectre 
est donnée sur Fig. 

Sur les branches prévues par théorie pour une transition 
ont été observées. L’allure générale bande 1-0 est trés semblable 
celle bande 0-0. Notons simplement que les raies branche 
ses satellites sont trés serrées chevauchent fortement, comme 
dans une branche type contraire, les raies des branches 
sont trés espacées. Ces deux effets sont plus marqués ici que dans bande 0-0. 

L’analyse bande 1-0 confirmé entiérement des raies 
bande 0-0 par Ishaque Pearse. 

structure théorique été décrite par ces auteurs qui 
donnent schéma des niveaux des transitions permises avec 
Nous reprenons les mémes notations nombre quantique 
qui sera désigné par conformément actuel. trouvera également 
dans une étude spectre par Dixon (1959) schéma trés clair 
type transition. 


utilisant des relations combinaisons entre les raies bande 1-0, 
est possible déterminer directement les constantes Mais Ishaque 
Pearse ont obtenu les constantes avec vraisemblablement une 
précision supérieure celle que peut obtenir sur car, par rapport 
bande 1-0, bande 0-0 est plus intense, les raies chevauchent moins 
sont pas déplacées par effet perturbation. est donc préférable 
les différences entre les nombres d’ondes des raies bande 1-0 
ceux des raies correspondantes bande 0-0, qui permet 
les différences entre les constantes des deux états vibration: 


effet, d’une formule donnée par Gilbert (1936), différence 
entre les sommes des termes pour des deux bandes 
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premier membre est determiné prenant moyenne entre les niveaux 
sera plus loin que les niveaux cet état sont assez fortement perturbés. 
Les valeurs des constantes ont été alors obtenues par méthode des moindres 
carrés ont été portées dans Tableau III, ainsi que les valeurs 
calculées partir des valeurs données par Ishaque Pearse. 

Afin rendre compte différence comportement des niveaux 
est montrée sur Fig. pour les niveaux courbe été 
calculée avec les valeurs Tableau 


° 
~ 
1833.00 
i 1 H i 
0 10 20 30 40 SO 60 7 80 90 


fonction séparément pour chacune des composantes 


suivante, due Gilbert: 


2 


| 
| 
1834.00 | 
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+40 
Les différences entre les termes rotationnels qui apparaissent premier 
membre déterminent combinaisons différences entre les 


nombres d’ondes des raies. moyenne ces différences pour les niveaux 
été portée dans cette équation, mais les valeurs correspondant aux 
niveaux perturbés ont été remplacées par des valeurs obtenues 
supposant les niveaux séparés des niveaux par méme intervalle 
que dans bande 0-0, qui revient supposer perturbation absente 
dédoublement identique dans les deux bandes. Les valeurs portées 
fonction J(J+1) sont montrées sur Fig. courbe titre indicatif, 


| 
4.50 
en 
4,00 


(1) Sans excluant les niveaux 
(2) Sans avec les niveaux 
(3) Avec excluant les niveaux 


porté également les valeurs déterminées utilisant les niveaux 
réels (courbe 2). remarquera que n’est pas constant, mais varie 
réguliérement avec Dixon, étudiant les bandes appartenant 
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systéme correspondant, montré qu’il était nécessaire tenir compte 
termes supplémentaires prévus par Mulliken (1930): 


convenablement constante obtient une valeur indépendante 
Ceci est atteint pour 


0.065 


Les valeurs sont portées dans Tableau III ainsi que les 
valeurs relatives bande 0-0 qui ont été recalculées 
les mesures Ishaque Pearse. Ces auteurs, n’ayant pas introduit con- 
stante yo, trouvent une valeur peu inférieure pour 


Chaque niveau est dédoublé sous-niveaux dont 


peut déterminer des relations 


des relations analogues qui trouvent aisément schéma 
transition. Hebb (1936) calculé théoriquement cet effet fonction trois 
constantes Co, partir ces résultats, Gilbert ajuster les trois 
constantes maniére rendre compte des résultats expérimentaux pour 
bande 0-0. Par contre, dédoublement dans bande 1-0 s’écarte considé- 
rablement des prévisions théoriques; pour les états dédoublement 
est trés peu différent d’une bande mais est considérable pour 
des trois constantes est impossible. Ceci suggére 
que moins des deux ensembles niveaux est perturbé dans 
les différences entre les raies correspondantes des deux bandes. par 
exemple, 


ainsi que d’autres relations méme type. Les valeurs obtenues sont portées 
dans Tableau II. L’importance dans est frappante. 


'S 
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TABLEAU 


Variation dédoublement entre bande 1-0 
bande 0-0 


0.04* 

0.07* 

0.43 0.05 
0.41 0.03 
0.56 0.04 —0.09 
0.43 0.18 
0.34 —0.10 
0.26 
0.26 
0.37* —0.08 


*Détermination douteuse. 


dans quelle mesure les deux niveaux sont également affectés par celle-ci, 
les différences des termes rotationnels entre les deux bandes ont été 
calculées utilisant les formules (1936) qui (Herzberg 


1950, 235) 


Ces formules sont qu’approchées tiennent pas compte dédouble- 
ment Néanmoins elles ont mérite d’une simplicité relative donnent une 
bonne approximation pour les faibles valeurs calcul été fait 
utilisant les valeurs des constantes rotationnelles données dans 
Tableau III. Plusieurs valeurs ont été essayées représenter 


TABLEAU III 
Constantes 


| 
| 


*D'aprés Ishaque et Pearse (1939). 
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= | 
0,50 
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4 A | 
2 a 7) 4 | 
re) 
4 a 
+ 5 6 7 8 9 10 " 12 13 


Fic. observées moins valeurs calculées différence entre les termes rota- 


composantes 


mieux les résultats expérimentaux. Pour 

(3) —15.257 

obtient résultat représenté sur Fig. les valeurs observées moins 
les valeurs calculées des différences, 


ont été portées fonction pour les niveaux pris séparément. Les 
points groupent assez bien autour d’ordonnée nulle, excepté ceux 
relatifs aux niveaux F3, qui écartent considérablement pour plupart. 
Toute autre valeur donne une dispersion plus importante. valeur (3) 
est différente valeur —15.241 calculée partir des valeurs 
Tableau III. Cependant cette différence reste dans marge d’erreur possible, 
l’on tient compte fait que les coefficients pas été introduits 
dans les formules 

L’existence d’une perturbation des niveaux semble assez 
fermement établie. perturbation des autres niveaux est, sinon nulle, 
moins assez faible est difficile prouver son existence avec certitude. 
effet faible variation dédoublement quand passe bande 0-0 
bande pour les états peut s’expliquer aussi bien par une per- 
turbation que par changement des valeurs des constantes intervenant dans 


LEGAY: SYSTEME 805 


théorie Hebb. précision des mesures n’est pas suffisante pour décider. 
dispersion autour zéro des points autres que ceux apparte- 
nant aux niveaux peut également s’interpréter par les raisons précédentes 
aussi par fait que des formules peut-étre 
pas suffisante. 

est difficile déduire des résultats précédents nature pertur- 
bateur. l’on admet que seule une des composantes des doublets est 
perturbée, parmi tous les états donnés par les deux premiéres configurations 
électroniques (voir Herzberg, 341), seuls sont possibles les états 
les régles sélection des perturbations. Cependant une 
perturbation par base semble peu probable, car nombre quan- 
tique vibration devrait atteindre des valeurs trop élevées. Une perturba- 
tion par reste une hypothése plausible, mais est difficile 
fermement données expérimentales plus 
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FORBIDDEN TRANSITIONS DIATOMIC MOLECULES 


THE ROTATION-VIBRATION SPECTRUM THE 
HYDROGEN-DEUTERIDE (HD) MOLECULE! 


ABSTRACT 


The 1-0, 2-0, and 3-0 bands have been observed 
with PbS infrared grating spectrometer and the 3-0 and bands have been 
photographed with concave grating spectrograph. From these spectra 
precise values the vibrational and rotational constants its electronic 
ground state have been determined. The variation with vis similar 
that recently established for that is, the values for and 
above the values corresponding linear extrapolation the subsequent points. 
This leads uncertainty the and values which much greater than 
that the and values. The (very low) intensity the 
spectrum close agreement with that predicted. The decrease 
intensity going from 1-0 4-0 strikingly slow, far slower than normal 
series fundamental and overtone bands. This also qualitative agreement 
with theoretical expectation for molecule like which does not have dipole 
moment its equilibrium position. 


INTRODUCTION 


The spectrum homonuclear diatomic molecule 
one the most strongly forbidden transitions. forbidden dipole 
radiation account (1) the absence permanent dipole moment the 
molecule, (2) the selection rule that electronic states the same parity cannot 
combine with one another (3) the selection rule that symmetric 
rotational levels not combine with antisymmetric ones. The spectrum can 
occur the free molecule only quadrupole radiation which many orders 
magnitude weaker than dipole radiation. This particularly long 
wavelengths since the intensity inversely proportional rather than 
The intensity the quadrupole spectrum was predicted the 
order the normal dipole intensity infrared spectra hetero- 
nuclear molecules. spite this low intensity, using extremely long 
paths (up atm), has been observed described paper this 
series (Herzberg 1949, 
homonuclear molecule the three selection rules just given are not 
independent. But heteronuclear molecule consisting two isotopes 
the same atomic number, such HD, while the rules (1) and (3) are longer 
applicable, the electronic selection rule (2) still applies and may considered 
the reason for the extremely low intensity the rotation—vibration spectra 
such isotopic molecules. However, rule (2) when not supported 
not absolute rule because the interaction electronic motion and 
received March 14, 1960. 
from the Division Pure Physics, National Research Council, Ottawa, 
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rotation and vibration. account this interaction, dipole transitions 
contradicting rule (2) will occur with small intensity. Alternatively, one may 
ascribe the occurrence these transitions the presence small dipole 
moment which arises when the molecule vibrating (i.e. account inter- 
action vibration and electronic motion) but absent the equilibrium 
position. These transitions belong therefore the group forbidden transi- 
tions which occur because the particular selection rule holds only certain 
approximation [type (1) (Herzberg 1950a)]. They are the same type 
for example transitions which are forbidden the electronic selection 
rules but which can occur very weakly account the interaction rotation 
and electronic motion. 

The simplest spectrum the type here considered that the mole- 
cule. estimate the intensity this spectrum was first made Wick 
(1935) and was first observed this laboratory (Herzberg means 
the second and third overtone the photographic infrared. Since then, 
measurements the fundamental and first and second overtone have been 
made with high resolution infrared spectrometer. the present paper all 
the measurements from 3450 13,680 are collected and evaluated. 


EXPERIMENTAL 


order obtain the photographic infrared spectrum White multiple 
reflection cell similar to, but smaller than, that used earlier papers (Bernstein 
and Herzberg 1948) was used. This cell was made glass tube long and 
wide and was fitted with set mirrors radius curvature. 
With 200 traversals the intensity was still sufficient obtain spectra 
with 21-ft concave grating spectrograph. The absorbing paths were thus 

For the studies the region similar absorption tube 
length and diameter but made brass was used. With photo- 
electric detection and electronic recording the number traversals that can 
used limited the fall the signal-to-noise ratio. The optimum number 
traversals dependent the reflectivity the mirrors and for the tube 
used was found traversals. Pressures atm were used, permitting 
maximum path 360 atm. 

The light sources used were zirconium concentrated arc the photo- 
graphic infrared and tungsten filament lamp the region. 

The photographic infrared spectrum was obtained the first order the 
N.R.C. 21-ft grating spectrograph I-N and I-M Eastman Kodak infrared 
plates. Exposure times ranged from hours. iron are was used the 
source comparison spectrum. 

The spectrum was obtained with the infrared spectrometer 
described Douglas and Sharma (1953), with 7200 lines/inch grating 
first, second, and third orders and with 15,000 lines/inch grating first order 
using suitable filters isolate the orders. Interferometrically determined Hg, 
Ne, Ar, and lines (Burns and Adams 1953; Burns, Adams, and Longwell 
1950; Meggers and Humphreys 1934) were used higher orders standards 
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for the lower-order infrared spectra. Interference fringes produced 
etalon continuous spectrum were used interpolate between 
these standards. The vacuum corrections Edlén (1953) were applied. Calibra- 
tion and spectrum scans were recorded simultaneously (with photo- 
multiplier detector for the former) illuminating different sections the 
entrance slit. Careful corrections were made for the shift produced the 
curvature the entrance slit image the spectrometer the exit slit (see 
Minkowski 1942). 

Measurements different plates the photographic infrared and 
different recordings the lead-sulphide region show agreement within 
given band +0.02 The possibility slight systematic errors remains 
due very slight deviations from the integral relations between different 
orders, the photographic region due slight shift between comparison 
and absorption spectrum account temperature barometric pressure 
changes. Thus there systematic difference —0.04 between the 
photographic and lead-sulphide wave numbers the 3-0 band. The absolute 
accuracy therefore estimated the order +0.04 

The hydrogen deuteride (HD) was prepared the action heavy water 
hydride according the method Wender, 
Friedel, and Orchin (1949) modified Fookson, Pomerantz, and Rich 
(1951). Mass spectrometric analyses several samples showed from 
95% HD. measurements with the brass tube was necessary prepare 
fresh frequently since and HDO lines gradually came and inter- 
fered with the measurements the lines. Preferential adsorption 
the metal wall and exchange reactions with the already present 
probably account for the appearance and HDO vapor the tube. 


OBSERVED SPECTRA 


each the 1-0, 2-0, and 3-0 bands seven lines were observed, the 
band six lines. The region the band freest from lines. 
spectrogram this band reproduced Fig. can seen, the lines 
form and branch with zero gap between, i.e. expected the 
band has the same structure that ordinary infrared band hetero- 
nuclear diatomic molecule. Fig. part recording the band 
shown including the R(1), R(2), and R(3) lines well number 
lines. these regions where many lines are present the lines 
can only recognized comparison with spectra taken with the 
absorption tube. Even some the observed lines are blends with 
and possibly HDO lines. 

Table the results the measurements the four bands are listed. 
For the 3-0 band the mean the measurements from the photographic 
infrared plates and from the infrared recordings given. These measurements 
are given separately Table order show the slight systematic difference 
between them. 

The sharpness the absorption lines all plates and recordings even 
pressures atm indicates that the slit width used was the limiting factor. 


4-0 band HD. 
Absorbing path 1000 atm, First order 21-ft grating. 
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TABLE 

Wave numbers lines observed bands 

1-0 band 2-0 band 
R(J) o-c o-c R(J) o-c P(J) o-c 
3874.374* +0.009 3450.440 —0.023 7306.541* +0.010 
4009.096 +0.002 

3-0 band band 


7 *Blends of doubtful accuracy. 


10445.466 +0.005 13551.066 


TABLE 


much the same way for the spectra HCI, HF, CO, and others except 
that here, because the high value, the branches are very short and therefore 
one place another the standard method must slightly modified. 

the line positions are calculated from the preliminary constants based 


Line Infrared Phot. infrared Average 
P(3) 10068 408 10068 
P(2) 178.400 178.382 
511.455 511.468 511.462 

R(2) 564.965 
Any pressure broadening the lines below the detectable limit. This 
result similar that for the quadrupole spectrum for which 
broadening was detectable atm (Herzberg Any collision induced 
dipole and quadrupole moments and are clearly extremely small. 

BAND ANALYSIS AND MOLECULAR CONSTANTS 
already mentioned, the observed band structure agrees with the expecta- 

tion for dipole radiation. Because the absence the distinction between 
symmetric and antisymmetric rotational levels, transitions with 
can occur and fact are the only ones occurring, just for ordinary hetero- 
nuclear molecules. Consequently, the analysis the band structure proceeds 
i 


the 3-0 and 4-0 bands and the Raman data Teal and MacWood (1935) 
systematic deviations from the present data (—0.7 for the and —0.3 


Comparison infrared and photographic infrared 
measurements 3-0 band (in 
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for the 2-0 band) are found. The present data allow much more precise 
evaluation themselves and can also compared and combined with the 
more recent accurate Raman data Stoicheff (1957). 

Since all the branches are very short not possible obtain reliable 
value for the rotational constant have therefore adopted for all 
vibrational levels the value obtained Stoicheff for 
the level from the rotational Raman spectrum. 

Unfortunately, only two values (for and are available. 
Subtracting from these values they give 


for and Therefore and can just determined. For this 
purpose the average for all four bands were used. The and values 
obtained the same way from the are much less accurate since only one 
set available for each and averaging possible. Therefore for the 
upper states, order take account all observed lines, the constants 
and were determined directly from the observed line frequencies adding 
obtained from the constants just determined and making least 
squares solution for vo, B’, and D’. The results are presented Tables III 
and IV. Table III, for comparison, the rotational constants for and 


TABLE III 


Rotational constants (in 


Present work 

0.0244, 

0.0231, 

TABLE 


Band origins and vibrational quanta (in 


obtained Stoicheff from Raman spectra are given. seen that the agree- 


ment most satisfactory. 


final check the wave numbers the band lines were calculated back 
from the constants. The deviations (o—c) obtained this way are included 
Table These are well within +0.02 for non-blended lines the 1-0, 2-0, 
and 4-0 bands. The only larger deviation the case the 3-0 band probably 
due unrecognized blending with lines the band 9400 
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Figs. and the values and less linear term are plotted 
against Just the case (Herzberg and Howe 1959) clearly seen 
that the first two points are appreciably above the (straight) line representing 
the remaining points. account this effect representation the 


0.00 
-0.10 
-0.20 


-0.30 


Fic. curve the ground state HD. 
(a) Plotted directly. Plotted after subtraction linear term. 


values formula somewhat ambiguous. three-parameter formula 
through the first three values 


But yields deviations (o—c) for and +0.0427 and +0.0446 
while four-parameter formula using the first four values 


45 
= 
1 2 3 4 Vv 
cm 
0.20 
By- [45.50 -1.80 
0.10 
~ 
3 4 Vv 
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all five values are used for five-parameter solution one finds 

while least squares solution with four parameters yields 


with residuals —0.0018, +0.0072, —0.0108, +0.0072, —0.0018. The values 
following from the last two equations are rather different, showing, did 
the work that the representation the curve polynomial 
not very satisfactory. seems clear from the shape the curve Fig. 
that the curve concave upwards near and therefore likely 
lead value that higher than that which obtained linear extra- 
polation from and that is, higher than 45.6317 This would 
eliminate the value from (3) and suggests that the from i.e. 


45.663 


close the true value although must realized that the other values 


represent the observed values just well even better. 
The values are plotted against Fig. Their accuracy not sufficient 


distinguish the resulting curve from straight line 


2 3 


Fic. curve the ground state HD. 


probable that for higher values the curve has minimum with 
subsequent steep rise does the corresponding curve for 

The band origins determined from the least squares solution for each 
band are given Table together with the first differences, i.e. the 
values. The first value, may compared with the value obtained 
Stoicheff from the Raman band: 3632.05; Our value higher 
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TABLE 


Rotational and vibrational constants the 
ground state 


6 Stoicheff gives 3632.053. 


0.096 When this discrepancy was first found amounted 0.125 
Subsequently the whole process determining the wavelengths was 
carefully scrutinized but error was detected the infrared measurements. 
Fortunately number lines are present the recordings the 
band which have recently been measured independently under high resolution 
Plyler and Tidwell (1957). When these lines were measured our 
recordings, wave numbers were found for them which differed only 
+0.029+0.010 from those Plyler and Tidwell. Even though the 
absolute accuracy Plyler and Tidwell’s data not known have corrected 
our data the 1-0 band the amount given. Such correction could 
ascribed the breakdown the integral relation between overlapping orders. 
The remaining discrepancy with Stoicheff still well outside the probable 
error each determination. independent either our 
least squares solution. This can seen the following combination relation 


which must hold rigorously between the lines the Raman band and the 
and lines the infrared band can readily verified. However, one 
finds for land for the left-hand side from the Raman data 7260.252 
and 7248.708 and for the right-hand side from the infrared data 7260.461 and 
7248.911 that is, for both cases the left-hand side 0.20 lower than 
the right-hand side, which means difference 0.10 for vo. The reason for 
this discrepancy not clear. The following vibrational formulae are based 
the infrared measurements. the Raman value for should turn out 
the correct one, will have reduced approximately 0.1 

term. While the four points are not sufficient establish the shape the 
curve higher values, the region investigated there again strong 
similarity H», that the first two values are appreciably above the 
straight line which represents the next values. For the same reasons for 
again difficult find unambiguous vibrational constants. 

three-parameter formula representing the first three values 
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° 2 3 4 


Fic. curve the ground state HD. 


outside any possible error. four-parameter formula derived from all four 
observed values 


Finally three-parameter expression obtained the least squares method 
from the four-observed values 


leading residuals +0.034, —0.101, +0.101, and —0.034 The dif- 
ferences between these formulae point the degree ambiguity the 
when higher values become known. may removed only when the 
initial curvature understood from theory. must emphasized, however, 
that except for the very small discrepancy between Raman and infrared 
values there ambiguity the relative position the energy 
levels. The ambiguity exists only with respect their absolute position with 
regard the potential minimum (i.e. the value the zero-point energy). 
Transforming equation (9) the vibrational energy levels obtain 


Table the rotational and vibrational constants the ground 
state are collected. The value has been obtained from neglecting the 
mass the two electrons while has been obtained from corrected B,* 
the way indicated Stoicheff (1957), and may compared with his 
results. These two values may compared with the corresponding values 
for which are, according Herzberg and Howe (1959), 0.74158, 
0.74116 gratifying note that the two values agree even 
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better than the the differences being well within the accuracy the deter- 
minations. The corrections included going from B,* are such that 
they should make r,* identical for isotopic molecules. Stoicheff (1957) the 
basis less complete data found greater difference between the than 
between the The improvement here obtained due the availability 
data for higher vibrational levels and the recognition the anomaly the 


INTENSITIES 


Wick (1935) has predicted approximate value for the absolute intensity 
the fundamental the infrared. From this figure path atm 
has been estimated the minimum observe this band. The actual path 
used here for the observation seven lines the band was 200 atm. 
The strongest line R(1) was also observed atm and had then inten- 
sity which suggested that the minimum was probably atm good 
qualitative agreement with Wick’s prediction. 

striking fact that, although the maximum path used was only 
times the minimum required for the band, overtone bands were 
observable. other words the intensity the overtone bands falls off far 
more slowly than ordinary heteronuclear diatomic molecules like 
HF, CO. CO, for example, the 4-0 band requires path about 
100 atm (Herzberg and Rao 1949) while the band can obtained 
with path less than 0.1 atm, that is, the intensity ratio greater than 
100,000:1. 

have attempted rough determination the relative intensities 
three the bands 1-0, 2-0, 3-0 measuring the integrated intensities 
selected lines these bands, selected the basis intensity and freedom 
from blends. These were R(3) and R(4) the band, P(1) and P(2) the 
2-0 band, and R(1) the band. Assuming that for these low intensities 
the equivalent width (i.e. the integrated intensity) proportional the true 
intensity and transforming means the theoretical intensity ratios 
R(1) each band one should obtain the same ratio independent whether 
R(3) R(4) and P(1) P(2) 2-0 are used. The considerable scatter 
that actually found (up may due either errors placing the 
base line 100% absorption failure proportionality equivalent 
width and true intensity. average ratio find, adding the estimate the 
intensity ratio the 3-0 and bands given Herzberg 


= ] 0.24 0.08. 


view the large scatter just mentioned not certain whether not 
the initial rise the intensity going from the the band real. 
any rate the numbers given clearly demonstrate that the intensity decrease 
anomalously slow. Qualitatively this anomaly can understood the 
fact that HD, unlike ordinary molecules like CO, the dipole moment 
produced only the vibration and increases with it. the equilibrium 
position negligibly small. illustrate this point Fig. shows four phases 
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Motion center positive (and negative) charges during vibration HD. 


the vibration HD: The center mass course fixed; the center 
the positive charges halfway between the and nuclei (indicated 
and therefore oscillates with the frequency vibration and with amplitude 
that increases with the amplitude the vibration. the electrons were 
following instantaneously the motion the nuclei their center would coincide 
with that the positive charges, but actual fact the electrons lag slightly 
behind the nuclei and therefore dipole moment arises. The magnitude 
this moment increases with the amplitude vibration since the velocity the 
nuclei increases and therefore the lagging-behind the electrons becomes 
more pronounced. quantum theoretical formulation was given few years 
ago (1952). According him the intensity ratio depends constant 
whose square essentially the ratio intensity Raman and Rayleigh 
scattering. Using the observed ratio for that 0.065 one finds from 
diagram 


4—0 


while for 0.2 one finds 


Within the rather large uncertainty. the experimental intensities and the 
approximations made the calculations the agreement must considered 
satisfactory. 

may perhaps significant that one matches the observed 
suitable choice then both and are larger than predicted. spite 
good deal searching the band was not found indicating that unlike 
the intensity 5-0 not high compared the predicted value. The 
intensity excess the even-numbered bands may due large quadratic 
term the development the dipole moment function displacement. 


CONCLUSION 


The observation electric dipole infrared spectrum establishes 
qualitative distinction from which cannot have electric dipole spectrum 
any approximation, and this distinction course conformity with the 
requirements quantum mechanics. The intensity this spectrum deter- 
mined the mass difference. For the intensity would appreciably 


Oh 
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greater, roughly factor (3/2)?. may noted that just 
would also have quadrupole absorption spectrum which would have 
intensity very similar that that is, path lengths times those 
employed had been available the infrared bands would have shown addition 
the and branches also and branches whose intensities would 
show different dependence than the and branches. 

spite the ambiguities the constants mentioned believe that they, 
together with the Raman data, represent the best values present available. 
The energy levels the level are now known with considerable 
precision and can used for precise calculations thermodynamic functions. 
However, the and AG, formulae cannot safely extrapolated higher 
values because was found for formula with three four parameters 
will represent the whole course the and curves. order obtain 
reliable values for the higher vibrational levels analysis the 
vacuum ultraviolet spectrum under high resolution will required. 
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SOME RECENT DETERMINATIONS ATOMIC MASSES 
THE STRONTIUM-ZIRCONIUM 


ABSTRACT 


large double-focusing mass spectrometer has been used obtain new values 
for the masses and Mass differences calculated from these values 
are found better agreement with nuclear transmutation information than 
were previous mass spectroscopically derived values. 


INTRODUCTION 

There has existed the past considerable disagreement between mass 
spectroscopic and nuclear transmutation determinations the masses 
certain nuclides. The most serious these, that the secondary standard 
atomic mass, has been resolved only the last few years. Another 
major discrepancy between the data from these two fields investigation 
has existed the region (Duckworth 1957). 
investigation large double-focusing mass spectrometer, recently completed 
McMaster University, was used attempt resolve this discrepancy. 


EXPERIMENTAL 

The Instrument 

The mass spectrometer used double-focusing instrument, similar 
many respects the type first used Dempster (1935). Figure 
schematic diagram the experimental arrangement. are formed 
electron bombardment source which was modified provide for the evapora- 
tion reasonably volatile solids, and are accelerated through potential 
difference kilovolts between and The angular spread the beam 
entering the electrostatic analyzer limited the slit which 0.12 in. 
width. The 90° electrostatic analyzer cylindrical condenser and 
employed symmetrically. this respect the arrangement differs from Demp- 
ster’s original arrangement. The energy spread the ions entering the magnetic 
analyzer limited .05-in. slit, Ions one particular value specific 
charge are deflected the magnetic field that they pass through the 
collector slit, Ss. Those that pass through this slit emerge from the magnetic 
field and are detected 14-stage Allen-type electron multiplier. 

The radius curvature the ion path both magnetic and electrostatic 
analyzers ft. With 0.002 in. and 0.003 in., the instrument has 
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ELECTROSTATIC 


COLLECTOR 


MAGNETIC 
ANALYZER 


Fic. Schematic diagram the mass spectrometer. 


theoretical resolution (at the base the peaks) about part 20,000. 
The experimental resolution obtained this investigation was about part 
18,000. 

The magnetic analyzer was used previously single-focusing instrument 
and its features described (Kerr, Bainbridge, Dewdney, and Duckworth 
1959). When the electrostatic analyzer was added, the positions the source 
and electron multiplier were changed the positions shown Fig. and the 
current through the magnet coils was reversed. 

Each plate the electrostatic analyzer made blocks, each about 
cross section, fitted together and machined give con- 
tinuously inner surface. The blocks are made Armco iron, are gold- 
plated, and sit alumina insulators, ground height 0.625 in., which, 
turn, rest thick steel base plate. The gap between the condenser plates 
in. Fringe blocks ground potential, located 0.255 in. from the ends 
the condenser plates, are used terminate the electrostatic field (Herzog 
1935). The steady deflecting voltage the electrostatic analyzer plates 
obtained from batteries. 
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The Peak-Matching Technique 

The method obtaining mass difference values variation the now 
well-known peak-matching technique introduced Smith and Damm (1956), 
and adapted for use deflection instrument Quisenberry, Scolman, and 
Nier (1956). 

wave form consisting 30-cycle square wave (derived from batteries 
and chopper relay) plus 60-cycle sawtooth (from oscilloscope) imposed 
the steady deflecting voltage the electrostatic analyzer. When the 
amplitude the square wave properly adjusted, the two members the 
mass spectral doublet under study are alternately transmitted the instru- 
ment. The sawtooth component serves sweep voltage. this way the 
doublet members are alternately swept across the collector slit and displayed 
peaks the screen oscilloscope (the same instrument which 
supplying the 60-cycle sweep mentioned above), connected the output 
the electron multiplier a-c. amplifier. The square wave amplitude 
adjusted carefully until the peaks appearing the screen are superimposed. 
The mass difference between the doublet members then given the 
expression 


V+AV 
where the amplitude (peak peak) the square wave, the total 
steady voltage across the electrostatic analyzer, and the mass the 
heavier member the doublet. The voltages and are measured with 
differential voltmeter. 

Many voltage measurements were made each doublet, and probable 
error was arrived combining the usual statistical error with the specified 
error associated with the voltmeter. 


RESULTS AND DISCUSSION 


Table summarizes the mass differences obtained during this investigation. 
Earlier related results are included for comparison. The doublet 


TABLE 


New mass differences (in milli-mass units) with 
values obtained previous work 


This work 


Doublet 


Previous work 
11.2390+30 11.244 +10 (a) 

11.2355 +6* 

(c) 
100.126 +25 (d) 
48.240 +12 previous value 
36.999 37.00 +18 (e) 


4.9484+25 previous value 


*Actual doublet observed was N:-CO. 
tActual doublet observed was CeHi4-Sr™. 
(a) Kerr and Duckworth (1958). 


(b) Scolman, Quisenberry, and Nier (1956). 
(c) Smith (1958). 

(d) Collins, Johnson, and Nier (1954). 

(e) Duckworth and Preston (1951). 
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was used check the performance the instrument during the course 
the investigation. The and ions were obtained evaporating 
strontium metal and zirconium tetrachloride into the ionization region the 


source. 
Table gives the new mass values directly calculated from these differences 


and compares them with previous work. 


TABLE 


New atomic masses (in atomic mass units) with 
values obtained previous work 


Nuclide This work Previous work 
14.00752658+9 
85.93667 85.93648 +10 (c) 
(d) 
(d) 
89.93350 +20 89.93292 +25 (c) 


89.93295 +36 (e) 


(a) Scolman, Quisenberry, and Nier (1956). 

(6) Smith (1958). 

(c) Collins, Johnson, and Nier (1954). 

(d) Duckworth and Preston (1951). 

(e) Duckworth, Preston, and Woodcock (1950). 


The mass values used determining these new values were: 1.0081451 
and 12.0038156+4 (Quisenberry, Giese, and Benson 1957), 
34.9799720+22 (Giese and Benson 1957), and 
(Duckworth 1957). Evidently, the accuracy the new mass value 
largely limited the error the comparison mass 

Comparison with Mass Differences from Nuclear Reaction Data 

The Q-values the reactions and have been 
published (Way al. 1955). These results are compared with mass spectroscopic 
data Table 


TABLE III 


mass difference 


Transmutation data 1.99698+6 (Way al. 1955) 
Mass spectroscopic Johnson, 
data and Nier 1954) 
(Duckworth 
and Preston 1951) 
1.99696+5 (This work) 


The results this investigation agree very well with the transmutation 
data, while the best previous mass spectroscopic data disagree about 
0.5 milli-mass unit. 

Two groups values are obtained from transmutation data for the 
mass difference, depending upon whether the route the 
route taken. The former gives weighted mean 
1.99932+10 atomic mass units, while the latter gives The best 


— 
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previous mass spectroscopic data (Collins, Johnson, and Nier 1954) give 
1.99896+27 atomic mass units, while this work gives 

discrepancy 0.39 milli-mass unit the closed loop 
lies the base the discrepancy the transmutation 
data. Our work supports the link. The reaction 
(Chidley, Katz, and Kowalski 1958) seems clear-cut case where the 
ground state was definitely observed. The Q-value the decay 
was reported (Peacock and Jones 1948) 3.97 milli-mass units. This 
value has since been estimated too high (Stirling and Goldberg 1956) 
0.27 milli-mass unit. Actually, the old Q-value better agreement with 
the other transmutation data and also with the mass values obtained this 
investigation. This discrepancy transmutation data deserves further 
attention. 


REFERENCES 


G., Katz, L., and 1958. Can. Phys. 36, 407. 

L., H., and 1954. Phys. Rev. 94, 398. 

1935. Proc. Am. Phil. Soc. 75, 755. 

1957. Progress nuclear physics, Vol. (Pergamon Press Ltd., 
London), pp. 138-161. 

and Preston, R.S. 1951. Phys. Rev. 82, 468. 

E., Preston, S., and 1950. Phys. Rev. 79, 188. 

and Benson, 1957. Bull. Am. Phys. Soc. II, 223. 

1935. Physik, 97, 596. 

T., BAINBRIDGE, R., W., and H.E. 1959. Advances 
mass spectrometry (Pergamon Press Ltd., London). 

Kerr, and 1958. Can. Phys. 36, 986. 

Peacock, and 1948. U.S. Atomic Energy Commission AECD-1812. 

QUISENBERRY, S., F., and 1957. Bull. Am. Phys. Soc. II, 
223. 

S., T., and 1956. Phys. Rev. 102, 1071. 

T., QUISENBERRY, S., and 1956. Phys. Rev. 102, 1076. 

1958. Phys. Rev. 111, 1606. 

and Damm, 1956. Rev. Sci. Instr. 27, 638. 

STIRLING, and 1956. Bull. Am. Phys. Soc. II, 291. 

Way, K., L., and vAN 1955. U.S. Atomic Energy 
Commission 


es 
- 
1S 
yn 
-88 


ABSOLUTE MEASUREMENT THE ACCELERATION 
DUE GRAVITY OTTAWA! 


ABSTRACT 


apparatus for determining the absolute value gravity measuring the 
distances through which rule falls discrete time intervals described. From 
the data associated with drops with two non-magnetic stainless steel rules 
vacuum, value the absolute gravity station Ottawa 980.6132 
with possible error +0.0015 has been obtained. This value 
13.7+2.0 milligal less than the Potsdam value that position. 


INTRODUCTION 


The determination the acceleration due gravity (g) matter con- 
siderable importance since its value involved the determination 
number derived physical, geophysical, and astronomical quantities. The 
value both position and, slight extent, time dependent; however, 
the variations involved can measured very accurately (to few parts 
10° that principle only single absolute determination required 
for the correct values assigned the entire gravity grid the earth. 
Reversible pendulums were used make series such absolute determina- 
tions Potsdam from 1898 1904 (Kiihnen and 1906) and 
from these are derived the ‘Potsdam values’ now used throughout the 
world. 

number later determinations have been made various places. 
These have demonstrated that the Potsdam system gives values that are 
too high amount that between and milligal (~8 and 
but the discrepancies among these measurements are such that they cannot 
provide acceptable alternative the Potsdam system, i.e. one that would 

the 1954 Rome meeting the International Union Geodesy and 
Geophysics was requested that the standardizing laboratories the 
participating countries attempt complete new absolute determinations 
which would enable the Potsdam system superseded the 1957 Toronto 
meeting the I.U.G.G. fact, sufficient data were not available the 1957 
meeting, and the request was renewed and emphasized with the 1960 Helsinki 
meeting the being proposed for consideration the results. The 
measurements reported here constitute part this co-operative effort. The 
methods, results, and errors are discussed the following sections. 


received April 1960. 
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PREVIOUSLY PUBLISHED RESULTS 


Although several basically different methods determining are available, 
practice only two them have been are prospect being used: these 
are the determination the periods reversible compound pendulum and 
the timing body free fall. Before the Second World War, inaccuracies 
the measurement short intervals time were such that the latter method 
could not compete with the former, and only compound pendulum determina- 
tions were made. From the large number these, some few, commonly 
conceded, stand out reason both the care with which the determinations 
were made and the known corrections computed, and the experience the 
experimenters. Among these are the experiment referred above Kiihnen 
and Potsdam the period 1898 1904, the experiment 
Heyl and Cook (1936) Washington, and the experiment Clark (1939) 
Teddington. The results these experiments were reworked Jeffreys (1948) 
give —10.7, —17.9, and —13.1 milligal the Potsdam system respectively. 
Postwar pendulum determinations Agaletzki (1956) Leningrad and 
Baglietto (1957) Buenos Aires gave results —12.1 milligal and —7.1 
milligal from Potsdam respectively. 

1946 Volet suggested that the timing rod free fall could principle 
give results that were comparable with those obtained from compound 
pendulums far accuracy was concerned, and that such measurements 
would not subject the systematic experimental errors that affected the 
pendulums. These errors may conceivably quite high (of the order 
milligals) and may have affected all pendulum experiments the same sense 
and approximately the same amount. contrast seems probable that 
the systematic errors free fall experiments will random sense from 
experiment experiment. Free fall experiments, other than preliminary 
and non-definitive one Volet (1952), have been made Agaletzki (1956) 
and Martsiniak (1956) Leningrad two slightly different methods giving 
results —7.5 and —9.3 milligal the Potsdam system, and Thulin 
(1958) Sévres giving —12.8 milligal the Potsdam system. addition, 
free fall experiments are underway Brunswick, Teddington, and Washington. 

must pointed out that the reduction all these results the 
Potsdam system introduces the present errors the gravity ties between the 
experimental stations concerned; these errors may some instances the 
order milligal. However, undoubtedly possible reduce these com- 
parison errors the order half milligal less. 


FREE FALL METHODS 


One the objects this determination was the avoidance any 
systematic errors that might exist previous determinations means 
compound pendulums. The only method measurement that did not appear 
though would probably introduce more systematic errors than removed 
was the determination the motion body free fall. was also necessary, 
satisfactory substitute for the Potsdam system was obtained, that 
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The simplest free fall experiment consists determining the time fall 
(t) and the corresponding distance fallen body that has been released 
from rest. The local value then, course, given by: 


(1) 


and see that and are the errors the time and distance measurements 
respectively that the error given by: 


(2) 


For purposes comparison will assume that meters (giving 
second), that lengths are measured with uncertainty part 
0.25 meter, whichever the greater, while time intervals are measured 
with uncertainty second. From equation (2) then obtain 
maximum possible value 0.26 milligal, 0.2 milligal arising from 6S/S, 
and 0.06 milligal from 26t/t. This satisfactorily low error, but unfortunately 
the experiment cannot performed this manner. impracticable 
release body from rest without imparting initial transitory acceleration 
it, and impossible prevent initial relaxation the body after its 
release. Position and time measurements cannot, therefore, commenced 
until some time after the release. 

practice, this leaves with two alternative methods measurement, 
involving either falling object rising and falling object. the former 
method, adopted this laboratory, two time intervals (¢; and and the two 
corresponding length intervals and all measured from arbitrary 
after the body released from rest. then given the formula: 


If, corresponding the previous case, consider the uncertainties 


retain the same accuracy measurement and length fall before, 
(and hence and this should done such manner that the maximum 
value minimized. The method doing this shown Appendix 
and from this and other considerations value was selected for 
giving 0.25 micron and 0.2 second. then have maximum value 
when and all have their maximum numerical values and 
signs such that the errors they introduce are additive, i.e., —0.25 micron, 
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Inserting the appropriate values S2, and then obtain: 


1.2 milligal. 


Thus can seen that although there decrease the accuracies 
measurement the time and distance intervals, yet this method has less than 
one quarter the precision the but impracticable, first case. is, 
however, sufficiently accurate for our purposes, even allowing for the effects 
certain additional errors. 

For completeness should pointed out that throwing the free-falling 
body vertically upwards and timing both rise and fall (i.e. using the second 
the alternative methods referred above), possible obtain the precision 
the method. fixed length this case and and are the 
time intervals between the two crossings the lower boundary and between 
the two crossings the upper boundary respectively have: 


(6) 


and 


(7b) for t, > te. 


For the same value and accuracy measurement before and the 
additional (and obvious) requirement that made very small obtain, 
the case, maximum value 0.26 milligal. Thus increase 
accuracy factor 4.6 (for comparable geometry) will result from the 
substitution the ‘‘rise and for the ‘‘free experiment. However, 
the additional complexity the “rise and experiment considerable, 
and value high meters for would probably impractical such 
appraisal the relative accuracies. consideration this and the relative 
complexities the experiments resulted our selecting the 
system. “rise and experiment currently underway Teddington. 

For the equations presented above was assumed that the local value 
was constant over the entire experimental region; general not. Above 
the surface spherically symmetrical sphere radius the variation 
given by: 


2dr 


(8) 


- 
(7a) 
9 
( 
CO 
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which corresponds gradient about —0.3 milligal for the case 
the earth. Below the surface such sphere having surface density 
have the expression: 


which for the earth gives dg/dr varying linearly with from —0.23 milligal 

can seen, therefore, that the variation will usually significant 
over vertical distance meters. However, such variation present 
approximately linear over the experimental apparatus (about meters 
the present case), the appropriate correction for can readily calculated 
(see Appendix I). 


APPARATUS 


indicated above, the value was determined this laboratory 
measuring two vertical intervals length and the two corresponding time 
intervals for rule falling freely through these vertical intervals. 

choosing the type rule dropped there are, broadly speaking, the 
alternatives of: standard meter rule carrying only two end plates, the 
upper one for suspension purposes, the lower one make contact with the 
deceleration assembly, more complicated assembly. 

Alternative (a) has the advantages employing rule ‘primary standard’ 
quality and being rigid construction that movements the graduations 
about the center gravity are virtually eliminated after the initial short 
delay period. is, however, limited ultimate accuracy its meter 
usable length; also difficult achieve adequate uniformity temperature 
along its length, while the difficulties verifying beyond doubt that this 
required uniformity exists are, owing the insulating properties vacuum 
and metallic point contacts, almost insuperable. For these reasons 
selected the second alternative (b); this enabled use 2-meter rule, 
provide radiation shield, measure temperatures the rule, measure 
tilt the rule during the drop, and provide adjustments for obtaining 
truly vertical hanging position. The resulting complex construction con- 
siderable disadvantage, but was probably justified. 

Each the two rules used had mounted three semitransparent scales. 
These were spaced that, the rule fell, the scales passed the axis 
camera synchronization with short-duration light source activated 
times per second. set fiducial cross hairs was incorporated the camera 
system. The intervals between the scales were calibrated terms the 
International Metre, and the time interval between the flashes was found 
photographing oscilloscope traces relating the light flashes standard 
frequency signal. appropriate interval before the first flash, electrical 
impulse released the rule from its topmost position (see Fig. 1); was 
mechanically arrested the end its fall. motor-driven carriage and reset 
mechanism then returned the rule its original position without the necessity 
breaking the vacuum seals. 
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résumé the main points considered the construction the apparatus 
follows: 

(a) Semitransparent glass scales, ruled intervals, were attached 
stainless steel rule. The relative positions the scales, the position each 
graduation the scales, and the over-all coefficient thermal expansion 
the rule were determined. 

(b) reduce the effect calibration errors, five stationary fiducial lines 
were used each determination position. The timing the first flash, and 
hence the rule positions, was varied from drop drop, and two sets (Marks 
and rules and scales, independently calibrated, were used the 
experiment. 

(c) Because the high coefficient thermal expansion 
the rule, was desirable know the temperature along the length the rule 
within approximately 0.005° This was accomplished standard 
temperature stabilizing and thermometry techniques; temperatures were auto- 
matically recorded. 

(d) Although the center gravity the rule assembly has acceleration 
other forces being absent, the scales can move slightly with respect the 
center gravity and lead erroneous determinations The rule was 
assembled reduce this relative movement much possible and 
provide high frequency oscillation and high degree damping the 
oscillations set the initial shock release. 

(e) The rule was dropped maintain nearly vertical was 
practicable throughout the drop. This served reduce focussing difficulties, 
cosine corrections length, and indeterminacies arising from the non- 
horizontality the rule the optical axis the camera. Since was not 
possible maintain the rule exactly vertical, the optical axis this point was 
lined close the horizontal was practicable. Determination the 
angles the fiducial lines with the horizontal were also made. 

(f) Non-gravitational forces affecting the rule are those due residual gases 
the column and the interaction the rule with magnetic fields; electro- 
static forces are certainly negligible. The pressure was the order 
measured ionization gauge. This was high enough impart 
appreciable deceleration the rule assembly and correction for this effect 
was determined experimentally. The entire rule assembly was made 
nominally non-magnetic materials that the only effect magnetic field 
would that arising from the introduction eddy currents from some 
ferromagnetic inclusions the rule assembly. approximate calculation 
the upper magnitude these effects due the original field and field gradient 
suggests that they are negligible. Nevertheless, the column was wound with 
series current carrying coils which reduced the vertical component 
the column. 

(g) Three seismographs were used monitor the movements the fixed 
camera mounting. This ensured that the position this was not affected 
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random microseisms or, more seriously, since this would constitute systematic 
instead random error, the release shock transmitted through the 
apparatus and the building walls. 

(h) The time intervals were measured determining the departure from 
the nominal 0.1-second repetition rate the illuminating flash. This departure 
did not, fact, exceed second. 

(t) the rule position and the timing are accurately and unambigu- 
ously determined, necessary have high-contrast, fine grain photo- 
graphic plate and short period flash with relation that 
approximately symmetrical about the peak intensity. Two types plates and 
two light sources (air spark and xenon flash) were used this experiment; 
these appeared approximate sufficiently closely the required conditions. 

Details the more important these considerations are follows. 


Construction and Calibration the Rules 

The two rules employed were machined from No. 304 (S.A.E.) stainless 
steel bar stock, and differed only minor details. Each rule (Fig. 
approximately meters long with channel section about 2.3 slotted 
accommodate five scales. 

The scales are aluminized give 30% light transmission, the graduations 
being ruled through the aluminum. Four these scales have range from 
mm, with 0.1-mm subdivisions; two longitudinal parallel lines 0.2 
apart define the calibrated portion the graduations while two diagonal lines 
serve identify the particular graduations which are being photographed. 
The scales are cemented stainless steel holders which, turn, were fixed 
the channel bar. They were aligned, using the standard optical techniques, 
colinear and coplanar (to within when the rule was unconstrained. 

The fifth scale auxiliary focussing scale (Fig. 2). has lines ruled 
0.1-mm intervals and has its plane set angle about degrees the 
plane the other four scales, its graduations being parallel their center line. 
When the delay circuit was adjusted that this scale was photographed 
during drop, were able measure the fore and aft deviation the rule 
during that drop. The accuracy this measurement depended the depth 
focus the camera well the horizontal separation the graduations; 
the system used was the order +40 

(Westinghouse type 14A, nominal resistance 10° 
5000 were held close thermal contact with the rule inserting 
them into aluminum holders which were themselves cemented and clamped 
the rule. The connecting wires (40 S.W.G. Tensolite cable) were cemented 
the rule along its length point where they were brought out the 
duralumin sheath, cemented it, and attached the 12-pole contactor 
shown 

The radiation shield duralumin tube diameter, 210 long) 
which clamped the collar supporting the rule its upper Airy point. 
the lower Airy point screws project within about 0.1 the rule when 
the assembly the vertical position. These limit whip the rule when 
strikes the arrestment and provide the conventional second support for the 
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PRESTON-THOMAS AL.: MEASUREMENT ACCELERATION 833 


rule when being calibrated the horizontal position. The assembly 
made hang vertically (from radius hemisphere resting plane) 
means the coarse and fine adjustments shown Fig. 

Calibration 

The intervals between the calibrated scales the rules were indirectly 
compared with the Canadian Standard Metre rule No. 306, the transfers 
being effected means two laboratory reference standards Nos. 169 and 
7556, and meters length respectively. Rule 306 has been compared with 
the standards the International Bureau Weights and Measures several 
different occasions since its acquisition 1921, the latest calibration being 
May 1957. 

The two reference standards are invar. Rule 169 was purchased from the 
Société Genevoise Physique 1913; 1951 was sent back 
the makers for refinishing and regraduating, and was received back the 
National Research Council 1953. Rule 7556, obtained expressly for the 
purpose calibrating the rules, was bought from the Société Genevoise 
1956. both rules the graduations are narrow and cleanly cut—excellent for 
pointing operations. 

The calibrations the two principal intervals each rule were carried 
out immediately before and after each definitive series drops. 

Since was impractical immerse the rules directly the water the 
comparator tank, was necessary encase them specially constructed 
watertight (copper) tubes, with dams around the viewing ports prevent the 
entrance the water. order that the reference rules viewed under the 
same conditions, i.e. through air rather than through water, similar enclosures 
had provided for them. 

The subdivisions (0.1 mm) the scales had been previously calibrated 
terms the over-all length each individual scale and these lengths, together 
with certain other principal intervals the scales, had been compared situ 
(by the Abbé method) with master scale. The corresponding intervals the 
master scale had been measured interferometrically and had been compared 
with subdivisions rule 169. Thus the final calculation each determina- 
tion position could referred the zero the appropriate scale. 

The short (i.e. the cm) intervals the rules were compared with rule 
169. These measurements were carried out four steps the following 
sequence: 

(a) Four different intervals, covering practically the entire operating range 
scales and (Fig. were measured using various intervals the region 
adjacent the zero end rule 169. 

(b) Rule 169 was then turned end-for-end and the procedure repeated 
using, this time, intervals the region adjacent its 100-cm end. 

(c) and (d) These calibrations were carried out after the series drops had 
been completed; the rule was turned end-for-end from its previous position 
and the procedures and repeated. 

The basic short interval used the subsequent calculations was the 
mean these measurements corrected the respective zeros their 
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scales. Table gives the results the calibrations the short interval 
the rule; the standard deviation the determinations 0.2 


TABLE 
Comparison the calibrations various lengths the short (27 cm) interval 


Distance between scale (20° cm) 


Interval (cm) November 1959 December 1959 
0.70 0.88 26.417 138 26.417 138 
169 70.0 
0.62 1.30 26.417 150 26.417 150 
169 97.5 70.0 
1.68 26.417 135 26.417 109 
169 97.5 70.0 
0.40 0.08 26.417 142 26.417 161 
169 70.0 
0.70 0.88 26.417 156 26.417 167 
169 0.0 
0.62 1.30 26.417 170 26.417 156 
169 0.5 28.0 
1.00 1.68 26.417 125 26.417 134 
169 0.5 
0.40 0.08 26.417 131 26.417 151 
169 0.5 27.0 
Means 26.417 143 26.417 146 
Mean all 26.417 145 
Standard deviation 


the measurement the long (i.e. the 199 cm) intervals the rules were 
compared with rule 7556 using procedure similar that described for the 
short intervals. addition the long intervals were built terms rule 
169 using intermediate scale (mid scale, Fig. the rules. The differences 
between the direct and the built-up calibrations for these intervals were 0.1 


for each rule. 


Optical System 
The optical system used photograph the scales illustrated Fig. Its 


component parts are: 

(a) flashing light source very short duration which includes: xenon 
tube (or the air gap spark used for the rule) flashing exactly 
intervals; condenser system focus the light from the flash 
onto the and infrared filter prevent undue heating the rule when 
for focussing purposes; 

primary lens system, having very carefully levelled optical axis, 
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beams 


Fic. system. Legend: (a) xenon flash, shutter, (c) and (d) condenser lenses, 
heat-absorbing glass, (f) window, (g) scale, plane parallel viewing window, defining 
aperture, primary camera lens, fiducial cross scale image plane, (m) secondary 
camera lens, rotating photographic plate and image plane. 


which focuses image the scale into the plane the set five fiducial 


cross hairs; and 
secondary lens system which focuses the combined image scale and 


cross hairs onto the emulsion photographic plate which rotating such 
way that the pictures from successive flashes the drop are located 


circle without overlapping. 
The xenon discharge tube and its associated circuit (see Fig. are mounted 


shielded chassis just front the illuminating window. The tube itself 


470k 


10k 


pulse 


+250V 


tance, wh. 


has discharge gap about capillary which almost filled 
the electrodes that the scrubbing action the expanding gas will help 
keep the glass clean. xenon pressure about has been found give 
good relation with little afterglow. Use striking voltage 
kv) several times the breakdown voltage the tube produces short 


and accurately reproducible initiation times. 

The image the flash focussed into the plane the scale about 
high and wide. subtends cone light large enough ensure 
that the transmitted beam through any used point the scale completely 
fills the camera optical This necessary order that the effective 


() | A aul 
' 
beam 
flash 
470 
20k 


836 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


optical axis the camera may made and kept accurately horizontal, even 
though occasional slight changes are made the position the condenser 
lens obtain optimum illumination. 

The optical axis defined circular aperture located immediately 
front the first camera lens and with its center exactly the same level 
the point the plane the scale which focuses the middle one the five 
cross hairs. 

Aperture size was chosen This value gives depth sharp focus 
about +0.0040 and the pictures will therefore indicate immediately the 
rule not dropping straight. larger aperture makes difficult illuminate 
the system uniformly from the flash source and could result some loss 
sharpness the picture for tolerable horizontal shifts the rule. The 
remainder the optical system relatively straightforward. There 
over-all magnification about 10, equally divided between the two stages. 
The scale brought into sharp focus the plane the cross hairs moving 
the first camera lens and defining aperture horizontal slide (without 
rotation which might disturb the optical axis) after preliminary adjustment 
the rule itself moving the whole column. Lateral adjustments can 
made sufficient accuracy moving the column until the single vertical 
hairline approximately midway between the images the two vertical 
rulings the scale. have found that visual focussing most convenient 
the scale under continuous illumination through blue filter (so that the 
eye has the same color sensitivity the photographic plate) and the com- 
bined image the last focal plane viewed with small magnifying glass. 

Before any optical alignment procedures were carried out the mechanical 
axis the whole assembly was aligned and levelled. The optical axis finally 
determined is, fact, almost identical with the mechanical axis. The final 
adjustments for levelling the optical axis were made means precisely 
levelled surface plate carrying microscope. The microscope was focussed 
right bisector cross hairs accurately locating the center the aperture; 
was then translated along the surface plate and sighted onto the image the 
fiducial cross hairs the scale plane. this way was possible check the 
optical axis with all components place except the viewing window. The 
optical neutrality this window has been established separately. 

The lenses and cross hairs are mounted heavy channel iron beam which 
turn firmly mounted two J-beams. These are completely 
isolated from the main column that the disturbances caused release 
the rule will not result shifts camera position. The plateholder with its 
motor-drive mounted separate /-beam isolz from the rest the 
optical system. check for possible shifts two quick response seismographs 
are connected either end the camera mounting. further slow response 
seismograph used ensure that long term seismic disturbances can cause 
error the measurement. These seismographs showed that shifts greater 
than 0.1 occurred prior the arrestment shock during any drop. 
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Light Source and Photographic Image 

are concerned with recording photographic plate the position 
the image moving line that has been illuminated flash with some 
particular relationship. our case, three identical flashes are 
produced the instants and and require that the subsequently 
deduced distance intervals derived from the corresponding photographs shall 
give the correct value This will occur only the measured positions 
the photographic images depart from some arbitrarily defined position 
the optical images (e.g. that the instant peak intensity the flash) 
amounts 6x2, and 6x; such that (see equation (3)): 

(10) te (6x2—6x1) = (6x3—6x}). 

This condition satisfied Cv, where the velocity the rule 

may assume that photographic image produced the plate all 
points where the integral equal greater than certain 
value (K). That say image will then produced point distance 
optical image the peak intensity the flash) relation the form: 


(11) 

t+(d/20) 
satisfied, where the thickness the optical image, its velocity, and 


will observed that the limits development, giving the width the 
i.e. the flash symmetrical time. Such flash, therefore, satisfies the 
conditions equation (10) and has the additional advantage giving equal 
density gradients both limits the photographic image, thus compensating 
some extent for the non-infinite contrast the plates. This lack infinite 
contrast, even the best fine grain high contrast plates, makes necessary 
employ very short duration flash that the boundaries the photo- 
graphic image may unequivocally determined. short duration flash also 
serves prevent undue elongation the image. 

After considerable experimental work obtained air spark and 
xenon discharge approximating the desired conditions short duration 


PHOTOMULTIPLIER 


Fic. Oscilloscope trace showing 100-kc reference pulse, 10-Mc timing signal, and 
photomultiplier output giving time-intensity relation xenon flash. 


} 
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symmetrical light flashes (see Fig. 5). The departure from symmetry consists 
long duration tail which impossible eliminate entirely. 
can shown graphically that for the tails associated with our light sources 
there significant departure from the condition equation (10) unless the 
peak intensities the flashes illuminating the three scales differ 25% 
more. guard against this possibility, the light source was always started 
several seconds before the drop took place, and the photographic record the 
photomultiplier output showed that the intensities the three flashes 
illuminating any plate were, fact, substantially equal. 

Since the flash duration not infinitely short, there will difference 
the appearance the photographs the three scales due the differing 
velocities the rule. This can seen the microdensitometer traces 
scales and shown Fig. 


Temperature Control and Measurement 

The two stainless steel rules used for the final drops had coefficient 
certainties were not produce appreciable errors (0.1 milligal more) the 
mean temperature the 2-meter interval had known +0.004° and 
that the 27-cm interval +0.01° This fairly difficult requirement 
meet for bar that vacuum-insulated and has the dimensions and thermal 
conductivity the rules. was accomplished controlling the temperature 
the column, putting duralumin radiation shield around the rule, and limit- 
ing the heat input the rule. 

The column was ‘water-cooled’ maintain its temperature within 
+0.05° and the duralumin radiation shield reduced the rate drift 
temperature any point the rule +0.02° per hour for difference 
between shield and rule. The heat sources inside the column were three electric 
motors (for the damping mechanism, the contactor, and 
mechanism), which were located lose their heat the column, and the 
arrestment, which absorbed about 150 joules kinetic energy each drop. 
practice, this latter quantity heat appeared not reach the rule 
amount joules for 0.01° rise over length) sufficient 
troublesome. 

Besides these heat sources, certain amount heat was introduced the 
light source used during the damping the oscillations the rule and for the 
focussing the camera before each drop. This heat input was reduced far 
possible the use heat absorbing filters and rapid damping and 
focussing; however, the cumulative effect was sometimes such require 
respites from dropping during day’s operation. The gradual temperature 
drift evident Fig. from this cause. 

Temperature measurement was means six type 14A thermistors 
thermal contact with the rule. These were used conjunction with motor- 
driven selector switch, resistance bridge, and Leeds and Northrup micro- 
volt d-c. breaker amplifier feeding into AZAR recorder. These last two 
were adjusted give sensitivity 0.01° per chart division, see 
Fig. The entire unit, consisting thermistors, switch, bridge amplifier, and 
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= 


Fic. profiles along rule 


recorder, was periodically calibrated. During four such calibrations over 
period years calibration drift greater than 0.002° was observed. 
was thus possible determine the temperature profile along the rule with the 
necessary precision. 

selecting the appropriate mercury-in-glass thermostats was possible 
adjust the rule temperature any desired value both during the calibration 
and dropping the rule. The temperatures selected were such that the 
calibrations (excepting those used determine the coefficient thermal 
expansion, but including the calibration the Canadian Primary Standard 
Metre, No. 306, Sévres) and all drops were made with the three standard 
and the two rules entirely within the temperature interval 19.7° 
20.2° any one drop, calibration, the temperature profile the rule 
was within temperature range 0.08° and its rate drift temperature 
did not exceed 0.05° per hour. 

These precautions resulted all corrections for thermal expansion being 
very small and also made certain that errors could not arise from uncertain- 
ties the determination the coefficients thermal expansion. 

Time 

The accuracy time measurement determined by: 

(i) the accuracy measurement cathode ray tube traces terms the 
time scale defined Mc/s and 100 signal; 

(ii) the knowledge the frequency the crystal oscillator. 

measuring the timing the light flash against the 100 pulse all 
trigger errors are eliminated. The pulse shapes repeat themselves with such 
regularity that any arbitrary point the pulse can measured against the 
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Mc/s signal about 10% the period the latter, i.e. second. 
The time intervals and can therefore measured second 
terms the time scale defined the crystal oscillator. 
The absolute frequency this time scale, which forms part the national 
frequency standard, known parts through comparisons with: 
(i) the time scale U.T.2 determined the Dominion 
(ii) the caesium-beam frequency resonator standard; 
(iii) the time scale U.T.2 the other national observatories. 
The maximum total error time measurement therefore second 
part for the shorter time interval. 


Electronic 

The controls are designed release the rule preset interval time before 
the first photograph, open the camera shutter and, after preset number 
flashes, cut off the high voltage supply from the light source. additional 
function monitor the exact timing and shape the light pulse. All the 
control and monitoring signals are obtained from precision Essen ring type 
quartz crystal oscillator. This oscillator forms part the national frequency 
standard, hence its frequency and history known high accuracy. 

The block diagram Fig. shows the principles the control circuit. The 
fundamental frequency (100 the oscillator divided down pulses 
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Fic. diagram electrical control and timing circuit. 


per second which are used for control; the 100 signal also multiplied 
toa signal that used for defining the time scale the oscillograph 
trace. 

small part the light from the spark picked photomultiplier 
and used monitor the flash, and the resulting pulse added the 100 
pulses and Mc/s signals obtained from the crystal oscillator. The composite 
signal displayed the cathode ray tube, successive light flashes being 
automatically displaced the vertical direction. single photograph then 
shows all the light flashes for one drop. Figure shows trace due one flash. 
The small amplitude sine wave the signal, the first pulse the 
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100 kc/s pulse and the second the photomultiplier signal. Departures from 
the nominal values (0.2 second and 0.6 second) and are shown 


variations this trace form. 
All time measurements are made terms the Mc/s signal and the 


end correction applied for the departure from nominal frequency the 
crystal oscillator. 


Vacuum System 

The column evacuated through 15-cm pumping line which leads directly 
the water-cooled baffle 15-cm oil diffusion pump. This pump has 
nominal 600 liter/sec pumping speed which allows maintain pressures 
below 0.1 even the presence substantial air leaks, and ensures that 
other vapor pressures, e.g. lubricants the motors and fluid from the shock 
absorbers, are held negligible values. The column pressure measured with 
ionization gauge which has been checked situ means Knudsen 


gauge. 


Contactor 

Electrical connections for the six thermistors are made means 
contacting ring the rule and corresponding outer ring the top the 
column. This outer ring lowered into contact position means small 
motor before the full weight the rule the support plane. Two pairs 
pins and sockets position the rule the support plane before the carriage 
lowered for the next drop. 


Damping 

Impulses imparted the rule during the positioning operation cause 
pendulate. The amplitude this pendulation reduced the order 
motor-driven forked cam which brought into contact with the assembly 
and then slowly withdrawn from it. 


Release 

The rule released current pulse. The ground and lapped, hardened 
steel plane, which the spherical bearing the rule rests, maintained 
position electromagnet. When this de-energized, spring lowers the 
plane acceleration several gravities and withdraws from the path 
the falling rule. Resetting means motor-driven cam assembly. 

was not possible withdraw the support plane without giving very 
slight lateral impulse the top the rule. prevent this impulse from throw- 
ing the scales out focus compensating impulse was applied for the first 
millimeter the drop means knife-edge which pressed with constant 
tension against short pin projecting from the top the rule. 


Arrestment 

The arrestment mechanism consists essentially cage suspended four 
coil springs. The falling rule assembly received the cage and brought 
rest average deceleration gravities. The energy the stretched 
springs given two oil-filled shock absorbers which allow the cage and 
rule assembly return slowly their rest position. The initial shock 


™ 


PRESTON-THOMAS AL.: MEASUREMENT ACCELERATION 843 


impact, during which high velocity imparted the cage, mitigated 
sponge rubber cushion; the peak deceleration, determined high-speed 
camera, the order gravities. 


Rule Return 

The rule returned its initial position means motor-driven carriage. 
Because the motor generates substantial amount heat placed outside 
the column and operates through O-ring vacuum seal. 


EXPERIMENTAL PROCEDURE 


The scales having been aligned the rule, the rule was adjusted hang 
vertically, using two plumb wires and travelling microscopes short focal 
length determine the relation the upper and lower scales the plumb 
The rule was then given rough preliminary calibration, inserted into 
the apparatus, and dropped about times. This was with the intention 
seeing any gross changes the relative positions the scales occurred 
during the first few decelerations (none ever did) and adjusting the release 
mechanism until the rule was dropping nearly vertically, i.e. with rotation 
0.2 minute less. 

After this procedure the rule was removed from the apparatus, calibrated 
described above, and replaced the apparatus. Before the definitive drops 
began the rule was brought within 0.1° the temperature which was 
calibrated and the ‘degaussing’ field was switched on. 

The sequence operations for each drop was: 

Thermistor contactor lifted, focussing scale illuminated, rule mechanically 
damped, camera focussed, air conditioning switched off (to reduce vibrations), 
plate holder inserted and set running, pressure read, flash started (with shutter 
closed), amount initial delay set, seismographs started, rule dropped. The 
system was then restored the ‘ready’ position, the whole cycle taking about 
minutes. When the drops were completed the rule was recalibrated. 

The plates were developed, washed, and dried the usual manner. Each 
plate was then read Zeiss 200-mm comparator three observers. Since 
there were five reference marks, this meant that there were determinations 
position for each image. 


Measurement Photographic Plates 

The distance through which the rule falls between the first and either the 
two subsequent photographic exposures was computed from the measured 
positions the five fiducial lines superimposed each the three scales, 
from the linear calibration the scale itself, and from the previously deter- 
mined interval between the scales concerned. The procedure 
measuring the photographic plates was point each the fiducial lines 
and the adjacent graduations. 

With perfect calibrations and ideal observers would have for any one 
drop the relations: 
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where the distance fiducial line read observer from the zero 
the scale designated the superscript, and and are two constants con- 
necting scales and and scales and respectively. assess the consistency 
the readings and the accuracy the calibrations, values and were 
determined for each the sets readings made for each drop. Table 
shows the results for typical* drop (number for rule II). 


TABLE 


Comparison the results the readings the three observers 
(A, and each the five fiducial lines for typical 
drop (drop 38, rule IT) 


3.7231 


Fiducial line 
Fiducial line 
Fiducial line 


Mean all 8.3852 3.7231 


Range from mean (+0.9, 
—0.3) (+0.4, —0.5) 


Fiducial line 


Fiducial line 


additional check the consistency the readings, the value deter- 
mined from the readings each observer for all drops was computed. These 
differed from the mean value less than 0.1 milligal. 

RESULTS 

preliminary set drops was made August 1957 with the rule 
and gave value for (employing one observer only) 980.6136 

definitive series drops was made August 1958 with the same rule 
(Mk after the scales had been realigned and minor improvements made 
the apparatus. The value obtained for was 980.6139 both series 
for this rule air spark light source was employed. This rule did not have 
auxiliary focussing scale. 

December 1959 another definitive series drops was carried out with 
the rule using xenon discharge tube light source. The value for 
obtained from these drops was 980.6124 

The individual values obtained from these drops are shown Table 
and Fig. 


*Of the sets values for and obtained with this rule, half have variations greater 
and half have variations smaller than those this plate. 


se 
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TABLE III 
Dec. 1959 Aug. 1958 Aug. 1957 
Drop Drop Drop 
6135 
.6122 
Mean 980.6139 
Std. dev. +0.0006 +0.0009 +0.0022 
Range 0.0029 0.0041 0.0100 


practice, about the drops involved some mechanical human 
malfunctioning, that position, time, and satisfactory seismograph records 
were not all available. drop where all three were available has been omitted 
from the above table. the two definitive sets drops drop had 
omitted due unsatisfactory seismograph records. The mean the definitive 

CORRECTIONS 
Extension Rule under Vacuum 

This was calculated from standard data for the steel used: the correction 

amounted +0.2 milligal. 


Extrapolation Vacuum Condition 

The definitive drops were carried out pressure about 0.07 Hg. 
The appropriate correction was found experimentally. The correction 
was high (+0.5 milligal) and considerable extent indeterminate; has 
therefore been discussed separately, see Appendix 
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Maximum spread 2.9 mgal 


Drops 


A 

Drops Maximum spread 10.0 mgal 


| | | 


or NM HO OH 
© 


Fic. Histogram the three sets drops with rules and 


Non-Horizontality the Fiducial Lines 

the drops made with the rule the fiducial lines departed significantly 
from horizontality. The rule was dropping with measurable lateral shift, and 
appropriate correction, averaging 0.3 milligal, was applied each drop. 


Departure Standard Frequency from Nominal 
This was parts for the drops the rule and —9.5 parts 


for rule (—0.16 and —0.19 milligal). 


Correction for Gravitational Gradient 

The values obtained from equation (3) apply position space 
approximately below the center gravity the rule assembly the 
start the time and distance intervals (see Appendix The local vertical 
The appropriate correction reduce the results the position the optical 


axis is, therefore, +0.12 milligal. 


ERRORS 


Table gives summary those errors which could significantly influence 
measured value Additional comments are: 

Rule calibration and plate reading errors are due chiefly uncertainty 

primary calibration the interval, and lesser extent pointing 


-Dec 
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TABLE 


Probable 
maximum errors 
(1) Rule calibration and Short interval +0.62 mgal 

plate reading Long interval +0.33 

(2) Time Short interval +0.12 
Long interval +0.10 

(3) Temperature +0.08 
(4) Correction for thermal expansion +0.03 
(5) Levelling optical axis +0.08 

(6) Tilt cross hairs +0.03* 


(7) Residual gas pressure 


*This applies for the Mk II rule; for the Mk I rule e was 0.18 milligal. 


errors the plates; these latter are presumably random and are averaged 


over number readings. 
Errors time measurement, discussed earlier, may amount 


second. 

Thermal shielding, thermal contact with the rule, thermistor power input 
wwatt) and reading accuracy were such that rule temperature the 
thermistors was known better than +0.003° The measuring circuit 
allowed monitoring for electrical leakages and readings any sort were 
taken when leakage could have contributed significant error. 

The error the correction for thermal expansion arises from +1% 
uncertainty the expansion coefficient, and will only this large the 
extreme temperatures encountered (19.7° and 20.2° C). 

The maximum tilt the optical axis minutes coupled with the 
greatest horizontal motion the third scale required give the quoted 
+0.08-milligal error. 

The tilt error arises from lateral shifts the scales with respect the 
vertical cross hair the rule drops, together with uncertainty about 
minutes the amount tilt. The error much greater for the rule 
because substantial lateral shift was present throughout the drops. 

The lower limit (—0.2 milligal) for the error the residual gas correction 
has been set coincide with the theoretical value the operating pressure. 
The higher upper limit was selected result the unsatisfactory agreement 
between theory and experiment. 

number other factors that might affect the results were considered. 
These were either eliminated appropriate design the apparatus, else 
measured and corrected for. Errors, after any appropriate correction has been 
applied, arising from the following sources are therefore considered 
negligible: 

Dilation the rule, vibration the rule assembly, seismic vibrations, 
cosine corrections for tip and curvature the rule, coriolis forces, diurnal 
variations vertical gradient asymmetry the flash, flash lifetime, 
differential shrinkage the photographic emulsion, and magnetic fields. 
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The sum the probable maximum errors about milligal, not inconsistent 
with the difference between the results the two rule assemblies. However, 
0.5 milligal (the uncertainty the pressure correction) presumably the 
same sense and the same order magnitude for each rule: the light 
this, the separation the two sets results seems fairly high. Taking this 
into consideration, believe that are justified assigning possible error 
+1.5 milligal from the mean the two sets results. 


DISCUSSION 


The absolute gravity station Ottawa located the subbasement the 
National Research Council building, Sussex Drive, Ottawa, and defined 
the intersection the optical axis and the geometric axis the vacuum 
column. Its geographical co-ordinates are: latitude, 45° 26’ longitude, 
75° 41’ altitude, 57.6 meters above mean sea level. The value 
this point 7.29+0.03 milligal higher than the value the national 
reference gravity station the Dominion Observatory, Ottawa. The accepted 
best value the reference station the Potsdam system 980.6196 
(Winter al.). This gives Potsdam value for the absolute gravity station 
980.6269 Our measured value 980.6132 the mean the 
two definitive sets drops, therefore 13.7 milligal below the Potsdam 
system. The uncertainty the gravity ties the order +0.5 milligal. 
From the above considerations conclude that the correction the Potsdam 
system required bring absolute gravity units 
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APPENDIX 


Time Ratios 

require the ratio the time intervals and such that the error 

Under uniform acceleration the intervals are related 
equation (3), i.e.: 
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Where are the uncertainties measurement, the uncertainty 

—te 


S2— 


determined from experimental conditions. 


1 aK+1 vot (vot2gti)K 


practice, the contribution the second term the above equation 
small, and the optimum may determined from the first term only making 
(15) 


Gravitational Gradient 
general where isaconstantand isa vertical displacement. 
define the position the center gravity 


where the time interval between the release the rule and the first flash 
(at 0), the equation motion given by: 


double integration this equation gives 


distances and and time intervals and obtain: 


The second term within the brackets therefore gives the fractional correction 


*Assumption equal magnitudes for the time errors consistent with the conditions 
our experiments. 


> 
J 
(17) 
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applied the value determined from equation (3) referred 
the position the center gravity For the constants our 
0.6 sec) the expression the bracket has the value 


Force Proportional Velocity 

retarding force proportional the velocity the rule arises from the 
action the residual gas the vacuum chamber. For such force have 


where the mass the rule. Then have 
and the previous case can write: 


For our case, 0.3 K/m. 


APPENDIX 


Pressure Effect 

The forces exerted the rule the residual gases the column can 
divided into ‘static’ and ‘streaming’ effects. pressures below only 
‘molecular flow’ the gas need considered. The static effect then arises 
from the average momentum imparted the moving rule the randomly 
directed molecules striking it. If, probable, the coefficient slip the 
molecules zero, that say that the direction departure molecule 
from the rule independent its original direction approach, then the 
net force the rule given by: 


(23) 


where the superficial area the rule assembly, the pressure, the 
velocity the rule, the average molecular weight, and and have 
their usual significance. The effect the measured value force that 
proportional the rule velocity for our case, see Appendix given by: 


K/m 


where the mass the rule assembly. 
This can rewritten, using equation (23), as: 


Bs 
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980.614 
DROPS 
980.613 


980.612 


980.611 


980.610 


980.609 


980.608 


microns 


Fic. curve between 0.01 and 1.0 Hg. The cross the average 
value nine drops this pressure. The dotted line gives the theoretical slope for the moleeular 
flow region ‘static’ air. 


The dotted line Fig. has this pressure dependence. The solid line 
Fig. closer fit the experimental points series drops various 
pressures. These drops were made after the definitive series, and with less 
attention paid the quality the plates, attitude the rule, accuracy 
reading, etc. Their scatter, therefore, greater than that the definitive 
drops. 

The departure the results shown Fig. from the theoretical line 
suggests that systematic departure from ‘‘static conditions the given 
pressure exists. This may consist of: inaccuracies factor the pressure 
measurement; heavier than air molecules quantities sufficient increase 

take these order: 

The ionization gauge used was calibrated situ against Knudsen (thermo 
molecular) gauge and found accurate except the lowest pressures 
where the ion gauge read high presumably due the 
presence hydrocarbon molecules. 

probable that the pressure the column was due air molecules both 
during the definitive drops and during the drops taken determine the 
pressure effect; the first case inadvertent and the second case 
deliberate air leak was present, while the accuracy the ionization gauge 
both cases strongly suggests that was measuring air molecules. 

The known position and size the leaks involved and the pumping 
orifice make possible estimate the effect stream gas. This given 
equation (23) where now the velocity the gas stream relative the 
rule and the area, not necessarily the total superficial area the rule 
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assembly, which the stream acts. Taking the position the two leaks into 
account does not seem possible that streaming would affect the result 
more than 0.1 milligal pressure dyne while such effect 
estimated exist would increase rather than diminish the slope the curve 
Fig. 10. 

attempt resolve this discrepancy between theoretical and experi- 
mental slope series drops were made higher pressures. the pressure 
molecular flow and viscous drag; still higher pressures the inertial drag 
the entrained air, the pumping effect the rule, and similar effects will enter 
the picture. the intermediate region the dependence pressure was 
about four times great theory (c.f., for example, Dushman (1949)) 
indicates. 

these rather unsatisfactory circumstances appears that the proper 
correction for the pressure effect that taken from the slope the curve 
Fig. 10; this amounts, for the operating pressure Hg, 
0.5 milligal. matter regret the authors that the experiment was 
not conducted pressure the order 0.01 Hg. 
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THE WAVELENGTH DISPLACEMENTS SOME INFRARED LINES 
BETWEEN LIMB AND CENTER THE SUN. 


LuIsE HERZBERG? 


ABSTRACT 


The differences between the wavelengths the solar limb and the center 
and regions the spectrum. The values the displacements 
(in km/sec) the lines the two wavelength regions studied are found 
the same those obtained Adam for neutral metal lines 
The displacements the lines are similar magnitude and 
the same direction those Although the data are insufficient 
decide the question term dependence the solar wavelength shifts 
any relation the shifts observed laboratory light source can excluded. 
For the lines and 48900 corresponding two different transi- 
tions, the limb-center displacements differ from each other both magnitude 
and direction. The displacements the lines are 
smaller than those the lines, while those the lines are 
significantly larger and direction opposite those observed for lines 

Where possible, comparison has been made between the wavelengths observed 
limb and center the disk and the solar wavelengths predicted General 
Relativity Theory. cases the wavelengths the limb were found closer 
the predicted values than the wavelengths measured the center the disk. 
While for the lines the predicted solar wavelengths and those observed 
near the limb good agreement, the wavelengths close 
the solar limb the lines are found significantly larger 
than those predicted relativity theory. 


INTRODUCTION 


the basis General Relativity Theory all solar absorption lines are 
expected shifted the red with respect those terrestrial vacuum 
source amount times the wavelength Measurement 
this gravitational red shift constitutes one the classical tests for the validity 
theory. Although spectroscopic techniques measure wave- 
length differences the required order have been available since before the 
problem first arose, there yet definite answer the question whether 
not the solar lines show the predicted shift. The reason for this that the 
relativitistic gravitational shift, present, masked other effects whose 
origin is, far, not understood, and whose influence the total observed 
wavelength shift cannot properly assessed. 

There are two major observational facts, mutually interdependent, which 
make obvious that the wavelength shifts solar lines with respect their 
terrestrial counterparts have component (or components) other than the 
gravitational shift and the Doppler shifts due solar rotation and relative 
motions earth and sun. These observations are: (a) The differences between 
the wavelengths measured the spectrum the center the solar disk and 
the spectrum standard light source are not strictly proportional 
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the wavelength. (b) For most lines the wavelength, properly corrected 
for the effect solar rotation, varies along the radius the solar disk, resulting 

Until fairly recently, mainly the basis work St. John (1928), 
was believed that the differences between the observed and predicted solar 
wavelengths could accounted for Doppler shifts due radial motion 
the photospheric gases. However, the very accurate data solar and 
terrestrial wavelengths, now accumulated Adam and co-workers 
University Observatory, Oxford, show that the effect much too compli- 
cated interpreted such straightforward manner. The new material 
has stimulated much theoretical work the processes which cause the 
observed solar wavelength displacements. far, however, 
treatment the subject has yielded results which are satisfactory agreement 
with the observations. 

The eventual identification the physical processes causing the observed 
solar wavelength displacements depends largely knowledge the latter’s 
dependence parameters other than the wavelengths, for instance, the 
intensity, the atomic energy levels involved the transition, etc. seems 
reasonable expect that such relations will become more obvious observa- 
tions are extended longer wavelengths, since the displacements due any 
cause, general, increase with wavelength. The Oxford data, far, refer 
lines the visible region the spectrum only. Observations the photo- 
graphic infrared could therefore expected yield useful results. 
liminary study the \10700 region, made the Dominion Observatory 
(Herzberg 1957), showed systematic differences the solar dis- 
placements two groups lines, belonging, respectively, two different 
multiplets. 

The work was therefore continued with survey the dis- 
placements suitable solar lines the region the spectrum. When 
the results obtained for the lines proved special interest, the 
region was added order include the investigation the intense 
lines this wavelength. The results obtained for lines 
and are presented below. 


OBSERVATIONS. REDUCTION PLATES 


The observations were made during the summer months 1958 and 1959. 
The equipment was the same that used the previous investigation: 
horizontal telescope producing solar image 230 diameter; image 
rotator which makes possible place section the solar image any 
desired direction over the slit the spectrograph; 20-ft spectrograph with 
15,000 lines/inch plane grating. 

spectrum from point near the solar limb and one from the center 
the disk were always taken immediately after each other without disturbing 
the plate. this way, the use suitable diaphragms over the spectrograph 
slit, were obtained which the spectrum from point near the 
limb bracketed the spectrum from the center (see Fig. 3). Pairs such 
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strips from opposite limb points have been evaluated units. They are re- 
Only those plates were accepted which both strips 
had been obtained within small time interval and with the same setting 
the image rotator; this way errors due inaccuracy the determination 
the direction and variations the size the solar image are 
automatically excluded. 

All exposures were made with the spectrograph slit perpendicular the 
solar diameter along which the observations were made. order obtain 
the maximum difference between center and limb spectrum which advanta- 
geous for accurate measurement (see below) the plates were general taken 
along the diameter. The guiding accuracy estimated have been 
+1mm limb distance and +2° latitude angle. 

The spectra were taken the third order the grating (dispersion 0.7 
A/mm) with slit width microns. The exposure times were the order 
seconds and minutes respectively, the 8500 and regions 
(Kodak spectroscopic plates IN, not hypersensitized). 

The plates were measured conventional measuring machine with 
field view large enough accommodate the whole width strip 
spectra. The cross hair was set successively the center section (limb 
spectrum) and the two outer sections (spectrum from center disk). 
difficulty was experienced focusing attention one spectrum time. 
therefore believed that appreciable errors have been introduced the 
measurements due the simultaneous presence the field view two 
different spectra. The wavelengths given Babcock and Moore (1947) for 
lines the center the disk and for terrestrial lines were used obtain 
the dispersion. This sufficiently accurate determine the small wavelength 
differences between center and limb spectra. 


RESULTS 
From every plate pair values 
can obtained for each solar line, corresponding the total wavelength 
difference between the two opposite limb points, respectively, and the center 
the disk. After the usual corrections for the motion the earth have been 


made, simple addition subtraction yields values for the dis- 
placement and the Doppler shift due solar rotation respectively: 


(limb-center displacement; limb shift) 
(Doppler shift due solar rotation). 


The results obtained are presented Tables and II. The values tabulated 
are averages results derived from measurements number individual 
plates, indicated. The slight fluctuation the number plates due 
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TABLE 


displacements and rotational Doppler shifts solar lines 
the region 


r/R 0.982 r/R 0.956 
(A) (A) plates (A) (A) 

8471.744 +0.009 0.044 +0.011 0.050 
8501 +0.009 0.046 +0.003 0.046 
8514.082 +0.011 0.055 +0.012 0.054 
8515.122 +0.012 0.051 +0.010 0.050 
8526 +0.012 0.059 +0.008 0.054 
8542.144 —0.018 0.054 —0.020 0.045 
8582.271 +0.012 0,054 +0.011 0.052 
8597 .059 +0.011 0.054 +0.007 0.048 
8598 836 +0.010 0.048 +0.009 0.051 
8611.812 +0.013 0.052 +0.012 0.052 


Note: N, hazy. 
*Center of solar disk, Babcock and Moore (1947). 
tReduced tor/R = 1. 


TABLE 


displacements and rotational Doppler shifts solar lines 
the 48900 region 


r/R 0.982 r/R 0.956 


platest (A) platest (A) (A) 
Fe 8920 036 3 6 +0014 0 040 3 +0009 0.051 
Ca II 8927 6 +0. 007 0.054 3 +0. 006 0.060 
Fe $943 05% 2 6 +0.012 0. 058 3 +0. 009 0 052 
Fe 8945. 198 5 6 +0. 012 0. 057 3 +0010 0 062 


Note: N, hazy. 

*Center of solar disk, Babcock and Moore (1947). 

TReduced to r/K = 1. 

Tesults from one plate, taken at latitude 440° (see text) not included in Adsot 
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the fact that occasionally, otherwise satisfactory spectrum taken near 
the limb, the contrast for certain lines was too low make measurement 
possible. 

The tables contain results obtained two different disk positions, 
r/R 0.982 and r/R 0.956 (r, distance from center; radius solar 
disk), corresponding limb distances the solar image and 
respectively. The values 5-mm limb distance had been obtained early 
stage the work when observations were made series different limb 
distances the hope determining for the different lines the wavelength 
variation along the solar radius. After was recognized that the method was 
not suitable for this purpose, further work was restricted 2-mm limb 
distance. However, the values obtained for 5-mm limb distance are recorded 
here with those for 2-mm limb distance because the degree consistency 
between the two sets data provides measure for the reliability the 
results general. 

far know (Freundlich, Brunn, and Briick 1930; Herzberg 1957), 
the limb-center displacement solar lines independent heliographic 
latitude. The observations were therefore not diversified position angle, 
but, for reasons explained above, were made along the east-west diameter. 
Results from two plates, taken position angle +30° owing over- 
sight, have been given full weight when forming the average values for the 
displacement but, because the uncertainty the 
law solar rotation they have been rejected the computation the 
rotational Doppler shift 

limits error are indicated the tables because there are not enough 
single measurements compute the conventional m.s.e. values. The accuracy 
the final results can estimated from the scatter the values obtained 
from individual plates, plotted Fig. for some the lines. This scatter 
due combination measuring error and guiding error. Since each 
the values for the displacement the rotational Doppler shift 
derived from measurements four different spectra, twice the error 
the single measurement has allowed for. 

Some the variations the values for the displacement 
different lines are clearly outside the limits accuracy. These will dis- 
cussed below. The fluctuation the values for the rotational Doppler 
however, does not seem significant, and mean values, taken over the 
different lines the same wavelength region, have been used determine 
the rotational velocity the sun. 

Values for the linear velocity the equator due rotation, have been 
computed both from the average rotational Doppler shifts all lines and 


from those the lines intensity The results obtained are 


follows: 


otf 
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Fic. Wavelength shifts due solar rotation (triangles) and dis- 
placement (circles and diamonds) for different solar lines obtained from measurement 
individual plates. (Upright triangles and circles refer measurements the 
diameter, inverted triangles and diamonds measurements diameter position angle 


+30°.) 


These values are within the range other spectroscopically determined 
values the solar rotational velocity, for instance 


1.956 km/sec region (Hart 1954)) 


and 
1.85 km/sec region (Herzberg 1957)). 


DISCUSSION 


(a) Lines Neutral Iron, (and Nickel, 

The values for the displacements the (and lines 
48500 and show, within their respective wavelength regions, 
noticeable variations with intensity, excitation potential, other properties 
the lines concerned. The average values for the two groups lines can 
therefore considered fairly representative for their wavelength regions, 
and can compared with the corresponding value for neutral lines (mostly 
Fe) 46100 (Adam 1948) order study the wavelength dependence 
the effect. 
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Excluding infrared lines intensity because their lower accuracy, 
the average values for the displacements the three regions 
the spectrum are: 


These three values, within the accuracy guiding, refer the same disk 
position (r/R 0.982 for and r/R 0.984 for A). 

are 


The slight increase these values with wavelength only the order 
magnitude the estimated error; 0.02 km/sec corresponds 0.0004 
6100 and 0.0006 9000 Within the limits the present observa- 
tional data therefore, the solar displacement neutral metal 
lines proportional the wavelength. 

This result makes not unreasonable expect similar linear dependence 
the wavelength exist also for the line shifts and, there- 
fore, for the differences between the observed solar wavelengths and those 
predicted General Relativity Theory 

For the lines the region, the relation between predicted and 
observed solar wavelength shifts has been investigated Adam (1948), who 
has measured directly the wavelengths different points across the disk and 
compared them with the wavelengths standard laboratory light source 
(vacuum arc). was found that the average wavelength near the limb 
(r/R 0.984) corresponded closely that predicted General Relativity 
Theory, while the average wavelength the center the disk was approxi- 
mately halfway between the wavelengths measured the limb and the 
laboratory. 

For the infrared solar lines, wavelengths the center the disk have 
been measured Babcock and Moore (1947). From these the wavelengths 
the solar limb can obtained adding the values for the 
displacements (Tables and Unfortunately, corresponding wavelengths 
laboratory light source vacuo are not available far. However, for 
approximately half the infrared lines considered, Meggers (1935) gives the 
wavelengths measured the spectrum arc air atmospheric pressure, 
together with corrections applied convert these wavelengths into 
those vacuum light source. one accepts the wavelengths computed 


laboratory wavelengths, one obtains: 


this way 
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for \8500 


(mean 
and 
for 
(mean 


These values, together with the corresponding data for the region 


are shown, the left, the diagram Fig. The wavelength shifts are given 
km/sec (see above). these units 


0.63 


The vertical solid lines connecting the points for and 
represent the values the limb-center displacements The dotted 


+4 
2 
= 
(predicted) 
» (lab.) 
= o o 


Fic. Average shifts with respect standard laboratory light source groups solar 
lines (see text) the center the disk and point near the limb. (For the lines 
8900 because lack good laboratory data, the shift between the wavelength the 
the solar disk and the laboratory has been arbitrarily assumed equal 
zero. 
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line connecting the point for with the abscissa has been 
drawn stress the fact that, this wavelength region, the difference 
between solar and terrestrial (vacuum arc) wavelength known the same 
accuracy the value for the displacement. 

glance the diagram shows that, for the lines the \8500 
region, the solar wavelength shifts (in km/sec) are indeed practically identical 
with those observed \6100 For the \8900 lines the shifts between 
solar and laboratory wavelengths are larger. However, for this group lines 
the laboratory wavelengths are uncertain that the values obtained for 
the solar terrestrial displacements (different from those for the 
displacements) have too wide margin error relevant for the present 
discussion. The evidence obtained therefore favors the conclusion that for 
lines iron and similar neutral metal atoms, the infrared well 
the region, the solar wavelengths observed point near the limb 
are close those predicted General Relativity Theory. 


(b) Lines Neutral Silicon, 

Tables and show the values obtained for the solar 
displacements the lines studied, together with the laboratory wave- 
lengths and other relevant data, arranged according multiplets. 

The laboratory wavelengths the lines the literature (Kiess 1938) 
have not been obtained under very good vacuum conditions and must 
assumed affected disturbing factors the light source. They are 
therefore referred Moore (1956) has found that the wavelengths 


TABLE III 
Wavelength displacements lines 
Part 

Excitation 

potential 


Note: N, n, hazy; w, widened, not shar 

*Center of solar disk, Babcock and ious (1947). 
tKiess (1938). 

tMultiplet No., assigned by Moore (1945). 
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TABLE 
Wavelength displacements lines 
Part 
Excitation 

(45)t 


Note: N, n, hazy; w, widened, not sharp; 1, unsymmetrically shaded to longer wavelengths; P, predicted. 
*Center of solar disk, Babcock and Moore (1947). 

tKiess (1938). 

tMultiplet No., assigned by Moore (1945). 


measured the center the solar disk provide more consistent system 
atomic energy levels than the present laboratory data and that therefore 
the solar wavelengths can taken the better approximation the 
wavelengths perfect vacuum laboratory light source. (The relativistic 
shift solar lines within the accuracy the laboratory wavelengths.) The 
values (A° can therefore considered measure for the magnitudes 
the wavelength shifts produced the laboratory light source essentially 
pressure effects. 

pressure effects had noticeable influence the wavelengths solar 
(photospheric) lines one would expect the effect greater the center 
the disk than the limb, because the decrease pressure with increasing 
height the photosphere. The values for the limb-center displace- 
ments would then analogous the differences and should, 
pressure were the main governing factor, show similar term dependence. 

Results obtained earlier investigation (Herzberg 1957) lines 
belonging multiplets No. and No. (Moore’s (1945) notation) seemed 
give indication that such effect might exist. The results for these 
multiplets are included Table can seen, the values for 
well those for are consistently larger for multiplet No. 
than for multiplet No. order see whether not this agreement 
more than coincidence, special attention was given lines the present 


investigation. 

The best test object amongst the lines studied provided multiplet 
No. 54, although unfortunately only one its lines intense enough 
measured disk position r/R 0.982. The lines this multiplet and those 
multiplet No. have the same lower energy levels, but their values for 
are opposite directions. the values for the dis- 


| 
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placement were really related the differences the lines multi- 
plet No. should show very small, perhaps even negative values for 
The results observation (see Table III) show that this not the case. The 
limb-center displacements measured for the lines multiplet No. are 
the same order magnitude and the same direction those the 
lines belonging multiplet No. 53. 

Another case interest the present context the unclassified line 
\8648.472 (Table IV). This line has exceptionally large shift 
the same direction lines multiplet No. (negative); addition 
asymmetrically broadened towards longer wavelengths. Even this line 
has been found have displacement average magnitude 
and the usual direction (positive). 

Although the present data are not sufficient exclude altogether the 
possibility that the displacements the solar lines are term 
dependent, they definitely show that there direct relation between these 
displacements and the (inferred) wavelength shifts observed the laboratory 
light source. 


(c) Lines Singly Ionized Calcium, 
The values obtained for the displacements the lines 
are presented Table together with other relevant data. Only for four 


TABLE 
Wavelength displacements lines 


Excitation 

(A) potential 


solar disk, Babcock and Moore (1947). 
(1937). 
Shenstone (1946). 


the five lines the plates has been possible derive value for 
the displacement; the fifth line, could not measured 
the spectrograms obtained points near the solar limb because blending 
with the line (see Fig. 4). 

electronic transitions. The \8500 lines have low excitation potential (1.7 ev) 
and are exceptionally intense and broad; the lines have high excita- 
tion potential (7.0 ev), and their intensity and general appearance that 
average line the solar spectrum. 

The differences between these two groups lines are reflected remark- 
able manner their limb-center displacements: the values obtained for the 
low excitation lines are more than twice large those for the high excita- 


864 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


tion lines, and the opposite direction. The good agreement between the 
results obtained from individual plates, plotted the diagram, Fig. 
shows that this effect definitely outside the limits error. Inspection 
the diagram and Tables and shows further that the dis- 
placement the high excitation lines the direction parallel that 
observed for the lines and Sil, while the displacement 
the low excitation lines ‘‘reversed’’. 

Although the very broad lines are not easy measure, 
asymmetry should have been introduced the measuring process, because 
all spectra have been measured both the direction increasing and 
decreasing wavelengths. Actually the behavior the low excitation 
lines sufficiently different from that the average solar lines that can 
noticed inspection the plates (see Figs. and 4). 

For the low excitation lines standard laboratory wavelengths 
have been provided Wagman (1937). From these, together with Babcock 
and Moore’s (1947) wavelengths for the center the solar disk and the 
values for the displacements given Table one obtains the 
following average wavelength shifts for the lines \8498.062 and 


lab 


The corresponding values for the high excitation lines cannot 
determined, due lack reliable laboratory data. 

the diagram Fig. the results obtained for are shown side 
side with those for (The value for the high excitation lines 
has been arbitrarily set equal zero.) 

seen that the magnitude the (average) displacement, 
given the length the solid vertical lines, relatively large for the low 
excitation, relatively small for the high excitation lines II. appears 
significant that the low excitation lines, addition the reversal 
the sign the limb shift, show shifts between the wavelengths the center 
the solar disk and the laboratory which are considerably greater than 
those predicted relativity theory. Therefore, for these lines well for 
the lines, the solar wavelengths measured the limb are 
closer the predicted wavelengths than those measured the center the 
disk. However the case the lines, different from that the 
lines, even close the limb 0.982) the solar wavelengths are con- 
siderably larger than the predicted wavelengths. 

The only case large limb shift that has been observed 
previously seems that the central absorption the doubly 
self-reversed resonance line \3933.684 (St. John 1910). This line 
and the lines have the upper energy level common. addition, 
all three lines are exceptionally intense the solar spectrum, and must 
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therefore formed high the solar atmosphere. Further observational data 
will needed decide which the foregoing facts are relevant explain 
the similar and exceptional solar wavelength shifts these lines. 
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distances from the east and west limb and the center the sun. The two lines, seen well 
separated the spectrum the center the disk, merge close approach the limb. 
This effect very likely due the fact that the limb shifts the two lines are opposite 
directions. The unmarked spectra are from the center the disk, those marked and 
(d) r/R 


NOTES 


GAMMA RADIATION THE DECAY 


KJELBERG, TANIGUCHI, AND YAFFE 


Gamma rays have been reported recently the decay 5.3-hour 
Alexander al. (1958) found gamma ray 0.31 Mev low abundance. 
Workers the Ames Laboratory (1959) observed 0.27- and 0.30-Mev 
coincidence. 

Coulomb excitation studies Wall (1954), Mark al. (1955), Temmer 
and Heydenburg (1955), Paulissen and Mark (1955), and McGowan and 
Stelson (1958) have established excited states 0.300, 0.582, and 
0.675 Mev, together with transitions 0.282 and 0.375 Mev. Alexander 
al. (1958) observed y-rays reasonable agreement with these the decay 
1.2-minute together with 0.14- and 0.56-Mev coincidence. 

Recent work this laboratory has indicated that the decay 
more complex than indicated previous investigators. 

was produced fission product bombarding thorium metal 
the McGill synchrocyclotron with 45-Mev protons for period minutes. 
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Silver was separated radiochemically method essentially due Glen- 
denin (1951). The use short bombardment times and rapid chemical 
separations, followed decay period hours, gave samples 
essentially free 3.2-hour (which semishielded its parent 
21-hour The gamma spectrum was obtained with 3in. 
crystal feeding into 100-channel pulse height analyzer. 

Fig. shown typical spectrum obtained this manner. The solid 
line indicates the raw data obtained after small correction for background 
had been made. The Compton contribution and amount backscattered 
y-radiation were determined using well-known y-standards under 
identical with those the original measurement. The photopeaks which 
resulted are indicated the dotted lines the figure. All y-peaks decayed 
with 5.3-hour half-life. small contribution 10%) the neighborhood 
the 0.30-Mev peak, due 7.5-day was measured after the decay 
and corrections were made indicated. 

The results, corrected for crystal efficiencies according Heath (1957), are 
given Table 


TABLE 
Gamma rays arising from the decay 
Energy (Mev) Relative abundance 

0.12+0.03 
100 

0.58+0.02 
0.88+0.02 4+1 
0.98+0.02 5+1 
1.18+0.02 4+1 


The y-energies observed are essentially agreement with those observed 
for transitions between levels except for the three highest energies 
which previously have not been found. There were also indications 
0.27-Mev y-ray low abundance. This has been included the 0.30-Mev 
peak for intensity calculations. The 0.38-Mev y-ray reported Alexander 
al. (1958) the decay and McGowan and Stelson (1958) 
was not observed. The latter found that 16% the transitions proceed from 
the 0.675-Mev level through the 0.300-Mev level, giving 0.375-Mev tran- 
sition. This intensity was just outside the limit detection the present 
experiment. 

gamma rays with energies above 1.18 Mev were observed. These would 
have been detected, present, about the 0.300-Mev y-ray. 

tentative decay scheme for using all presently available data, 
given Fig. This preliminary work has not established whether the 
0.88 Mev y-ray arises from level 0.88 Mev whether represents 
transition between the 1.18- and 0.30-Mev levels. Since the necessary spectro- 
meter and coincidence equipment are not available our laboratory the 
moment proceed the further elucidation this scheme, these data are 
now being published. 
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NOTE THE HYPERFINE STRUCTURE THE STATE 
BORON AND 11* 


attempt detect “hyperfine structure anomaly” for boron and 
the hyperfine structure the ground states this pair iso- 
topes has been measured the atomic beam magnetic resonance method. 
Preliminary measurements (Title 1957 quoted Brouwer and Petch 1958) 

*Issued N.R.C. No, 
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suggested that appreciable anomaly existed. true, this would contradict 
theoretical expectations. The present, more precise measurements show, how- 
ever, that the anomaly zero within parts 10°. 

The apparatus and techniques used were essentially the same those 
described Wessel and Lew (1953) but with certain improvements instru- 
mentation facilitate more precise measurements. The oven was, the 
case silver and gold, graphite tube directly heated high current 
temperature about 2000° 2100° However, was found impractical 
and unnecessary use thoria crucible contain the boron. 

For the study boron lines, the naturally occurring isotopic mixture was 
used. For the study the boron lines, the separated isotope was used. 

the numerous Zeeman lines that were observable, three lines each 
isotope were chosen for careful measurement. These were, the (F, 
designation: 


(2,1) 


9 


seen from these relations that the first line each group field-inde- 
pendent first order and that the other two lines have equal and opposite 
field dependences the same degree approximation. Thus the average 
the latter two gives directly except for small second-order correction. 
fields the order gauss, the second-order correction the order 
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.0005 Mc/sec for and .001 Mc/sec for Hence not necessary 
know accurately obtain accurate value from the observations. 

The majority the measurements the h.f.s. intervals were made 
fields less than gauss. The widths the lines observed were from 0.040 
Mc/sec 0.050 Mc/sec. The signal-to-noise ratio ranged from 

For the determination only was studied. series measurements 
were made field nominally gauss. The low frequency line was 
measured alternately with the low frequency line The line widths 
were about 0.070 Mc/sec. The magnetic field was calculated from the 
frequency using the constants given Kusch and Hughes (1959), 
viz. Av(Na) 1771.6262 Mc/sec, 2.002309. The nuclear magnetic 
moment used was the diamagnetically corrected value tabulated Walchli 
(1953) +2.217529 nuclear magnetons. These measurements could 
summarized the pair corresponding frequencies, 35.694+0.005 
48.145+0.007 gauss. 

The results the observations are follows: 


The disagrees with the earlier preliminary measurement Wessel 
(1953), who obtained value Mc/sec. The observed value 
agrees with the theoretically expected one 0.66589 for pure 
Saunders coupling with the anomalous magnetic moment the electron taken 
into account. 

The ratio the magnetic dipole interaction constants 


This compared with the ratio the nuclear g-factors found 
Brouwer and Petch (1958) 


Thus the ‘‘h.f.s. defined 


expected from the reduced mass correction but this beyond the 
accuracy the measurements. 


and 1958. Can. Phys. 36, 632. 

and 1959. Handbuch der Physik, Vol. (Springer- 
Verlag, Berlin). 

1957. Quoted Brouwer and Petch 1958. 
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INTERNAL PAIR ANGULAR CORRELATION SPECTROMETER FOR 
DETERMINATION Y-RAY MULTIPOLARITIES* 


BARTHOLOMEW 


The angular correlation between the positive and negative electrons pro- 
duced y-ray internal pair conversion index the multipolarity 
the y-ray transition. the purpose this note point out the possible 
advantages possessed certain design 180° flat field pair spectrometer 
for measuring this correlation under conditions good energy resolution (of 
the order 1%). 

Direct measurements angular correlations internal pair electrons, with 
electron counters but without energy analysis, have been employed Devons 
and Goldring (1954) for determining multipolarities isolated y-rays. Similar 
measurements have been made Gorodetsky (1956), who achieved 
some energy resolution using scintillation detectors and summing the voltage 
pulses for coincident electrons and positrons. Internal pair angular correlations 
have also been measured with the aid cloud chamber Harries (1954). 

The ‘type-3’ slits described Bartholomew, Campion, and Robinson 
(1960) effectively define conical acceptance solid angles for the electrons and 
positrons transmitted the detectors from the source which, Fig. 
assumed for simplicity localized point midway between the de- 
tectors. The acceptance solid angle for each detector approximately 
steradians where the ratio the aperture gap width, the plane sym- 
metry (the plane the page, Fig. 1), the distance from the inner jaw the 
slit the source. The line profile obtained with these apertures determined 
largely the geometry the apertures and source and only small ex- 
tent the electron pair distribution. For point source the theoretical base 
width the line equal units the y-ray energy, and the resolution 
(full width line half maximum intensity) approximately equal 
6/2. Clearly rotating the detectors about the source, possible record 
the coincidence counting rate between the detectors function the angle 
between the principal axes the acceptance cones. 

order estimate the degree discrimination, for different multipoles, 
the electron angular correlations observed with arrangement such 


*Issued A.E.C.L. No. 1008. 
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DETECTOR DETECTOR 


Fic. Schematic geometry the spectrometer. The ‘type-3’ detector apertures are 
shown section the plane symmetry. They define conical acceptance solid angles for 
particles transmitted from the source the detectors. The magnetic field the 
plane the paper. rotating the detectors about the source, the variation the coinci- 
dence counting rate function the angle between the acceptance cones may obtained. 


shown Fig. use the theoretical expressions given Rose (1949). 
compute, for each multipole the intensity, cos W,, pairs 
per quantum with the angle between the directions emission the particles 
equal and equal positron and electron total energies, W_. The 
correlation observed detecting pair coincidences counter angles 
(For simplicity the assumption made that the emission 
angles all pairs transmitted the counter angle are equal exact 
calculation the correlation would include integration over all values 
allowed the acceptance solid angles.) The ratios and 
for various are plotted function y-ray energy Fig. The ratios, 
which compare the angular correlations for M1, E2, and 
radiation with that for radiation are nearly independent energy; values 
from Figs. and that experimental determination for given y-ray 
should provide practical method discriminating between the various 
multipoles provided the maximum counter angle not too small. practice, 
the finite acceptance solid angles presented the detector apertures, the 
finite area the source, and the scattering the pairs the source would all 
tend blur the angular correlation that the ratios would less than 
shown Fig. and the discrimination between multipoles would probably 
not pronounced shown Fig. However, not expected that 
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Y-RAY ENERGY Mev 


Fic. The ratio the number pairs with emission angle the number with 
function y-ray energy for M1, and radiation and 90°, 60°, and 
45°. The multipole order and maximum emission angle are given the curves. 


these effects will large entirely vitiate the method. remains 
shown that such spectrometer would have acceptable counting efficiency 
near resolution. 

order compute the absolute efficiency the spectrometer for 
require the angular correlation function for the internal pair particles with 
respect axes fixed the spectrometer. Using equations and the paper 
Rose (1949), transformation angle variables followed integration 
over the azimuthal angles leads the following expressions for the number 
pairs, with W_, per quantum energy the interval cos 
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Fic. The discrimination factor, for various multipoles function 
emission angle, 


cos where and are the positron and electron polar angles 
measured with respect the Z-axis Fig. For radiation: 


and for radiation: 


2 3 2 2 $ 


4 


2 


and 


1/137. 
Following Bartholomew a/. (1960) the line profile (pair counts per y-quantum 
vs. strength magnetic field) produced the type-3 apertures given 
value may obtained from the expressions: 


1 1 (k°—4) 


where the momentum parameter, measures the departure 
the momentum either particle from the momentum equal energy divi- 
sion the magnetic field parameter, measures the difference 
between the magnetic field, and the maximum field which coincidences 
can occur, Ho; and where the lesser and computing the line 
profiles for various values the peak counting efficiency curves shown 
Fig. are obtained. Mev and resolution these efficiencies are com- 
parable those obtained experimentally with the internal pair spectrometer 
lens type (Alburger 

The line profiles for cannot easily computed. However, shown 
Bartholomew al. (1960), the line profile depends strongly the aper- 
ture geometry and only weakly the electron angular distribution. Therefore, 
there seems reason expect that the line shape for would 
differ radically from that for Moreover, the position the high-field 
end point the line, Ho, should suffer appreciable shift varied. 
Two causes such shift are (a) the recoil energy imparted the nucleus 
varies with and for source finite area and the sum the 
distances from the point origin pair the detector apertures depends 
the (Fig. the point origin and, general, greater 
than that for Since the momentum imparted the nucleus the 
order moc the first effect negligible even for light nuclei. The second effect 
may shown produce change the end point less than 0.1% (and 
may therefore also ignored resolution) the source width less 
than about one tenth the distance separating the detectors. shown 
Bartholomew al. (1960), restriction the source small dimensions 
also necessary obtain good resolution with apertures this type. 


*It assumed that the distribution function which holds for close approxi- 
mation the true distribution over the range momentum division implied the limits 
integration for 

efficiencies are expressed counts per transition for the 6.06 Mev nuclear 
pair line. compare with the present efficiencies his efficiencies must multiplied the 
internal pair conversion coefficient (Rose 1949). 


4 

4 

7 
NOTES 875 


876 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


An 
MI-0.5% @ 
10-9 


ENERGY Mev 


values the resolution. Curves for both and radiation and resolution and 
three points for radiation higher resolutions are shown. 


The present note intended only demonstrate the feasibility this 
method measuring y-ray multipolarities. mentioned above, account 
has been taken the effects finite solid angle, finite source area, and 
scattering electrons and positrons the source, and attempt has been 
made calculate the exact line profile and peak intensity for the 
uncertainties the expected values the correlations (Figs. and resulting 
from neglect these effects must added the uncertainty arising from 
approximations the theory internal pair production (Rose 1949). Further 
investigation the practicability the technique and empirical calibra- 
tion the correlations, using transitions known multipolarity, can 
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carried out experimentally. Modifications existing pair spectrometer 
(Kinsey and Bartholomew 1953) are presently being made for this purpose. 
These modifications include the installation mechanism for changing 
the angle between the detectors and the introduction shielding between 
the source and detectors prevent detection y-rays emitted the source. 
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ALBURGER, 1956. Rev. Sci. Instr. 27, 991. 

BARTHOLOMEW, A., CAMPION, J., and 1960. Phys. 38, 194. 
ARMBRUSTER, R., and CHEVALLIER, 1956. radium, 17, 548. 
1954. Proc. Phys. Soc. (London), 67, 153. 

1949. Phys. Rev. 76, 678. 


RECEIVED 1960. 

Puysics 

ATOMIC ENERGY CANADA LIMITED, 
RIVER, ONTARIO. 


decays principally orbital-electron capture followed the 
emission single y-ray with energy about 0.67 Mev (Wapstra al. 
1953; Robinson and Fink 1955; Bhatki 1956). number other y-rays 
have been observed, but their intensities are less than that the 
0.67-Mev radiation (Bhatki a/. 1956). The half-life has been measured 
6.2+0.2 days (Robinson and Fink 1955) and 7.1 days (Camac 1944). 

Because the energy the intense y-ray lies very close that the 
accurately known 0.6616-Mev y-ray lends itself precise deter- 
mination the relatively simple techniques scintillation spectroscopy. 
the same time, the predominance this line the spectrum provides 
means measuring the half-life without interference from possible 
impurities. the work reported here these two measurements have been 
combined. 

Ten grams high-purity (99.99%) cesium chloride were irradiated for 
hours with 30-Mev rays from the Queen’s University synchrotron. this 
bombarding energy the principal reaction cesium 
Thirty-three-minute Cl* also produced and must given time decay 
away before measurements are started. 

The source was placed jig which located reproducible distance 
in.) from in. Nal crystal mounted DuMont 3-in. photo- 
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multiplier tube. The resolution the counter 0.662 Mev was 8.6%. The 
pulses from the counter were amplified and fed into Baird-Atomic 20-channel 
analyzer. The gain the spectrometer was adjusted that the photopeak 
the 0.67-Mev y-ray was distributed over about channels, and the counting 
rates were determined better than the peak. source com- 
parable activity was counted alternately with the sequence measure- 
ments consisted two groups three, arranged minimize the effects 
drift: 137-132-137. Between the two groups the 0.84-Mev 
photopeak from was scanned with the spectrometer give second 
point for the energy calibration the bias control. The above sequence was 
repeated approximately 2-day intervals over period almost three half- 
lives. 

For the half-life determination the relative activity the source the 
time each sequence measurements was obtained comparing the areas 
under the photopeaks with those under the photopeaks. (Since the 
half-life about years, provides effectively constant reference 
peak.) The ratio the mean area for the three runs the mean area for 
the three runs gave point the decay curve corresponding the mean 
time the observations. The decay curve (Fig. shows departure from 
exponential variation with time. least-squares fit the data gave 
half-life 6.48+0.03 days. 

The energy the y-ray was derived from the same measurements. For 
each sequence measurements the difference between the average pulse 
heights for the and photopeaks was determined and compared with 
the difference between the Mn* and peaks. (Since the energy 
difference about times the difference, not necessary 
know the energy the y-ray accurately the energy.) The two 
most precise determinations the y-ray energy are 0.66160+0.00014 Mev 
(Muller al. 1952) and 0.66165+0.00015 Mev (Lindstrom 1953). 
mean value 0.66162+0.00010 Mev was assumed for the purposes the 
present experiment. The best published values for the energy the y-ray 
are 0.840+0.007 (Maeder al. 1954) and 0.835+0.015 Mev (Deutsch and 
Elliott 1944). The same transition follows the decay and the energy 
this case has been measured Mev (Schardt and Dropesky 
1956). The mean value used here Mev. the basis the 
and energies, the energy the photopeak was found ex- 
ceed that the photopeak 0.00623+0.00019 Mev. (This the 
average nine determinations, each based three measurements each 
peak.) This energy difference will sensitive the presence weak photo- 
peaks lying under the main peak, which could shift the observed position 
the latter. coincidence investigation this energy region revealed the 
presence line 0.62 Mev with intensity between 0.5 and that 
the main line. The superposition these two peaks leads downward 
energy shift 0.0001 Mev the apparent position the main peak. 
The energy difference between the and y-rays therefore 
0.0003 Mev. The resulting energy for the y-ray 0.6679+0.0004 Mev. 
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DAYS 


Fic. curve for 


Previous measurements have yielded Mev (Wapstra 1953), 
0.673+0.006 Mev (Finston and Bernstein 1954), and Mev 
(Robinson and Fink 1955). 

practical matter connected with may worth noting. view 
the probable use large CsI crystals for the measurement high-energy 
experimenters may find that they have produced this radioisotope 
within the crystal volume. The energy the observed line this case will 
include that the rays following the electron capture process. the case 
K-capture this will add another 0.0346 Mev, and the observed photopeak 
will occur 0.702 Mev. 


The authors wish thank the National Research Council and the Atomic 
Energy Control Board for financial support. 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month issue. proof will sent the authors. 


Frequency Measurement Standard Frequency 


Measurements are made Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
reference standard (with assumed frequency 192 631 770 c.p.s.). Frequency deviations 
from nominal are quoted parts per negative sign indicates that the frequency 


below nominal. 


Date, 
March 


& 


31 
Midmonthly 
mean 


Midmonthly 
mean of WWV 


GBR, 16 kc/s 
MSF, 6-hour 24-hour WWVBt 
60 kc/s average* average 60 kc/s 
— 184 —174 —174 -78 
—181 —174 —176 —78 
— 187 —179 —30 
N.M. —176 —178 ~~84 
—179 —176 -177 N.M. 
—179 -—178 -179 N.M. 
‘“—170 —178 -178 
—190 —180 N.M. —90 
— 184 N.M. N.M. -79 
N.M. — 183 N.M. N.M. 
—188 N.M. N.M. —82 
—188 —189 — 187 N.M. 
—190 —189 — 188 N.M. 
— 186 —192 — 187 —80 
—185 -—191 — 187 —83 
N.M. —188 — 188 —80 
—193 —191 —191 —80 
— 166 —-161 — 163 —81 
—159 — 162 — 160 N.M. 
—158 —159 —158 N.M. 
—159 —161 —157 —82 
—152 —157 —155 —81 
— 148 —154 — 153 —82 
—145 —154 —152 —80 
—154 -—151 -151 — B84 
—152 —152 —149 N.M. 
— 129 —149 —145 N.M. 
—135 — 148 —141 -77 
N.M. —145 —143 —82 
— 143 — 139 —146 — 84 
N.M. —161 —156 —78 
— 168 — 169 — 168 —-79 


Note: N.M. No measurement. 
*Time of observations: 00.00 to 06.00 U.T. 
+Call signals for station KK2XEI have been changed to WWVB. 
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Measurements the Last Few Periods Sputnik III 
Radio Direction Finder* 


direction finder Ottawa was used track Sputnik III during the last weeks its 
life. The frequency the transmission from the satellite was 20.0045 megacycles per second. 
The tracking data were primarily used for another experiment and also provide immediate 
information the position the satellite aid possible visual sightings. Data this 
nature, however, could used with the more accurate but fewer visual observations give 
fuller picture the falling-in process. 

Observations were made automatically facsimile recorder (McLeish 1958) 
cessed give the time bearing (Bo), and distance nearest approach (Wolfe 
1958). The time crossing the 45th parallel was calculated plotting the tracks 
map. The mean period was derived from the values successive 16th crossings. 

The initial values and were derived assuming the eccentricity and the 
period 90.0 minutes for all crossings. The distance nearest approach was therefore 
recalculated using the approximation that the angular velocity the satellite its orbit 
where the angular distance (in the orbital plane) between the 
satellite and the position perigee. 

The values and were derived from data published the Smithsonian Observatory. 
The first estimates were used get which differed from only +1.5% 

New values 74s were derived the same manner before and are listed Table 
along with the time nearest approach Ottawa 75°.34’W.). All the observed 
crossings were from the north the south. 


TABLE 
Date (U.T.) Tas (U.T.) 

March 1810.0 1808.5 

1810.4 

1638.7 

1506.9 

1637 1636.8 

1502.5 

1457 1458.2 

1450.9 1451.3 

1442.5 1442.6 

1431.4 1431.5 

1400.3 1400.1 

1314.7 1314.7 

1244.9 1245.1 


observation was made April because equipment maintenance. 
The calculated mean periods and the corresponding times are given Table II. 


TABLE 
Date Time (U.T.) Mean period 

March 24 0610 90.03 
0440 89.95 

25 1510 89.92 

26 0440 89.88 

0305 89.77 

27 0435 89.76 

28 0300 89.68 

29 0255 89.57 

30 0245 89.46 

31 0235 89.31 

April 1 1415 89.02 
3 0150 88.69 

4 0125 88.47 

5 0100 88.15 


The estimated accuracy the time measurements Table seconds for and 
seconds for 
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Apparent Solar Periodicity Radio Star Scintillation 


From the beginning December 1957 the end November 1958, the radio sources 
Cassiopeia-A, Cygnus-A, Taurus-A, and Virgo-A were observed upper culmination from 
the Radio Observatory Queen’s University, Kingston, Ontario. The scintillation charac- 
teristics each the sources were studied Mc/s and significant correlations were 
obtained with local magnetic indices. 

generally accepted that the geomagnetic field exhibits cyclic changes which are related 
the period rotation the photosphere (Chapman 1940). The high correlation between 
scintillation and magnetic activity obtained the authors and others (Dagg 1957) would 
suggest similar periodicity scintillation data. attempt discover such periodicity 
exists the present records, autocorrelation function was calculated for the scintillation 
index for Cygnus-A, the source showing highest correlation with magnetic activity. The 
scintillation index used was measure both amplitude and rate scintillation. 

This function shown the accompanying figure, for intervals sidereal days. 
The limits shown are probable errors based solely the number observations. The solid 


a 
VU 
= 
z 
w 
a 
U 


SIDERIAL DAYS 


Fic. 1. 


line 5-day running mean indicate the trend the data. can seen that significant 
correlation obtained after about days. quasi-period 27.6 sidereal days has been 
estimated from the running mean, but this subject considerable error. 

These results strongly suggest that solar activity contributes directly the ionospheric 
processes responsible for the scintillation radio stars, least locations near the 
auroral zone. 


CHAPMAN, S. 1940. Geomagnetism, Vol. I (Clarendon Press, Oxford), pp. 370-384. 
Dacc, M. 1957. J. Atmospheric and Terrest. Phys. 10, 194. 
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Observations Macromosaic Substructures Lead 


Previous investigations macromosaic substructures called (Teghtsoonian 
and Chalmers 1951, 1952; Atwater and Chalmers 1957) have illustrated the effects freezing 
velocity, impurity type and content, and crystal orientation upon the formation and charac- 
teristics these dislocation arrays. The present investigation illustrates the effects produced 
the temperature distribution the liquid and solid the interface. 

Single crystals lead and lead-silver alloys the range at.% were grown 
horizontally graphite-coated lavite boats the zone-melting technique illustrated 
Fig. The unique feature this study that the crystals were grown under varying ambient 
temperatures illustrated. varying the ambient temperature, relative the melting 
temperature, the effect the temperature distribution array formation could 
studied. 

Two sizes crystals were grown, each in. wide and in. long, some in. deep and others 
in. deep. The thicker crystals were solidified with the interface perpendicular 
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Temperature Distribution 


SY 
Lavite Boat Liquid Zone 


Vycor Tube Zone Heater 
Furnace 


Specimen Arrangement 


Fic. 1. Zone melting apparatus and temperature distribution throughout specimen. 


the crystal axis, whereas was not always possible maintain this interface orientation 
for the thinner crystals the highest ambient temperature. 

The main body experiments were carried out freezing velocity, mm/min and 
the results are presented Table each ambient temperature two crystals were grown. 


TABLE I 
Striation formation 
Conditions Purity prwegel concentration at.% % Ag 

Note: — indicates no striations visible; + indicates striations visible. 


*This lead was zone-refined from a starting material of 99.9999% quoted purity. 


additional experiment was performed with and 100°C using 
various Two crystals were grown atv mm/min yielding visible arrays but much narrower 
striations than 2mm/min. Two crystals were grown 20mm/min and 
striation arrays were evident even though well-developed cellular substructure (Rutter 
and Chalmers 1953) was visible. 


None 

c 

c 

@ . . 

Striations 

a 

None 


Fic. 2. Schematic illustration of the dependence of striation formation on the solute concentration, Co, and 
the ratio of temperature gradient, G, to growth velocity, v. 
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Since decreasing decreases the temperature gradient the solid the interface and 
consequently the gradient the liquid the interface, the results Table can 
illustrated schematically Fig. With very pure material striations are produced 
even with low whereas, for crystals containing some measurable solute concentration, 
striation formation inhibited only under extremes high low However, was 
noted under such extremes well-developed cellular substructure was observed the 
crystals. 

was also observed that (i) the incubation distance for striations, the amount 
crystal grown before striation arrays appear, increased decreased, the width the 
striations decreased G/v decreased, and (iii) the misorientation decreased G/v decreased. 

These results suggest that, high purity lead, insufficient dislocations were grown into 
the crystals cause array formation and lead-silver alloys, the cellular substructure 
well developed, the solute segregation decreases the mobility dislocations the point 
that array formation cannot occur. 


Atwater, H. A. and CHALMERS, B. 1957. Can. J. Phys. 35, 208. 
Rutter, J. W. and CHALMERS, B. 1953. Can. J. Phys. 31, 15. 
TEGHTSOONIAN, E. and CHALMERS, B. 1951. Can. J. Phys. 29, 370. 
—— 1952. Can. J. Phys. 30, 388. 
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